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PREFACT 


This  project  on  sonic  rock  cutting  was  sponsored  by  the 
Advanced  Research  Projects  Agency  and  was  monitored  by  the  U.  S. 
Bureau  of  Mines,  Twin  Cities,  Minnesota  under  Contract  No. 
H0220037.  James  T.  Olson  of  the  Twin  Cities  Mining  Research 
Center  was  the  ARPA  Program  Coordinator,  with  Dr.  Patrick  J.  Cain 
of  the  Twin  Cities  organization  acting  as  Project  Officer.  Tl.e 
work  was  an  extension  of  earlier  work  under  Contract  No.  H0210010 
and  was  initiated  on  February  23,  1972  and  completed  on  April  22, 


The  work  was  performed  in  the  Engineering  Mechanics  Depart¬ 
ment  and  Sonic  Pow^r  Laboratory  of  The  Ohio  State  University 
under  Research  Foundation  Project  RF  3396.  Dr.  Karl  F.  Graff 
Piofessor  of  Engineerii  g  Mechanics  was  project  supervisor. 

Charles  C.  Libby,  Assistant  Professor  of  Welding  Engineering 
was  in  charge  of  the  layer  cutting  phase  of  work  from  February  23 
through  September  30,  1972.  Peter  Au  and  Gary  Streby  worked  on 
the  development  of  layer  cutting  techniques.  Long-chin  Shieh 
carried  out  experimental  work  in  rock  drilling  while  Ting-yu  Lo 
conducted  the  analysis  of  the  drilling  process.  The  analysis  of 
sonic  transducers  was  done  by  Ching-chio  Feng,  as  was  the  analy¬ 
sis  of  layer  cutting.  The  report  on  a  pulse-reflector  for  sonic 
transducers  was  prepared  by  Nelson  Ma.  Henry  A.  Bobulski,  Sonic 
Laboratory  Manager,  carried  out  instrumentation  and  assisted  in 
■the  development  phases  of  all  laboratory  work. 
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SUMMARY 


The  utilization  of  high  power,  high  frequency  sonic  transducers 
for  rock  cutting  has  been  in  progress  at  Ohio  State  University  for  two 
years  with  the  second  year  of  work  presented  in  this  report.  Two  modes 
of  rock  cutting  were  investigated;  the  first,  simple  drilling  and  the 
second,  layer  cutting.  An  impact  coupling  technique,  which  utilized  a 
small,  bouncing  cutting  tool  at  the  tip  of  the  transducer,  was  used  in 
both  modes  of  cutting.  The  overall  study  was  divided  into  a  system 
development  phase,  where  laboratory  cutting  apparatus  was  tested,  and 
a  system  analysis  phase,  where  basic  studies  of  sonic  processes  and 
devices  were  carried  out. 

Considering  first  the  analysis  phase  of  work,  basic  studies  of  the 
transducer,  the  drilling  process,  layer  cutting  and  of  a  pulse  reflector 
were  carried  out.  In  the  development  phase,  the  work  on  drilling  was 
continued  from  the  first  year  into  new  tool  designs  and  hard  rock 
cutting.  In  layer  cutting,  which  was  new  to  the  second  year,  apparatus 
was  constructed,  tested  on  limestone,  and  then  applied  to  cutting 
granite. 

While  drilling  and  layer  cutting  in  Indiana  limestone  was  success¬ 
ful  ir  the  sense  of  simply  being  able  to  remove  rock,  difficulties  were 
encountered  in  penetrating  hq^d  rock.  Numerous  variations  of  cutting 
parar.rters  were  tested,  until  some  success  in  layer  cutting  of  granite 
was  achieved.  Nevertheless,  results  were  marginal.  While  numerous 
cccplex  explanations  ure  available  regarding  the  difficulties,  it  is 
possible  to  sumr arize  them  quite  simply.  The  impact  coupling  technique 
of  energy  transfer  presently  represents  a  poor  impedance  match  between 
the  transducer  and  the  load. 
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The  possibility  of  using  ultrasonics  for  rock  cutting  was  noted  by 
Mauer  [l].*  The  technique  proposed  earlier  was  essentially  an  extension 
of  the  ultrasonic  slurry  drilling  technique  that  has  been  successfully 
used  for  many  y.-ars  for  cutting  extremely  hard  materials,  such  as 
tungsten  carbide.  However,  the  relative  inefficiency  and  slow  cutting 
rate  of  the  slurry  drilling  process  seemed  to  preclude  use  of  ultra¬ 
sonics  for  cutting  rock. 

However,  as  has  been  noted  in  an  earlier  report  [2],  developments 
in  power  sanies  at  Ohio  State  University  indicated  that  sonic  rock 
cutting  on  a  large  scale  might  be  feasible.  An  impact  coupling  means 
of  energy  transfer  had  been  developed  that  enabled  significant  amounts 
of  energy  to  be  transferred  to  a  work  surface.  Large  sonic  transducers, 
rated  at  10  kilowatts  had  been  developed.  Crude  experiments  on  drilling 
and  cutting  concrete  had  been  carried  out. 

Accordingly,  a  two-year  study  of  sonic  rock  cutting  has  been  con¬ 
ducted  of  which  this  report  covers  the  second  yen.  The  work  has  been 
divided  into  systems  development  and  systems  analysis  phases.  In  the 
development  aspect,  laboratory  rock  cutting  hardware  has  been  constructed 
and  tested.  In  the  analysis  aspects,  various  fundamental  characteristics 
of  sonic  devices  and  processes  have  been  investigated,  utilizing  both 
theory  and  experiment. 

The  work  in  the  second  year,  reported  herein,  was  aimed  at  utilizing 
two  modes  of  cutting  to  penetrate  hard  rock.  The  f  xrst  mode  was  the 
simple  drilling  of  holes  of  one  to  orie-and-a-quarter  inch  in  diameter. 

The  second  mode  was  layer  cutting.  In  this  work,  a  sonic  transducer 
was  traversed  across  the  rock  surface.  The  transducer  and  impact  tool 
were  so  arranged  as  to  cut  a  gouge  in  the  surface  of  one-quarter  to 
one-half  inch  in  depth.  Such  an  arrangement,  it  was  felt,  could  be 
extended  to  a  multiple  array  of  transducers  on  a  tunneling  machine. 

Fundamental  studies  tha-  accompanied  hardware  development  included 
analysis  of  the  transducer,  analysis  of  the  drilling  process  and  develop¬ 
ment  of  a  stress  pulse  reflector  for  the  transducers.  The  work  on  the 
transducer,  which  included  both  th' ory  and  experiment,  was  directed  at 
under  s  t  an  ding  the  influence  of  prectress,  voltage  and  temperature  on 
transducer  frequency.  The  drilling  analysis  also  comprised  theory  ana 
experiment  and  was  aimed  at  predicting  drilling  rate,  with  effects  of 
tool,  static  forces  and  rock  type  included  in  the  analysis.  The  pulse 
reflector  development  arose  from  the  need  to  mitigate  the  stress  pulses 
transmitted  into  the  transducer  during  impact  coupling. 


♦Numbers  in  brackets  refer  to  references  at  the  end  of  this  section. 


A  major  objective  of  the  second  year  of  work  was  cutting  of  hard 
rock,  and  as  basalt,  quartzite  and  granite.  Previous  work  had  been 
done  mainly  on  Indiana  limestone.  In  the  development  of  the  layer 
cutting  technique,  which  was  completely  new  in  the  second  year,  cutting 
tests  wore  carried  out  on  limestone  until  the  apparatus  was  perfected. 

In  terms  of  simple  ability  to  cut  rock,  results  were  acceptable  while 
working  in  the  softer  rock.  Considerable  difficulty  was  encountered 
in  cutting  hard  rock,  such  as  granite,  however.  This  also  held  true 
in  the  drilling  work.  While  some  success  was  achieved  in  cutting  hard 
rock,  the  results  were  quite  marginal,  both  in  terms  of  rate  of  rock 
removal  and  in  terms  of  specific  energy  of  removrl. 

Wh^le  numerous  complex  explanations  could  be  put  fo  *th  as  to  the 
diffioJLty  in  cutting  hard  rock  in  a  significant  way,  it  is  possible  to 
summarize  the  problem  quite  simply.  It  was  one  of  mechanical  (impedance) 
mismatch  between  the  transducer  and  the  load.  The  transducer  is  a 
powerful  motor,  known  to  be  capable  of  transmitting  10  kilowatts  of 
power  into  a  matched  load.  However,  the  impact  coupling  technique  is 
a  very  poor  matching  system  between  the  transducer  and  load.  The 
bouncing  tool  was  found  to  transmit  only  a  small  fraction  of  the  energy 
in  the  transdicer  to  the  rock  in  an  impact.  Further,  the  rate  of  tool 
impact  in  rock  cutting  is  on  the  order  of  hundreds  of  impacts  per 
second,  versus  transducer  vibration  frequencies  of  10, OCX)  Hz.  Thu , 
the  transducer  was  operating  effectively  unloaded  most  of  the  time. 

The  energy  transfer  process  was  thus  found  to  be  both  inadequate  and 
inefficient. 


Refers  aces 

1.  Mauer,  W.  C.,  Novel  Drilling  Techniques,  Pergamon  Press,  1968, 
pp.  44-49. 


2.  Graff,  K.  F.,  "Fundamental  Studies  in  the  Use  of  Sonic  Power  for 
R> ok  Cutting."  Annual  Technical  Report,  Bureau  oi  Mines  Contract 
No.  H0210010,  Dec.,  1971. 
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1.  SONIC  LAYER  CUTTING  OF  ROCK 


1-1  INTRODUCTION 

The  use  of  sonic  energy  for  cutting  layers  or  kerfs  in  reck  is  not 
new.  In  earlier  work,*  the  model  P-7  transducer  was  used  to  cut  kerfs 
in  concrete,  as  shown  in  Fig.  1-1.  In  this  illustration  the  transducer 
delivers  reciprocating  force  to  the  pivoted  blade  by  means  of  inter¬ 
mittent  or  bouncing  contact  at  the  upper  end  of  the  blade  or  tool. 

The  lower  end  of  the  tool  cuts  a  groove  in  the  limestone  rock  as  the 
transducer-tool  combination  is  traversed  from  right  to  left.  In  this 
assembly  the  traverse  is  made  parallel  to  the  rock  face  and  its  velocity 
is  determined  by  the  rate  at  which  a  control  wheel  is  rotated  manually. 

Rock  grooving,  by  this  means,  requires  relatively  large  amounts  of 
energy  per  cubic  volume  removed.  This  is  because  the  rock  is  highly 
pulverized  by  the  particular  tool  configuration  used.  Very  few  chips 
of  any  significant  size  are  produced.  The  groove  is,  therefore,  smooth¬ 
sided  and  sharp-edged.  The  process  thus  met  the  requirements  for  kerf 
cutting  in  highways  of  minimum  chipping  and  a  sharp-edged  groove. 
However,  in  layer  cutting  for  excavation  purposes,  maximum  volume  of 
rock  removal  and  rapid  cutting  speed  is  of  greater  interest.  In  par¬ 
ticular,  it  becomes  desirable  to  produce  considerable  chipping  from  the 
region  of  the  cut. 

The  basic  concept  of  layer  cutting  in  the  present  work  was  somewhat 
different  than  in  Fig.  1-1.  Basically,  the  approach  was  to  traverse  a 
sonic  transducer  across  the  face  of  the  rock,  using  a  small  impact  tool 
at  the  tip  of  the  transducer  instead  of  the  lengthy,  massive  tool  of 
Fig.  1-1.  By  slightly  inclining  the  transducer  relative  to  the  rock 
face,  it  was  possible  to  apply  sonic  impact,  via  the  impact  tool, 
directly  to  the  rock,  instead  of  indirectly  by  the  lever  action  of 
Fig.  1-1. 


The  use  of  a  small  impact  tool  and  the  impact  coupling  means  of 
energy  transfer  is  common  to  many  sonic  processes.  Basically,  in  the 
intermittent  or  noncontinuous  type  of  impact  mode  rock  cutting,  mechani¬ 
cal  energy  is  withdrawn  from  the  transducer  in  discrete  quantities 
during  the  intervals  between  impacts  of  the  metallic  tool  by  the  elec¬ 
trical  power  supply  system.  The  duration  of  tool  impact  is  very  short, 
a  small  percent  of  the  relatively  longer  time  between  impacts;  impact 
duration  is  fairly  constant,  whereas  the  duration  of  the  time  interval 
between  impacts  varies  inversely  with  the  amplitude  of  the  externally 
applieu  static  force.  This  static  force  is  applied  in  the  direction  of 
the  traverse  of  the  rock  cutting  assembly. 


*Work  sponsored  by  The  Ohio  Department  of  Highways  and  the  Federal 
Highway  Administration  under  Ohio  State  Engineering  Experiment  Station 
Project  EES  220,  "Development  of  Sonic  and  Ultrasonic  Power  Devices 
for  Application  in  Highway  Engineering." 
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In  the  amplication  of  sonic  energy  to  layer  cutting,  three  phases 
of  development  occurred.  In  the  first,  rather  crude,  manually-operated 
apparatus  was  used  to  study  various  principles  of  cutting.  In  the 
second  phase,  a  moderately  elaborate  piece  of  apparatus  was  developed 
to  carry  out  extensive  layer  cutting  tests  in  Indiana  limestone.  In 
the  final  phase  of  work,  layer  cutting  of  hard  rock  was  studied. 


1-2  PRELIMINARY  APPARATUS 

The  preliminary  apparatus  on  which  tests  were  made  on  the  principles 
of  layer  cutting  of  rock  by  sanic  impact  means  is  shown  in  Fig.  1-2. 

In  this  view,  a  surplus  machine  shop  grinder  has  been  utilized  to 
support  and  traverse  the  rock  relative  to  the  Model  P-11  transducer* 
and  its  tool  assembly.  The  block  of  limestone  rock  illustrated  hare  is 
bolted  down  to  the  work  surface  of  the  grinder.  The  handwheels  shewn 
in  the  lower  portion  of  the  pictures  are  utilized  to  raise  and  lower 
the  transducer  assembly,  and  to  traverse  the  assembly  horizontally. 

The  smal 3  impact  tool  at  the  tip  of  the  transducer  is  evident.  A 
second  tool  configuration,  detached  from  its  holding  fixture,  is  also 
being  shown  by  the  researcher  in  Fig.  1-2.  The  impact  tool  is  supported 
by  a  simple  pivot  about  two  inches  above  the  impact  point.  A  cantilever, 
attached  at  the  nodal  location  of  the  transducer,  holds  the  pivot 
assembly  fixed  relative  to  the  transducer.  A  closer  view  of  the  cutting 
apparatus  is  shown  in  Fig.  1-3,  where  a  partially  complete d  cut  has  been 
made.  Note  in  the  figures  that  the  transducer  is  slightly  inclined 
relative  to  the  rock  surface  and  that  the  line  of  transducer  action  is 
nearly  direct  into  the  cutting  surface. 

The  tool  is  free  to  swing  in  an  arc ,  with  its  motion  limited  at 
one  extreme  by  a  flat-faced  contact  with  the  transducer  tip,  and  at  the 
other  extreme  by  impact  with  the  rock.  As  the  transducer  is  moved 
parallel  to  the  rock  face,  the  tool  is  pushed  against  the  tip  of  the 
then  excited  transducer.  The  longitudinal  vibration  of  the  trfnsducer 
tip  applies  a  force-impulse  to  the  rool  during  a  short  interval  of  con¬ 
tact  time.  This  force-impulse  sends  the  tool  away  from  the  transducer, 
moving  in  an  arc  about  its  pivot,  and  into  impact  with  the  rock  surface. 
The  kinetic  energy  of  the  tool  is  converted  into  work,  breaking  the  rock 
surface.  After  the  tool  strikes  the  rock  face,  it  rebounds  toward  the 
transducer  tip,  which  is,  of  course,  also  moving  slowly  toward  the  tool 
due  to  the  traverse  motion  of  the  assembly.  The  process  of  impact  with 
the  transducer  and  the  rock  is,  therefore,  repeated  rapidly  by  the  tool, 
and  a  stable  impact- frequency  develops. 


*The  P-U  transducer  has  been  described  in  the  final  report  on 
Contract  H0210010. 
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Several  basic  parameters  controlling  the  layer  cutting  process 
were  identified  during  work,  with  the  preliminary  apparatus.  The  mass 
and  geometry  of  the  tools  were,  of  course,  parameters  of  the  process, 
as  they  are  for  drilling  or  other  impact  coupling  processes.  The  depth 
of  cut  and  speed  of  traverse  were  parameters,  as  was  the  angle  of 
inclination  of  the  transducer  relative  to  the  rock  face.  The  voltage 
to  the  transducer  was  a  parameter.  Finally,  the  material  properties  of 
the  rock  directly  influenced  the  cutting  process. 


1-3  iayer  CUTTING  APPARATUS 

The  design  and  assembly  of  a  larger  device  for  layer  cutting  of 
rock  proceeded  in-  °dir  ,ely  after  tests  in  the  preliminary  apparatus 
were  completed,  ine  resulting  apparatus  consists  essentially  of  a 
large,  welded,  steel  frame  with  a  track  of  approximately  eight  feet  in 
length,  along  which  a  carriage  may  be  traversed.  The  carriage,  on 
which  is  mounted  the  sonic  transducer  and  cutting  tool,  is  driven  by  a 
lead  screw  which  is  in  turn  chain  driven  by  an  adjustable -speed  gear 
motor . 

Figure  1-4  shows  an  overall  view  of  the  equipment,  with  frame, 
carriage,  transducer,  screw-drive  and  motor  all  shown.  Means  for 
supporting  and  clamping  rock  specimens  also  is  shown.  The  welded  steel 
frame  is  shewn  schematically  in  Fig.  1-5.  A  detail  of  the  electric 
motor  and  chain  drive  is  shown  in  Fig.  1-6.  By  varying  the  speed 
setting  of  the  motor,  as  well  as  by  interchanging  gears,  carriage 
traverse  speeds  from  0.25  to  4.75  inches  per  second  could  be  attained. 
The  speed  calibration  chart  for  the  system  is  given  in  Fig.  1-7. 

The  mounting  that  holds  the  transducer  at  the  nodal  point  and  in 
turn  attaches  to  the  carriage  is  designed  to  allow  the  transducer  to 
be  aligned  in  any  position  from  vertical  to  horizontal,  in  10°  incre¬ 
ments.  The  angle  adjustment  points  show  on  the  carriage  in  Fig.  1-4 
as  the  series  of  holes  along  a  circular  arc.  Provisions  for  load  cells 
were  also  incorporated  in  the  carriage  and  show  as  the  four  short, 
vertical  elements  between  the  wheeled  portion  and  transducer  holding 
portion  of  the  carriage  in  Fig.  1-4.  However,  actual  force  readings 
were  not  done  during  the  project. 


1-4  LAYER  CUTTING  IN  INDIANA  LIMESTONE 

In  carrying  out  tests  of  layer  cutting  in  Indiana  limestone,  the 
previously  mentioned  parameters  of  tool  mass  and  geometry,  depth  of 
cut,  speed  of  traverse  and  transducer  voltage  and  angle  of  attack  were 
varied.  Same  of  the  dependent  variables  of  the  process  included  power 
input  to  the  transducer  and  volume  of  rock  removed,  as  well  as  such 
items  as  frequency  of  tool  impact  and  forces  developed  in  the  supporting 
structure. 
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Apparatus  for  sonic  laye 


In  carrying  out  typical  layer  cutting  tests,  attention  was  given 
to  the  following  items : 

(1)  Transducer  Resonant  Frequency.  In  advance  of  any  series  of 
test  runs  utilizing  a  fixed-frequency  power  supply  with  a  sonic  trans¬ 
ducer,  it  is  desirable  to  first  establish  the  initial  resonant  frequency 
of  the  transducer.  In  this  work,  the  initial  resonant  frequency  of  the 
transducer  was  measured  at  no-load  on  the  Dranetz  impedance -admittance 
meter  at  the  point  of  minimum  impedance  while  a  variable  frequency, 

3  volts  excitation  level  is  applied  to  the  transducer.  In  actual  oper¬ 
ation  the  no-load  resonant  frequency  of  the  P-11  transducer  is  found  to 
vary  inversely  with  the  applied  terminal  voltage  as  previously  described.* 

Consequently,  it  was  necessary  to  time  the  transducer  to  a  higher 
frequency  than  the  power  supply  frequency,  to  account  for  the  frequency 
drop  due  to  high  voltage.  In  a  typical  test  run,  the  P-11  transducer 
had  an  initial  no-load  resonant  frequency  of  10,l40  Hz,  the  power 
supply  delivered  a  fixed  frequency  of  10,054  Hz  and  the  excitation 
level  measured  at  the  transducer  terminals  under  load  was  about  1000 
volts. 

(2)  Speed  of  Travel.  In  the  tests  using  the  manual  traverse 
mechanism,  it  was  not  possible  to  measure  speed  of  traverse  directly. 
Measurements  were,  therefore,  made  of  the  time  of  traverse  from  one 
face  of  the  limestone  block  to  the  other  face.  On  the  other  hand, 
tests  using  the  large  console  apparatus  were  made  with  predetermined 
speed  settings.  However,  no  estimate  of  the  possible  speed  reduction 
of  the  motor  due  to  load  was  included. 

(3)  Depth  of  Cut.  In  initially  establishing  the  depth  of  cut  to 
be  made,  measurements  were  made  of  the  interference  between  the  tool 
and  the  rock  face.  After  the  cut  was  made,  further  measurements  were 
taken  of  the  depth  to  which  rock  was  actually  removed.  The  difference 
between  these  measurements  gave  an  indication  of  the  relative  stiffness 
of  the  structure.  The  depth  measurements  of  the  cut  in  the  rock  were 
recorded  for  c; Isolating  the  volume  of  rock  removed. 

(4)  Indication  of  Overload.  It  was  not  possible  during  these 
actual  test  runs  to  monitor  the  frequency  of  impact  or  to  measure  the 
static  or  average  forces  required  to  drive  the  sonically  excited  tool 
across  the  face  of  the  limestone  rock.  As  mentioned  in  section  1-1 
the  frequency  of  impact  is  inversely  proportional  to  the  static  or 
average  force  applied  to  draw  the  transducer  and  its  tool  across  the 
rock  face.  A  very  high  frequency  of  impact  may  indicate  the  process 
is  fully  loaded. 


♦Annual  Technical  Report  on  previous  Contract  H0210010,  "Fundamental 
Studies  in  the  Use  of  Sonic  Power  for  Rock  Cutting,"  December,  1971- 


13 


In  sonic  impact  rock  cutting,  an  excessive  traverse  velocity  or  an 
excessive  depth  of  cut  may  be  defined  as  one  in  which  the  transducer 
cannot  maintain  impact-coupling  with  the  tool,  but  instead  becomes 
coupled  solidly  and  continuously  with  it.  Therefore,  the  point  at  which 
the  transducer  fails  to  maintain  the  impact  mode  of  operation  is  noted 
during  each  test  run  as  a  part  of  the  test  data.  This  point  corresponds 
to,  and  is  analogous  to,  a  point  in  a  motor  driven  process  at  which  the 
drive  motor  stalls.  To  all  intents  and  purposes,  tbn  process  under 
investigation  ceases  to  function  at  this  juncture. 

(5)  Volume  of  Rock  Removed.  Tne  volume  of  rock  removed  was  based 
on  the  dimensions  of  the  groove  formed  in  a  flat  surface  of  the  rock. 

The  cross  section  of  this  groove  is  generally  trapezoidal  in  shape. 

The  width  of  the  narrow  base  of  the  trapezoid  (w)  corresponds  to  the 
width  of  the  impact  tool.  The  height  of  the  trapezoid  (h)  corresponds 
to  the  depth  of  the  cut  in  the  rock  by  the  impact  tool.  The  width  of 
the  cut  at  the  surface  of  the  rock  (v)  varies  rather  widely  where  large 
chips  are  removea.  It  is  in  the  need  to  estimate  dimension  v  which 
results  in  a  reduction  of  the  accuracy  of  volume  measurement  to  ±  5$. 

The  area  (A)  of  the  trapezoidal  section  of  the  cut  is  calculated 
from  these  dimensions;  A  equals  one-half  the  sum  (w  +  v)  multiplied  by 
the  height  (h).  The  volume  of  rock  removed  is  then  calculated  by 
multiplying  the  area  A  by  the  effective  length  of  the  cut  (z).  The 
length  of  the  cut  must  be  estimated  due  to  edge  effects.  The  impact 
tool  tends  to  break  out  fairly  large  portions  of  the  rock,  both  at  the 
beginning  and  at  the  end  of  a  block  it  is  grooving.  In  estimating 
volume  of  rock  removed,  the  break-out  volume  has  been  ignored. 

The  accuracy  of  volume  estimates,  based  on  these  measurements  has 
been  found  to  be  within  %  of  the  actual  volume  of  rock  removed, 
measured  by  displacement  means. 

(6)  Angle  of  Attack  and  Face  Angle.  The  angle  of  attack  of  the 
transducer  may  be  defined  as  the  inclination  of  the  transducer's  center- 
line  in  reference  to  the  horizontal  line  of  the  work  surface.  By  in¬ 
creasing  this  angle  it  is  possible  to  increase  the  downward  component 
of  externally  applied  static  force.  The  increased  angle  of  attack  has 
been  assumed  to  be  required  for  harder  reck. 

It  must  be  recalled  that  the  impact  tool  position  is  fixed  to  the 
transducer-carriage  assembly.  The  angle  between  the  face  of  the  cutting 
tool  and  the  rock  is  theref  re  not  only  dependent  on  the  tool  design, 
but  on  the  angle  of  attack  of  the  overall  transducer.  Thus,  if  the 
angle  of  attack  is  varied,  generally  the  tool  must  be  redesigned  to 
account  for  the  variation  in  face  angle  relative  to  the  rock. 


1-5  RESULTS  FOR  LIMESTONE 


Various  tool  designs  were  utilized  in  the  limestone  cutting  work. 

A  typical  impact  tool  is  detailed  in  Fig.  1-8.  The  drawing  illustrates 
the  70°  face  angle  of  the  tool  and  the  carbide  insert  in  uhe  flat  cutting 
face  of  the  tool.  The  cutting  face  of  this  3/4"  wide  tool  is  held 
perpendicular  to  a  horizontal  work  surface  when  the  transducer  is 
inclined  at  20°  frcm  the  horizontal.  With  the  transducer  at  a  30°  tilt, 
the  tool  face  is  inclined  at  about  80°  from  the  horizontal  work  surface. 
Other  tool  designs  can  be  used  to  modify  the  relationship  between  the 
transducer  position,  or  its  attack  angle,  and  the  angle  which  the  tool 
face  makes  with  the  work  surface. 

Figure  1-9  illustrates  a  second  generation  tool  design  in  which 
the  impact  face  of  the  tool  has  been  modified  into  a  plow-shape  with  a 
130°  included  angle  between  the  faces  of  the  tool.  The  cutting  angle 
of  the  tool  is  Approximately  the  same  as  in  the  flat-faced  tool  in 
Fig.  1-8,  although  the  tool  width  is  reduced  to  l/2".  This  narrower 
tool  should  allow  faster  cutting  speeds  with  the  same  rock,  or  perhaps 
cutting  of  harder  rock  at  the  same  traverse  speed  as  the  3/4"  flat¬ 
faced  tool.  The  cutting  face  of  this  tool  is  held  perpendicular  to  the 
work  surface  with  the  transducer  inclined  at  15°  fran  the  horizontal. 

(A)  Results  from  Preliminary  Apparatus 

Figure  1-10  illustrates  the  first  data  obtained  on  layer 
cutting  of  limestone.  The  results  were  obtained  using  the  manually 
traversed  apparatus  of  Fig.  1-2.  The  data  presents  the  rate  of  rock 
removal  (in.3/sec)  versus  cutting  depth.  Because  a  manual  traverse 
was  used,  velocity  determinations  were  only  approximate. 

Figure  1-11  presents  the  preceding  data  in  a  different  form, 
showing  the  specific  energy  of  rock  removal  versus  cutting  '•'"pth.  The 
major  point  of  interest  is  the  existence  of  a  pronounced  minimum  at 
about  0.5-inch  cutting  depth. 

(b)  Results  from  Large  Cutting  Apparatus 

An  illustration  of  the  cutting  action  in  limestone  during  a 
typical  test  is  shown  in  Fig.  1-12.  Fairly  large  rock  fragments  are 
shown  being  dislodged.  A  large  breakout  area  at  the  end  of  the  block 
often  existed  and  was  discounted  in  all  data  reduction. 

Figures  1-13  and  1-14  present  data  on  volume  of  rock  removal 
versus  cutting  speed  for  two  different  depths  of  cut  (l/2"  and  l/4"). 

The  volume  of  excavation  shown  on  the  ordinate  of  Fig.  1-13  is  for  a 
6"  traverse  distance  in  limestone  with  a  3/4"  wide  tool.  The  angle  9 
is  the  attack  angle  of  th >  transducer,  that  is,  its  angle  of  in"1  ination 
from  the  horizontal.  Note  that  dotted  lines  indicate  the  point  at  which 
speed  of  travel  exceeded  the  level  at  which  impact  cutting  was  possible. 
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Fig.  1-10  Rate  of  rock  removal  for  early  apparatus 
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Fig.  1-12  Layer  cutting  action  on  large  apparatus 


Fig.  1-13  Volume  of  rock  removal  versus  cutting  speed 
for  a  0.  5  inch  depth  of  cut 
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Fig.  1-14  Voluire  of  rock  removal  versus  cutting  speed  for 
a  0.25  inch  depth  of  cut 


The  dotted  lines  indicate  possible  extrapolation  into  an  area  of  overload 
on  the  transducer.  Figure  1-13  indicates  that  the  impact  mode  of  oper¬ 
ation  could  not  be  maintained  with  a  l/2"  cut  depth,  and  a  20°  angle 
of  attack,  for  more  than  2.25  in. /sec  traverse  speed.  At  speeds  higher 
than  this,  the  tool  maintained  constant  contact  with  the  transducer, 
impact  cutting  ceased,  the  transducer  detuned,  its  power  input  dropped 
and  sonic  cutting  was  brought  to  a  halt . 

Figure  1-14  shows  that  with  the  same  tool  and  half  the  depth 
of  cut,  the  velocity  of  traverse  can  be  almost  doubled  when  compared  to 
Fig.  1-13.  The  point  at  which  impact  cutting  ceases  with  a  20°  attack 
angle  and  a  1/4"  depth  of  cut  is  4.75  in. /sec  compared  to  2.25  in. /sec 
with  the  1/2"  cut  depth.  It  is  to  be  noted  in  Figs.  1-13  and  l-l4  that 
the  volume  of  rock  removal  tends  xo  increase  for  a  shallower  angle  of 
attack  (exception:  the  9  =  10° ,  20°  curves  of  Fig.  1-14). 

In  Figs.  1-15  and  1-16,  data  on  the  specific  energy  of  rock 
removal  versus  cutting  speed  is  given.  The  data  is  for  1/2"  (Fig.  1-15) 
and  1/4"  (Fig.  1-16)  depths  of  cut  and  for  various  angles  of  attack. 

It  is  first  particularly  significant  to  note  the  rnnid  decrease  in 
specific  energy  requirements  for  increasing  cutting  speed  for  both 
depths  of  cut.  It  is  also  to  be  noted  that  the  specific  energy  of 
removal  is  rather  insensi  vve  to  the  attack  angle. 

It  is  of  particular  interest  that  the  specific  energy  curves 
are  "leveling  off"  at  the  higher  cutting  speeds .  Although  a  clear , 
lower  level  asympote  does  not  exist  .in  either  Fig.  1-15  or  1-16,  it 
does  not  appear  that  specific  energies  of  much  less  than  500  W-sec/in.3 
would  occur  for  the  l/2-inch  cutting  depth  or  less  than  1000  W-sec/in.3 
for  the  l/4-inch  cutting  depth. 

Figure  1-17  presents  some  of  the  preceding  data  in  terms  of 
specific  energy  for  various  traverse  speeds  for  a  fixed  angle  of  attack 
of  10°  .  The  specific  energy  is  thus  seen  to  decrease  with  increases  in 
both  traverse  speed  and  with  depth  of  cut.  Both  of  these  are  independent 
variables  and  each  tends  to  increase  the  loading  on  the  transducer. 

The  specific  energy,  therefore,  is  improved  as  the  transducer's  load 
increases . 


1-6  LAYER  CUTTING  IN  GRANITE 

The  ultimate  objective  of  the  work  on  rock  cutting  was  to  effec¬ 
tively  penetrate  hard  rock,  such  as  Barre  Granite,  Dresser  Basalt  and 
Sioux  Quartzite.  After  some  success  in  cutting  limestone,  the  effort 
was  made  to  transition  the  same  cutting  techniques  to  hard  rock.  All 
work  reported  in  the  following  was  carried  out  on  granite  specimens 
obtained  locally. 
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Fig.  1-15  Specific  energy  of  rock  removal  versus  cutting 
speed  for  0.5  inch  cutting  depth 
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Fig.  1-1/  Specific  energy  of  removal  versus  cutting  speed  for 
two  depths  of  cut  (attack  angle  at  10°) 


(A)  General  Difficulties  Encountered 

Difficulties  were  immediately  encountered  in  attempting  to 
cut  hard  rock;  the  tools  and  operating  conditions  sufficient  to  penetrate 
limestone  were  simply  not  adequate  for  cutting  granite.  It  was  found, 
for  example,  that  cutting  depths  of  l/2-inch  were  not  possible.  Nearly 
all  subsequent  work  was  restricted  to  l/l6-  to  3/l6-inch  attempted 
depths  of  cut.  Secondly,  the  previous  voltage  level  applied  to  the 
transducer  of  1200  volts  was  inadequate.  All  reported  tests,  therefore, 
were  run  at  2500  volts,  which  was  close  to  the  maximum  voltage  rating 
of  the  transducer.  Finally,  traverse  speeds  of  several  inches  per 
second,  possible  in  limestone,  were  not  possible  in  granite.  Speeds 
of  approximately  1  in. /sec  were  used  on  most  tests,  with  higher  speeds 
being  attained  only  near  the  end  of  the  research  work. 

(B)  Attack  Angle 

Shallow  attack  angles  of  10  to  20°  were  found  most  successful 
in  limestone.  However,  it  was  found  that  in  granite  shallow  attack 
angles  generally  did  not  yield  satisfactory  cutting.  What  usually 
occurred  was  that  the  cutting  tool  would  "ride  up"  on  the  rock.  That 
is,  instead  of  maintaining  a  given  depth  of  cut,  the  cut  depth  would 
diminish  until  the  tool  was  skidding  along  the  rock  surface. 

The  fact  that  the  tool  was  able  to  ride  up  was  a  result  of 
some  elasticity  and  backlash  in  the  system.  The  cantilever  tool 
support  was  the  major  source  of  elastic  deformation.  Some  slight 
looseness  between  the  wheels  of  the  carriage  and  the  track  was  the 
major  source  of  system  backlash. 

One  approach  to  the  ride  up  problem  was  merely  to  increase 
the  depth  of  cut  beyond  the  deformation -backlash  range  of  the  system. 

A  depth  of  l/2-inch  was  sufficient  for  this.  However,  as  previously 
indicated,  the  transducer  system  was  generally  incapable  of  penetrating 
the  granite  at  these  depths. 

The  approach  generally  used  was  to  increase  the  angle  of 
attack  from  30  to  40° .  This  increased  somewhat  the  downward  cutting 
action  of  the  tool  and  thus  tended  to  maintain  the  depth  of  cut. 

Probably  equally  important,  it  decreased  the  bending  on  the  cantilevered 
tool  support  arm  and  thus  reduced  the  elastic  deformation  of  the  system. 

Additionally,  stiffening  was  added  to  the  tool  support  arm.  Shims 
were  placed  between  the  carriage  wheels  and  the  track  to  reduce  system 
backlash.  Even  with  these  modifications  and  the  change  in  attack  angle, 
ride  up  of  the  tool  continued  to  occur  occasionally. 

A  last,  more  or  less  impromptu,  measure  taken  near  the  end  of 
work  was  to  manually  force  the  tool  downward  by  pressing  against  the 
cantilever  tool  support  fixture.  This,  in  conjunction  with  the  other 
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measures  taken  enabled  fairly  consistent  cuts  to  be  made  in  granite  near 
the  end  of  the  work  period. 

(c)  Tool  Designs 


Cutting  tool  widths,  3/4-inch  to  1  inch,  were  used  in  limestone. 
In  granite,  tool  widths  were  decreased  to  l/2-inch  to  increase  the  con¬ 
tact  stress  against  the  rock.  Initial  tests  were  run  using  the  plow¬ 
shaped  tool  previously  shown  in  Fig.  1-9.  With  this  tool,  it  was  found 
that  the  carbide  tip  was  repeatedly  fractured  during  cutting.  Sharp 
comers  were  broken  away.  Fractures  through  the  carbide  near  the  inter¬ 
face  between  steel  and  carbide  were  noted.  Finally,  fractures  along 
the  narrow  neck  of  the  tool  occurred.  Some  of  these  failures  are  shown 
in  Fig.  1-22. 

One  of  the  tool  redesigns  was  a  heavier,  plow-shaped  tool, 
shown  in  Fig.  1-18.  It  is  to  be  noted  that  the  neck  region  of  the  tool 
is  considerably  thickened  over  that  of  Fig.  1-9.  Additionally,  it  was 
felt  that  a  more  massive  tool  would  increase  the  energy  transfer  during 
the  impact  coupling  process. 

Another  design  was  the  saw-edged  or  stepped  tool,  shown  in 
Fig.  1-19.  The  idea  was  to  modify  the  plow  tool  to  obtain  a  stepped 
cutting  action,  with  the  cutting  steps  on  the  order  of  the  cutting 
depth  of  a  single  edged  tool.  Accordingly,  two  l/8-inch  steps  were 
placed  in  the  tool. 

A  chisel-shaped  tool,  shown  in  Fig.  1-20  was  also  tested. 

This  tool  used  a  slightly  more  complicated  pivoting  action  than  the 
simple  pivots  of  the  previous  tools.  In  addition,  it  had  an  even 
narrower  cutting  edge  (l/4").  During  test,  it  was  found  that  the 
looseness  of  the  pivot  arrangement  caused  the  tool  to  ride  up  quite 
easily. 


A  two-piece  round  impact  tool  was  designed,  as  shown  in 
Fig.  1-21.  The  cutting  tool  itself  had  a  rounded  transducer  impact 
point  to  improve  energy  transfer  to  the  tool  during  impact  coupling. 

The  cutting  edge  of  the  tool  had  no  particularly  new  design  features. 
The  second  part  of  the  tool  was  the  holder,  designed  to  pivot  in  the 
manner  of  earlier  tool  designs  and  to  rigidly  hold  the  tool.  The 
overall  intent  of  this  design  was  to  enable  a  number  of  easily  replace¬ 
able  tool  designs  to  be  tested.  In  operation,  it  was  generally  found 
that  the  tool  could  r.ot  be  rigidly  retained  in  the  holder  during  im}  act 
action. 


Several  of  the  tool  configurations  are  shown  in  Fig.  1-22.  Mainly, 
the  figure  depicts  some  of  the  fractures  and  failures  that  occurred  in 
the  tools.  Thus,  carbide  inserts  are  seen  to  be  fractured  or  missing 
in  two  of  the  tools.  A  fracture  across  the  neck  of  one  of  the  tools 
also  is  evident.  A  heavy  plow-shaped  tool  is  shown  as  is  the  stepped 
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Fig.  1-19  Stepped  cutting  tool 
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Fig.  1-21  Two-piece  impact  tool 


illiaairm 

iSjiiJUriMj 

Fig.  1-22  Various  layer  cutting  tools 


cutting  tool,  although  the  stepped  edges  are  not  revealed  clearly  in 
the  latter. 

(D)  Spring-Return  System 

An  important  step  of  the  impact -coupling  process  is  the  phase 
of  tool  rebound  from  the  load.  In  particular,  if  tool  rebound  is  too 
low,  energy  pickup  from  the  transducer  is  also  low,  leading  to  rather 
weak  impact  against  the  load.  It  was  felt  that  a  spring-assisted 
return  of  the  tool  against  the  transducer  would  enhance  the  cutting 
process.  A  spring  was  installed  in  front  of  the  tool  to  increase  tool 
rebound  velocity. 

Mechanical  springs  of  different  stiffhesses  were  tried. 
Eventually ,  a  small  block  of  silicone  rubber,  yielding  a  fairly  stiff 
spring,  wan  found  to  work  best.  The  configuration  of  spring  and  tool 
is  shown  in  Fig.  1-23. 
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(E)  Transmission  Line  Optimization 

Two  fixed  frequency  motor-generators  were  available  for 
supplying  high  frequency  power  to  the  transducer.  The  output  frequencies 
of  the  generators  were  slightly  different,  one  being  9>720  Hz,  the  other 
10,054  Hz.  For  optimum  action,  it  was  naturally  desirable  that  the 
transducer-transmission  line  resonance  be  as  close  as  possible  to 
generator  frequency  at  operating  voltages.  However,  the  operating  con¬ 
ditions  of  high  voltage  and  temperature  increases  tend  to  lower  trans¬ 
ducer  resonance  from  the  no-load,  low  voltage  conditions  under  which 
it  is  tuned.  Consequently,  it  was  felt  desirable  to  conduct  a  series 
of  tests  to  determine  th'.  optimum  transmission  line  length  to  achieve 
tuned  operation  under  load. 

The  tests  started  with  a  transducer  fitted  with  a  long  trans¬ 
mission  line,  so  that  the  system  resonant  frequency  was  considerably 
lower  than  the  supply  frequency.  A  layer  cutting  run  was  then  performed, 
with  the  transducer  input  power  being  noted.  The  transmission  line  wa.? 
then  trimmed  by  l/l6  inch  and  the  cutting  process  repeated.  As  the  line 
was  shortened,  the  power  input  increased,  np  to  a  point,  and  then  started 
to  diminish.  The  process  was  continued  until  it  was  obvious  that  the 
optimum  length  had  been  reached  and  passed. 

Figure  1-24  shows  the  low  voltage,  no-load  resonant  frequencies 
for  various  transmission  line  lengths.  The  slope  of  the  average  straight 
line  is  25  Hz  drop  per  l/l6  incii  increase  in  line  length.  In  Fig.  1-25, 
the  inprut  power  to  the  transducer  under  cutting  conditions  for  various 
transmission  line  lengths  is  shown.  Obviously,  a  very  pronounced  peak 
is  reached  at  a  length  of  10-1/16  inches.  Referring  to  Fig.  1-24,  it 
is  seen  that  this  length  has  a  no-load  resonant  frequency  of  10,125  Hz. 

1-7  RESULTS  FOR  GRANITE 

Although  many  difficulties  were  encountered  in  layer  cutting  of 
granite,  some  successful  cuts  were  accomplished.  Figure  1-26  shows  a 
cut  being  made  with  the  heavy,  plow-shaped  tool  of  Fig.  1-18,  with  the 
spring-assist  device  of  Fig.  1-23  also  being  utilized.  Figure  1-27 
gives  a  general  indication  of  the  depth  and  profiles  of  seme  of  the 
cuts . 

Figure  1-28  shows  data  obtained  on  the  volume  of  rock  removal  for 
various  cutting  speeds,  at  two  different  angles  of  attack.  The  specific 
energy  of  rock  removal  versus  cutting  speed  is  shown  in  Fig.  1-29  for 
the  two  different  angles  of  attack.  The  data  at  the  40°  angles  of 
attack  is  "well  behaved"  in  both  Fig.  1-28  and  1-29,  in  that  it  resem¬ 
bles  the  general  results  obtained  in  limestone  cutting.  The  results 
for  30° >  however,  indicate  a  marked  decrease  in  removal  efficiency  for 
speeds  beyond  about  l/2  inch  per  second.  This  result  mainly  reflects 
the  tendency  of  the  tool  to  ride  to  the  surface  at  lower  angles  of 
attack . 
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1-24  No-load,  low-voltage  resonant  frequency  change 
change  in  line  length 
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Fig.  1-25  Input  power  versus  line  length  under  cutting  conditions 
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Fig.  1-29  Specific  energy  of  rock  removal  versus  cutting 
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TECHNICAL  LITERATURE  APPLICABLE  TO 
ANALYSIS  OF  LAYER  CUTTING 


2-1  INTRODUCTION 

Sonic  layer  cutting  or  rock  utilized  a  transducer  vibrating  at 
high  frequency  (10  kHz  for  the  present  application)  as  the  energy 
source.  The  mechanical  vibrational  energy  is  transmitted  from  the 
transmission  line  of  the  transducer  to  the  rock  by  a  bouncing  tool  via 
an  impact  coupling  process  .  Although  significant  amounts  of  energy  are 
transmitted  to  the  work  piece  by  the  rapid  impacts,  the  energy  transfer 
involved  in  a  single  blow  is  relatively  small  compared  to,  say,  conven¬ 
tional  percussive  drilling.  In  order  to  induce  failure  of  the  rock 
through  many  small  increments  of  energy  is  a  different  process  than  a 
low  frequency,  high  energy  per  blow  process;  therefore,  it  is  most 
important  to  understand  the  effects  of  process  parameters  such  as  depth 
of  cut  and  tool  geometry  on  rock  fracture.  A  better  understanding  of 
the  rock  fracture  mechanics  was  required.  With  this  in  mind,  a  brief 
review  on  the  failure  theories  of  rock  is  given  before  proceeding  to 
the  technical  literature  applicable  to  the  layered  cutting  process 


2-2  FAILURE  THEORIES  OF  ROCK 

It  is  well  known  that  rocks  generally  are  highly  inhomogeneous, 
anisotropic  and  porous  .  The  constituents  and  structures  vary  from 
rork  to  rock.  It  is  indeed  very  difficult  to  give  a  universal  failure 
theory  to  cover  all  rocks.  Nevertheless,  some  o  nervations  have  been 
made  on  the  general  failure  modes  of  most  rocks .  Rocks  are  said  to  be 
brittle  under  low  confining  pressure.  A  number  of  failure  theories  have 
been  proposed  and  reviewed  in  several  papers  [l],  [21.*  The  various 
theories  can  be  divided  into  three  categories:  (A)  the  empirical 
criteria,  (b)  the  phenomenological  failure  theories,  and  (c)  the 
statistical  failure  theories.  In  the  following,  a  brief  review  will 
be  presented  on  the  above  three  categories. 

(A)  The  Empirical  Criteria 

This  category  includes  the  maximum  shear-  stress  theory  and 
the  Coulomb -Navier-Kohr  theory.  The  maximum  shear  theory  states  that 
failure  occurs  at  the  point  where  the  maximum  shear  stress  becomes 
equal  to  a  critical  value  called  the  shear  strength  of  the  material 
which  is  independent  of  the  stress  state.  If  oi  >  cr2  >  03  denote  the 
three  principal  stresses,  this  theory  predicts  that  failure  will  occur 
on  the  plane  inclined  to  the  maximum  stress  cq ,  and  whose  normal 
bisects  the  (oi,a3)  axis.  According  to  this  theory,  the  uniaxial 
tensile  and  compressive  strengths  shou'.d  be  the  same .  Ductile  materials 
(such  as  metals)  appear  to  be  the  only  materials  which  ^ail  in  this  way. 


*Numbers  in  brackets  refer  to  references  at  the  end  of  this  chapter. 
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However,  for  most  rock  materials,  the  strengths  are  usually  dependent 
on  the  stress  state  in  that  their  strength  increases  .with  an  increasing 
confining  pressure.  Moreover,  the  compressive  strength  of  a  rock  is 
usually  much  greater  than  its  tensile  strength.  The  failure  planes  are 
inclined  less  than  45°  (about  30° »  with  a  tendency  to  increase  slightly 
as  the  confining  pressure  increases)  to  the  direction  of  maximum  com¬ 
pression.  Therefore,  it  is  not  proper  to  apply  the  maximum  shear  theory 
to  rocks. 

To  overcome  some  of  the  deficiencies  of  the  maximum  shear  criterion, 
the  Coulomb -Navier  Criterion  was  proposed  by  Navier  (see,  for  example, 
King  [3])  by  utilizing  the  Coulomb  friction  concept  (see,  for  example, 
Nadai  [4]).  This  criterion  may  be  stated  as  follows:  Suppose  that 
0i  >  03  are  the  principal  stresses  (compressive)  in  two  dimensions. 

The  normal  compressive  stress  cr  and  the  shear  stress  a  on  a  plane 
inclined  to  <73  by  an  angle  0  (see  Fig.  2-1)  may  be  written  as 

0  =  \  (01  +  03)  +  |  (01  -  03)  cos  20 

*  -  !  (cr,.  -  03)  sin  29  (2-1) 


(o)  (b) 


Fig.  2-1  A  two-dimensional  stress  state 
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Assume  that  when  failure  occurs,  element  A  (Fig.  2-1  (b))  will  slip 
along  the  failure  plane  0  =  0O,  as  the  shear  stress  along  the  failure 
plane  exceeds  the  resisting  strength  T,  which  was  proposed  by  Navier 
to  be 


T  =  x0  +  pc 


(2-2) 


where  t0,  a  material  constant,  is  called  the  intrinsic  shear  strength 
of  the  material  and  4  is  the  internal  Coulomb  friction  coefficient. 
Physically,  4  is  the  coefficient  of  sliding  friction  of  one  portion  of 
the  rock  against  the  other.  Incorporating  eq.  (2-2)  into  the  assumed 
failure  criterion,  the  failure  condition  is  obtained  as 


M  > 


+  (!CT 


(2-3) 


(2-4) 


Thus,  failure  occurs,  when 
To  =  M  -na 

=  “  |  +  cr3)  +  |  (ctx  -  or3)(sin  20o  -  M  cos  20o) 

Frcm  eq.  (2-3)  it  is  seen  that  for  a  given  t0,  oi  and  ct3,  the  potential 
failure  plane  can  be  decided  by  the  angle  0O  which  makes  the  right  hand 
side  of  eq.  (2-3)  a  maximum  or  it  is  the  angle  0O  which  makes  the  stresi 
state  satisfy  eq.  (2-3),  whichever  occurs  first.  This  0n  is  found  to 
be 


tan  20n  ~  — 

Substituting  eq.  (2-5)  into  eq.  (2-4)  gives 

°i  (  ?  +  1  -  4)  -  a3(  n/JP" TT  +  4)  =  2tq 

or. 


(2-5) 


(2-6) 


2t0  = 


tan120o  +  ^  +  -  CT3) 

This  may  be  written  as 

oi  =  2t0  tan  0O  +  c r3  tan2  0O 


2  0r 


(2-7) 


(2-8) 


From  eq.  (2-8),  the  relation  between  oi  and  ct3  is  seen  to  be  a  straight 
line  (see  Fig.  2-2).  Denoting  20o  =  n/2  +  <p,  then  sin  20o  =  cos  cp, 
cos  20o  =  -sin  cp  and  thus  4  =  tan  qp  =  -1/tan  20o,  so  that'  eq.  (2-4) 
becomes 
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illka 
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tan  a=  tan2 8Q- 


i  .  Sin 


A 

A 


I— Sin 


r0  ton  0O=  Ojo 


Fig.  2-2  The  (Xj.  -03  failure  boundary 


o_  _  1  -  sin  cp  1  +  sin 

tin  =  <Ti  - 31  -  O’ 3  - 

u  A  COS  cp  J  cos  cp 


(2-9) 


Tm  =  t0  cos  <p  +  am  sin  cp 


(2-10) 


where  Tm  =  \  (o’!  -  cr3)  is  the  maximum  shear  and  crm  =  i  (oq  +  ct3)  is  the 
mean  stress.  From  eq.  (2-5),  it  is  seen  that  0O  >  n/4  which  is  in 
accord  with  the  experimental  observations.  From  eq.  (2-6),  it  is  seen 
that  the  tensile  strength  cr^o  is  given  by  letting  q  =  0  so  that  the 
equation  becomes 


°to  -  ~&3  ~ 


'[l"  +  1  +  q 


(2-11) 


The  compressive  strength  crco  is  given  by  substituting  a3  =  0,  or 


°co  -  °i  = 


+  1  -  u 


(2-12) 


Therefore,  crto  4  crco,  and  their  ratio  is 


°~tn  _  VU  +  1  -  U 
aco  u 


(2-13) 


which  is  less  than  1  and,  thus,  is  in  accord  with  the  experimental 
observations.  Equation  (2-4)  can  also  be  expressed  as 


t0  =  t  |  -  tan  cp  •  a 


(2-110 


which  is  plotted  in  Fig.  2-3,  and  is  again  a  straight  line.  Note  that 
the  circles  tangent  to  the  straight  lines  are  the  Mohr's  circles  which 
give  the  stress  states  at  failure.  The  failure  line  is  also  called  the 
Mohr  envelope. 


Fig.  2-3  Mohr  envelope  of  failure 

Several  shortcomings  ail  so  exist  with  this  criterion,  such  as; 

'i)  the  theory  cannot  predict  the  value  of  the  principal  stresses  at 
i  'lure,  and  (ii)  u  is  difficult  to  obtain.  However,  the  failure 
envelope  does  have  fairly  good  agreement  with  the  observed  results  for 
rock-like  materials.  Moreover,  the  relative  simplicity  of  the  criterion 
aids  in  making  it  one  of  the  most  widely  us.id  failure  theories. 

Later  development  in  the  area  of  shear  failure  followed  closely 
the  work  in  the  theory  of  plasticity .  Usually ,  the  failure  theories 
are  expressed  in  terms  of  stress  invariants.  For  example,  Von  Mises' 
criterion  [5],  [6] 

j2  =  ^  [(cl  -  03)2  +  (03  -  cri)2  +  (ai  -  0's)2 ]  =  (2-15) 

where  J2  is  called  the  second  stress  invariant  and  k  is  a  material  con¬ 
stant  associated  with  the  yield  stress  in  simple  shear.  Other  criteria 
are  associated  wi*-h  Tresca  [5],  [6],  Freudenthal  [7]  and  Bresler  and 
lister  [8] 

(B)  Phenomenological  Failure  Theories 

This  category  includes  Griffith  theory,  McClintock  and  Walsh  s 
modified  Griffith  theory,  and  Murrel's  generalized  Griffith  theory  in 
three  dimensions . 
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Griffith  suggested  that  there  are  a  great  number  of  micro-cracks 
existing  in  a  brittle  solid.  Failure  occurs  when  the  existing  micro- 
cracks  spread.  The  condition  of  crack  extension  was  proposed  to  depend 
on^energy^balance ;  i.e.,  if  an  increase  in  crack  sice  (length  leads  to 
a  diminution  of  the  total  free  energy,  the  system  becomes  ^stable  and 
the  crack  spreads.  The  uniaxial  critical  stress  a0  (tensile), 

to  be 


for  plane  stress 


mt  > 

If  \  ,  for  plane  strain 

Vrtcci  -  vj;  ' 


(2-16) 


where  v,  Y,  are  the  Poissons'  ratio  and  Youngs'  modulus  of  the  material, 
7  is  the  surface  energy  of  the  material,  and  2c  is  the  leng  o 
crack.  The  maximum  stress  at  the  crack  tip  is  found  o  e 

(2-17) 


ext  =  Ob 


where  p  is  the  radius  of  curvature  at  the  crack  tip. 

The  Griffith  theory  has  been  extended  to  the  case  of  a  biaxial 
stress  field  with  an  elliptical  flaw .  It  takes  the  form 


(cij.  -  cx3 )2  -  8cto(ai  +  ^3) 


(2-18) 


when  ax  +  3cx3  >  0  and  cxto  is  the  tensile  strength  of  the  material.  For 
ax  +  3ct3  <  0,  the  failure  criterion  is 

03  -  ato  <  0  (2-19) 

The  Mohr  envelope  for  eq.  (2-l8)  can  be  shown  to  be 

T2  =  4cxto(cr  +  CTto) 

which  is  a  parabola. 

Under  compression,  the  Griffith  cracks  can  be  closed.  Gri^-ik 
theory  was  modified  by  McClintock  and  Walsh  to  account  for  the  fnc  ion 
along  the  crack  surfaces.  The  failure  criterion  takes  the  form 

(2-21) 


(2-20) 


-d'(°i  +  °3)  +  (°1  "  a3^  s/i  +  d'-  -  ^tc 


where  p’  is  the  coefficient  of  friction  along  the  crack  surfaces.  This 
can  also  be  written  as 


d!  = 


—  + 


-J 1  +  p'2  +  p1 


>/l  +  p’2  -  P*  Yl  +  p 


7? 


cr3 


(2-22) 


P 


Comparing  eq.  (2-22)  with  (2-6),  it  is  easily  seen  that  the  two  results 
will  be  the  same  if 

h'  =  M  >  2crto  =  T  (2-23) 

Therefore,  if  again  we  let  n  =  tan  cp,  the  modified  Griffith  failure 
criterion  can  be  e.cpressed  as 

t  =  2at0  +  her  (2-24) 

Graphically,  the  Griffith  theory,  eq.  (2-19)  and  the  modified  Griffith 
theory,  eq.  (2-24)  are  shown  in  Fig.  2-4,  where  n  is  assumed  to  be 

0.8  <  \i  <  l.o. 


Fig.  2-4  Griffith  Theory  (G.T.)  and  Modified  Griffith 
Theory  (M.G.T.) 


From  the  above,  it  is  seen  that  the  Griffith  theories  can  predict 
the  tensile  strength  of  a  brittle  material  in  terms  of  crack  length, 
elastic  constants  and  surface  energy.  It  is  capable  of  predicting  the 
stresses  for  failure  subject  to  a  two-dimensional  stress  state.  How¬ 
ever,  the  Griffith  theories  cannot  evaluate  oco  and  erto  directly;  thus 
the  theory  xs  not  applicable  to  the  failure  of  rock  material. 

Murrell  [9]  extended  the  Griffith  theory  to  three  dimensions.  The 
failure  surface  was  found  to  consist  of  a  portion  of  a  paraboloid  of 
revolution,  symmetrical  about  the  line  01  =  02  =  03  and  passing  through 
the  o  igin.  The  equation  is 


T—f-  ■  J'VP.4  -■pT  --  JWLMU»1»JUPI W1  11  1  J  IJ  IUUL.  i  II  H  I  JMWII Mil'll .  Wl  «BW»pWWfJp 


(a2  -  a3)?  +  (a3  -  ax'f  +  (oi  -  a2)?  =  24atota.  +  a2  +  a3)  ^  ^ 

(C)  Statistical  Failure  Theories 

The  statistical  models  propose  that  the  strengths  of  materials 
are  decided  by  the  weakest  link  of  the  random  distribution  of  Griffith 
flaws  in  the  materials.  Thus  there  will  be  a  distribution  of  strengths 
in  a  given  specimen  because  a  different  stress  is  required  to  fracture 
a  specimen  at  one  or  another  point.  Several  papers  based  on  this  con¬ 
cept  are  reviewed  in  Brady  [2],  and  will  not  be  further  discussed  here. 

2-3  LITERATURE  APPLICABLE  TO  LA/ER  CUTTING 

There  is  considerable  literature  available  that  applies  to  the 
problem  of  rock  drilling  mechanics  in  both  the  theoretical  and  experi¬ 
mental  areas  (see  for  example,  Maurer  [10 ],  Dutta  [11],  Mahban  [12], 
Hustrulid  [14],  Cheatham  and  Gnirk  [14],  Pariseau  and  Fairhurst  [15  ], 
etc.).  Most  of  these,  however,  relate  to  percussive  drilling  configu¬ 
rations  and  have  little  application  to  layer  cutting.  Literature 
directly  concerning  layer  cutting  of  rock  is  rather  sparse.  There  have 
been  numerous  studies  in  metal  cutting  problems,  however,  that  can 
provide  a  starting  point  for  some  analysis .  Literature  concerning  the 
metal  cutting  mechanics  will  be  reviewed  first .  Then  the  literature 
pertaining  to  the  layer  cutting  mechanics  of  brittle  materials  will  be 
reviewed. 

Merchant  [16]  considered  simple  models  of  metal  cutting  processes, 
which  included  the  "orthogonal  cutting,"  with  a  straight- edged  cutting 
tool  moving  relative  to  the  work  piece  in  a  direction  perpendicular  to 
its  cutting  edge,  and  the  "oblique  cutting,"  with  the  cutting  edge 
oblique  to  the  direction  of  relative  motion  of  tool  and  work  piece. 

Several  assumptions  were  made  for  the  orthogonal  cutting:  (i)  the 
chip  is  continuous  with  mass  density  p  unchanged  in  the  cutting  process; 

(ii)  shear  failure  occurs  on  a  shear  plane  As  at  an  angle  cp  (see 
Fig.  2-5  for  the  geometry,  and  Fig.  2-6  for  the  shear  failure  model); 

(iii)  the  chip  is  in  static  equilibrium,  as  shown  in  Fig.  2-7;  and 

(iv)  the  cutting  force  required  R,  is  equal  to  the  shear  strength  of 
the  material  t0,  multiplied  by  the  area  of  the  shear  plane  As. 

To  attack  this  problem,  first  the  chip  thickness  ratio  tc/t0  was 
geometrically  related  to  the  shear  angle  \|r  and  the  rake  angle  a.  The 
shear  strain  in  the  chip  can  also  be  related  to  vjr  and  a.  Finally  the 
chip  velocity  Vc  may  be  related  to  the  tool  velocity  VD.  Applying 
force  equilibrium  to  the  chip,  Merchant  was  able  to  obtain  the  various 
force  relations.  He  obtained: 

|  =  tan(cp  -  a)  (2-26) 
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MM 


0  =  tan-1  (i 

q  =  friction 

coefficient 


Fig«  2-7  Merchants'  force  components  on  a  chip 
where  Q  =  thrust  force  and  P  =  cutting  force.  Also, 

p 

To  =  [sin  cos  cp  -  tan(cp  -  a)  sin‘°v] 


(2-27) 


or 


(2-28) 


P  „  TqAo  -^-.c.°s(°  -  & - 

sin  ^  cos  (cp  +  ijr  -  a ) 

where  tq  is  the  shear  strength  of  the  material. 

From  these  results,  one  can  see  that  since  P,  Q,  Aq,  cp,  and  a  can 
he  measured  experimentally,  t0  can  be  determined  from  eq.  (2-27). 
Alternately,  if  t0  is  known,  and  if  oc,  <p,  'j/  and  A^  are  measured,  the 
cutting  force  required,  P,  can  be  calculated  from  eq.  (2-28).  Figure 
2-8  shows  I/2tqAq  versus  cc  using  <p  -  tan-1n  as  a  parameter. 

From  Fig.  2-8,  one  can  see  that  the  cutting  force  required  to 
break,  a  layer  of  metal  is  directly  proportional  to  the  area  of  cutting 
and  the  shear  strength  of  the  material.  It  also  shows  that  the  cutting 
force  required  is  larger  for  p  smaller  rake  an„lc  a.  Furthermore, 
friction  between  the  tool  and  the  chip  also  increases  the  cutting  force 
needed. 
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Fig.  2-8  Cutting  force  P  versus  rake  angle  Qt 
relations  from  Merchants'  theory 


There  axe  certain  deficiencies  in  the  theory.  For  one,  it  is  .unpxied 
that  the  cutting  force  P  is  independent  of  the  cutting  speed  Vc,  which  at 
least  implies  that  the  dynamic  effect  in  metal  cutting  is  not  considered. 
Moreover,  the  assumption  of  a  single  shear  plane  does  not  agree  with  ex¬ 
perimental  observations.  The  model  proposes  an  external  surface  which  is 
not  smooth  at  the  junction  of  the  chip  and  the  original  surface.  Theo¬ 
retically,  an  infinite  shear  gradient  exists  in  this  shear  plane  and  an 
infinite  acceleration  of  the  chip  is  involved,  which  is  physically  im¬ 
possible.  A  transition  zone  of  plastic  deformation  (chip  formation  zone) 
is  definitely  necessary  to  modify  the  original  theory. 

Zorev  [17]  described  several  modified  models  which  included  the 
transition  zone.  He  discussed,  for  example,  the  models  of  Briks  [l8] 
and  Krivoukhov  [19].  He  analyzed  a  modified,  model  with  chip  deforma¬ 
tion  zone  which  provides  for  smooth  joining  of  outer  boundaries  of  the 
cut  layer  as  shown  in  Fig.  2-9. 

The  chip  formation  zone  takes  the  shape  of  a  fan.  The  stressed 
state  arising  is  characterized  by  the  fact  that  the  lines  of  the  con¬ 
stant  tangential  stresses  are  a  fan-shaped  pattern  of  straight  lines 
passing  through  the  cutting  edge.  The  initial  and  final  boundaries  of 
the  plastic  zone,  OL  and  OM,  are  determined  from  the  plasticity  condi¬ 
tions.  It  is  assumed  that  the  transformation  from  the  layer  being 
removed  into  a  chip  occurs  by  successive  shear  along  a  family  of  slip 


51 


Chip 


Fig.  2-9  Chip  formation  zone 


lines  passing  through  the  cutting  edge  0.  Shear  along  the  second  family 
of  slip  lines  is  considered  to  be  comparatively  small,  and  therefore 
is  neglected.  The  direction  of  the  shear  lines  is  taken  to  coincide 
with  the  direction  of  the  lines  of  maximum  tangential  stress. 

As  in  Merchant's  treatment,  again  the  shear  train  is  related  to 
geometry.  The  shear  stress  acting  on  a  shear  plane  is  related  to  the 
resultant  force  on  the  chip  by  considering  the  equilibrium  of  the  chip, 
neglecting  the  moment  acting  on  the  chip.  Experiments  have  shown  that 
the  shear  stress  obtained  for  the  initial  chip  formation  boundary  is 
2QPjo  too  large.  Zorev  modified  the  shear  angle  y  to  \|r*  to  offset  the 
discrepancy.  The  necessity  for  this  modification  is  that  the  actual 
slip  lines  are  not  necessarily  straight  as  proposed  earlier.  Actually, 
they  are  more  like  those  shown  in  Fig.  2-10. 

By  these  procedures,  Zorev  was  able  to  obtain  a  more  accurate 
cutting  force,  the  final,  the  outer,  and  the  initial  boundaries  of  the 
chip  formation  zone,  the  principal  axes  of  deformation  and  the  true 
deformation  of  the  chip.  The  resultant  force  on  the  tool  has  the 
expression: 

R  =  - t - , - ?QtJQta -  (2-29) 

cos(|i*+  w)  sin  ^^(0.9  -  0.5  cot  Vi) 

where  tyi*  is  the  modified  angle,  which  is  about  0.8  ^ . 
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Fig.  2-10  Lines  of  slip  in  the  chip  formation  zone 


Evans  [20]  considered  the  problem  of  plowing  brittle  materials  (coal) 
with  a  symmetric  shaped  wedge.  A  scheme  similar  to  Merchants'  metal 
cutting  mechanics  was  proposed  in  that  the  chip  is  in  static  equilibrium. 
However,  chips  formed  in  plowing  brittle  material  are  discontinuous  and 
the  fracture  mode  is  assumed  essentially  tensile  in  nature. 

The  model  to  be  treated  is  a  right-angled  quarter-space  of  brittle 
material,  as  shown  in  Fig.  2-11.  The  breakage  curve  is  assumed  to  be  a 
circular  arc  initiated  from  the  tip  of  the  wedge.  The  chip  is  acted  upon 
by  three  resultant  forces.  The  tensile  force  T  is  the  resultant  of  the 
tensile  stress  acting  along  the  circular  arc.  The  force  R'  making  an 
angle  cp  with  the  normal  of  the  tool  surface  is  exerted  by  the  tool. 
Thirdly,  a  concentrated  force  S,  acting  on  point  D,  produces  static 
equilibrium  of  the  chip. 


Fig.  2-11  Evans'  plowing  pioblem 


Immediately  before  the  fracture  occurs,  the  stress  acting  on  the 
circular  arc  is  assumed  to  reach  the  tensile  strength  of  the  material, 
(rto^-*  Thus,  T  can  be  expressed  as  a  function  of  deoth  of  cut  d,  and 
angle  0.  Equilibrium  of  the  chip  gives  P,  the  horizontal  thrust  as  t\ 
function  of  9  and  a,  multiplied  by  thus 

P  =  2R*  sin  a 


_ 2crtnd.(l  +  cot^a)  sin2(a  +  cp)  _  (2-30) 

h/d  +  2  cot  9  sin(a  +  9)  cos(a  +  cp)  -  2  sin- (a  +  cp) 


The  assumption  of  minimum  work  done  by  the  horizontal  thrust  P, 
i.e.,  oP/dO  =  0  is  employed.  In  addition,  the  findings  of  Evans  and 
ftirrell  [21]  are  used;  namely,  that  the  force  required  to  push  a  wedge 
in  normal  to  a  face  can  be  accounted  for  in  terms  of  a  "penetration 
resistance"  closely  allied  to  the  compressive  strength  <rc0  of  the  coal. 
The  resulting  two  supplemental  equations  are  then 

||  =  cot20  sin(a+  qp)  cos(a+qp)  +  cot  9  [|  2  sin2(a+q>)] 

-  sin(a+  qp)  cos(a  +  <p)  =  0  (2-31) 

P  =  2acoh  (2-32) 

From  the  last  three  equations,  one  can  solve  for  P  as  a  function  of 

a,  with  qp  and  K  =  (K  10  for  brittle  materials)  as  parameters. 

The  results  are  plotted  in  Fig.  2-12.  The  angle  0  can  also  be  calcu¬ 
lated.  From  Fig.  2-12,  it  is  readily  seen  that  the  effect  of  parameter 
K  is  small  (for  qp  =  0° ) .  This  ,  in  turn,  implies  that  the  penetration 
of  the  wedge  before  breakage  takes  place  is  small.  If  one  assumes  that 
the  penetration  is  negligibly  small,  the  expression  for  P  is  simplified 
to  be 


P  -  2a-fc0d 


sin(q  +  g?) 

1  -  sin(q  +  cp) 


(2-33) 


Equation  (2-33) 
20  and  30° . 


was  used  to  plot  the  results  in  Fig.  2-12  for  <p  =  10, 


a  (degree) 


Fig.  2-12  The  horizontal  thrust  force  P  versus 
the  half-wedge  angle  a 


Experiments  were  carried  out  by  Evans  on  several  coals.  The 
results  compared  favorably  to  his  theoretical  results.  However,  a 
discrepancy  was  noted  in  the  theoretical  and  experimental  tensile 
strengths.  Evans  attributed  this  discrepancy  to  the  assumption  that 
the  mode  of  failure  is  a  sudden  uniform  disruption,  which  does  not  take 
into  account  the  crack  propagation  effect.  He  proposed  a  way  to  com¬ 
pensate  for  this  effect  by  assuming  that  the  tensile  stress  varies  along 
the  circular  arc  according  to  the  law  ato(x/£)n.  Comparing  with  the 
experimental  results,  he  found  that  n  =  8  gave  a  good  result. 

_  ^sns  t22]  continued  his  study  on  plowing  to  consider  the  effect 
°  1Jn^eSS  of  the  wed€e-  The  geometry  of  the  problem  as  shown  in 
tig.  2- 13.  He  assumed  that  the  compressive  force  on  the  half-band  of 

the  wedge  needed  to  initiate  the  crack  is  of  the  f6rm  gp  (^fd,  where 

ac  is  the  compressive  stress  and  m  is  a  generalized  parameter  to  be 

^!™ined  ("hich  was  found  to  be  in  the  vicinity  of  ».  Again,  the 
cheme  of  a  tensile  failure  along  a  circular  arc  was  assumed.  The 
ip  equilibrium  condition  plus  the  minimum  work  hypothesis  enabled 

c^lci°a^dntai  f°^6  °n  the  W6dge  needed  to  “itiate  fracture  to  be 
initiation;  he^aLed'6  °f  penetrQtion  cf  the  before  crack 


=  2at0d  £ 


sin(a+cp) 


2  sin  9  cos((9  +  a+qj) 


1  +  sxn  fl  sin(o;+  cp)J 


— •->  -A-  ^  t  ^  I  Vj. 

cos(t)  +  U  +  cp) 


(2-34) 


Fig.  2-13  Evans'  blunt  wedge  plowing  model 


where  the  angle  0  is  given  by 


cos  (20  +  a  +  <p)  om-i/crc.Vb\m  , 

3  -  cos  20  ~  ^  fejVdj  cos(a  +  (2- 

Experimental  results  on  Cwintillery  coal  were  given  for  comparison. 
He  found  that  if  m  =  2/3,  his  theoretical  results  were  very  close  to 
the  experimental  ones,  as  shown  in  Fig.  2-l4. 


d 

Fig.  2-lU  Comparison  between  experimental  and 

theoretical  results,  continuous  cutting 
with  blunted  picks 


Recently,  Nishimatsu  [23]  has  treated  a  two-dimensional  rock  cutting 
model.  He  observed  that  the  failure  process  of  rock  cutting  is  as  given 
in  Fig.  2-15.  The  crushed  zones  can  be  divided  into  two  categories: 

(i)  the  major  crushed  zone  and,  (ii)  the  secondary  zone.  The  tool  first 
chips  off  several  small  pieces  (minor  chips),  then  a  large  fragment  of 
rock  is  fractured  (major  chip).  This  process  repeats  for  the  entire 
cutting  process. 


With  the  above  observation,  the  author  proceeded  to  analyze  the 
major  chip  crushing  model  as  shown  in  Fig.  2-16. 


Fig.  2-16  Stress  distribution  and  cutting  forces  for 
orthogonal  rock  cutting 


La‘r^  ^rtoVetrlbnti0n  ^  ■*"«« 


P  =  P0 


(db  '  *)' 


(2-36) 


where  (A, 0)  is  the  polar  coordinate  of  a  point  p  on  line  OB  ardn  is 
2^edethatdi^ribf i0n  fSCt0r  t0  te  FuX^oChe 

assumed  that  the  stress  p  has  a  constant  direction  along  AB  and  the 
reactant  force  E  will  be  In  equilibrium  with  the  acting  tecfp  toon 
the  tool,  in  no  doing,  p0  in  eq.  (2-36)  can  be  found  to  be  ’ 

sin  0\n+1  „ 


Po  .  -(n  +  1)(^) 


(2-37) 

S^ce^r^r^6  ?■,  the  tanSential  normal  components  along  AB. 
?a^Ure  wm  be  ^tiated  at  the  maximum  stress  concentration 

*t-res  ’  he  ^ avaer -Ooulomb -Mohr  criterion  can  be  applied  to  the  above 
stress  components  at  point  0  to  obtain 


F  =  - 


n  +  1  sin  e  tan  k  sin(0  -  a  +  cpj  -  cos(e  -  a  +  <p) 


(2-38) 


where  k  is  the  internal  friction  angle  of  the  material  and  -v0  is  the 
ahear  strength  of  the  material.  The  direction  of  the  line  0?  on  Sich 
e  allure  stress  acts  should  be  the  line  which  makes  the  cutting  force 
a  minimum.  Thus  the  resulting  cutting  force  is  found  to  be 


Tnt 


cos  k 


n  +  1  1  -  sin(k  -  a  +  qp) 


(2-39) 


Fc  1S  a  ^inear  Action  of  the  depth  of  cut  t.  The  horizontal 
cutting  force  P  and  the  vertical  thrust  Q  are 

P  =  Fc  cos  (qp  -  a) 

Q  =  Fc  sin(cp  -  a)  (2-t*o) 

Nishamatsu  also  gave  experimental  values  of  P  and  Q  versus  rake  angle 
a  for  sandy  tuff  and  cement  mortar.  The  forces  were  measured  with  a 
circular  plate-type  dynamometer.  He  found  that  qp  is  a  linear  function 
rake  angle  a  as  shown  in  Fig.  2-17.  The  experimental  cutting  force 

t[flc^ion  of  b  as  Predicted  by  the  theory.  The  stress 
distribution  factor  n  is  found  to  have  a  value 


n  =  11.3  -  0.18  a 

where  a  is  expressed  in  degrees. 


(2-4l) 
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^8*  2-17  The  angle  of  friction  of  rock  cutting 


vensson  and  Freshwater  [2U]  gave  some  experimental  studies  on  rock 
cutting  by  employing  a  photoelastic  technique  on  Catalin  800  and  poly- 
styrene.  Fringe  patterns  for  various  cutting  processes  rere  obtained 
and  he  was  able  to  deduce  the  stress  distribution  during  cutting.  From 
his  studies,  he  made  the  following  conclusions: 

v  -+4.?  Phol,oelastic  stress  analysis  of  chip  formation  in 

^^erialS’  Catalin  800  311(1  Polystyrene,  has  established 
that  chip  formation  is  a  process  of  crack  initiation  and  propa- 
gation  due  to  tensile  stresses  normal  to  the  crack  path.  No 
justification  could  therefore  be  found  for  the  application  of 
any  analysis  based  on  a  shear  theory  of  failure. 

The  .results  show  that  fracture  begins  at  a  very  small  tool 
penetration  and  the  crack  propagates  slowly  whilst  the  cutting 
force  is  increasing.  At  sane  critical  crack  length,  propaga¬ 
tion  becomes  more  rapid  and  the  cutting  force  decreases 
rapidly.  The  photoelastic  analysis  shows  there  is  an  intense 
concentration  of  tensile  stress  in  the  immediate  vicinity  of 
the  incipient  crack.  This  stress  rapidly  decays  and  over 
he  latter  portion  of  the  completed  crack  path,  the  normal 
stress  becomes  compressive.  The  stress  distribution  along 
'rf.  cfack  Path>  except  for  regions  near  the  tool  tip,  agrees 
with  the  theoretical  results  obtained  by  Shepherd. 
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The  crack  path  is  not  a  principal  stress  trajectory 
and  shear  stresses  of  the  same  order  of  magnitude  as  the 
normal  stresses  act  along  most  of  its  length.  Near  the 
tool  tip,  however,  the  stresses  normal  to  the  crack  path 
are  principal  stx esses.  Equilibrium  of  the  incipient 
chip  is  therefore  maintained  by  the  tool  force  acting  on 
the  chip  and  the  normal  and  shear  stress  distribution 
along  the  crack  pith.  This  stress  distribution  is  far 
more  complex  than  that  assumed  by  Evans.  Other  observa¬ 
tions  made  during  the  investigation  showed  that  the  crack 
path  was  not  a  circular  one.  As  the  crack  length  in¬ 
creases,  the  bending  effect  of  the  tool  force  becomes 
more  pronounced.  This  causes  a  rotation  of  the  principal 
stress  directions  such  that  the  radius  of  curvature  of 
the  resulting  crack  diminishes  as  it  penetrates  further. 

There  are  two  aspects  which  are  of  importance  in 
studying  the  machining  of  brittle  materials. 

(a)  The  maximum  force  required  to  initiate  the  crack, 

(b)  The  energy  required  to  remove  a  unit  volume  of 
material. 

From  the  results  of  these  investigations  it  appears 
that  the  maximum  force  is  largely  dependent  of  the  depth 
of  cut  but  would  possible  depend  upon  depth  of  cut  in 
addition  to  tool  angle. 

Further  work  is  necessary  to  establish  the  mechanism 
of  crack  initiation,  the  concept  of  a  critical  crack 
length  (which  governs  the  maximum  force)  and  the  effect 
of  tool  angle  on  these. 

Other  methods  of  rock  cutting  which  may  be  applied  to  layer  cutting 
analysis  include  crater  indexing  with  a  sharp  wedge  or  with  a  blunted 
tool.  Garner  [25]  gave  possible  fractures  for  crater  indexing  (Fig. 
2-l8)  and  Cheatham  [26]  gave  possible  slip  lines  for  indexing  in  a 
plastic  rock  (Fig.  2-l§). 
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Possible  tensile  fracture 


Possible  shear  fracture 


Fig.  2-19  Cheathams'  possible  slip  lines  for 
indexing  in  a  plastic  rock 
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3.  SONIC  ROCK  DRILLING 


Study  of  the  basic  drilling  configuration  was  begun  in  the  first 
year  of  sonic  rock  cutting  and  has  been  previously  reported.*  In  the 
current  phase  of  work,  the  basic  sonic  drilling  techniques  previously 
developed  were  utilized,  with  certain  modifications  and  refinements. 

The  basic  technique  was  to  force  a  P-11  transducer  vertically  downward 
with  a  static  force.  A  small,  bouncing  mass  drilling  tool  at  the  tip 
of  the  transducer  performed  the  cutting  action  through  repeated  impacts 
with  the  rock  surface.  This  impact  coupling  drilling  technique  also 
has  been  previously  described. 

The  basic  parameters  of  the  drilling  process  have  been  identified 
as  the  transducer  drive  voltage,  the  mass  and  geometry  of  the  impact 
tool,  the  static  force  and  the  rock  properties.  A  limited  number  of 
tests  were  conducted  to  establish  optimum  static  force  values.  Partic¬ 
ular  attention  was  given  in  this  phase  of  work  to  various  tool  designs. 

The  parameter  of  drive  voltage  was  held  constant  during  most  of  the 
drilling  tests.  Some  attention  was  given  to  using  variable  frequency 
versus  fixed  frequency  drive  to  the  transducer.  This  phase  of  work 
continued  drilling  tests  in  Indiana  limestone  and  initiated  drilling 
tests  on  granite. 

3-1  APPARATUS 

The  basic  drilling  apparatus  has  been  described  and  shown  in  an 
earlier  report.**  Some  modifications  were  incorporated  in  this  apparatus. 
First,  the  external  static  force  applied  at  the  nodal  point  of  the 
transducer  by  dead  weights  was  replaced  by  a  pneumatic  cylinder  capable 
of  applying  static  forces  of  up  to  200  pounds.  Secondly,  an  air-vacuum 
system  has  been  incorporated  for  removal  of  rock  debris.  This  system 
includes  an  air  hose  on  the  one  side  of  a  cylindrical -shaped  plastic 
container  at  uhe  transducer  tip,  and  a  vacuum  system  on  the  other  side. 
The  advantages  of  this  set-up  are  (l)  increased  drilling  penetration 
rate  due  to  debris  removal  from  the  cutting  surface,  (2)  considerable 
noise  reduction  in  the  drilling  process,  (3)  prevention  of  dust  from 
escaping  to  the  laboratory,  and  (4)  removing  some  of  the  heat  from  the 
sonic  tool.  Although  items  (2)  through  (t)  are  desirable  benefits,  the 
main  objective  of  the  air -vacuum  system  was  for  rock  debris  removal.  A 
schematic  of  the  modified  system  is  shown  in  Fig.  3-1* 


*  See  Annual  Technical  Report  on  previous  Contract  H0210010,  "Fundamental 
Studies  in  the  Use  of  Sonic  Power  for  Rock  Cutting,"  December,  1971. 

-**Ibid.,  p.  37- 
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Fig.  3-1  Drilling  apparatus 
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nf  *  was  devised  for  automatically  measuring  and  recording  denth 

of  penetration  versus  time  during  drilling  Tho  aePtn 

motion  was  indicated  by  the  !  ^stance  of  transducer 

This  signal  was  in  l  axactrical  output  of  a  ten-turn  potentiometer, 

there  to\n  xy  Lnrn  *  ^  t0  ^  electronic  Processing  unit  and  from 
siZtl  foTthi  ™ ;  The  processing  unit  also  generated  a  time  base 
3^  recorder .  A  simple  schematic  of  the  system  is  shown 


3-2  TOOL  DESIGN 

A  number  of  tool  configurations  were  developed  during  the  course 
of  work.  It  will  be  noted  that  tools  used  in  the  earlier  phase  of  work 
on  rock  cutting*  had  large  amounts  of  flat  surface  area  on  the  impact 
face  of  the  tools.  Such  designs  were  highly  inefficient,  since  rock 
penetration  occurred  mainly  by  crushing  action.  Tool  designs  during  the 
current  phase  of  work  were  aimed  at  reducing  flat  areas  and  increasing 
edges  where  stress  concentrations  could  occur. 

dhown^Uf  f,SeVSa  °f  the  l0?15  tested  for  Intone  drilling  are 
Umpc+o  j  3‘  ^  was  Generally  found  that  finely  pulverized 

limestone  would  compact  in  the  recesses  and  material  egress  holes  in 
^hese  configurations  Photographs  of  a  second  series  ft  tools  is  shown 
best  Jt  W?S  found  that  the  configuration  of  Fig.  3-4(b)  yielded 

LTmhL?/?  tS  111  limestone-  Finally ,  a  configuration  des^ned 

S^e  Sd  not  ne^t62'^1^  "t™**  "  ShOWn  in  Fig*  ^  UnfortunSety, 
p  rnu  cutting  tests  to  be  conducted  using  this  tool. 

3-3  DRILLING  RESULTS  IN  LIMESTONE 

..  A  series  of  drilling  tests  in  limestone  were  carried  out  to  assess 

Z  Zli%VT^°z  1vtatic  force  on  «“*«■  ZTls 

Fig.  \ ^  Fi6*.3-3.  The  results  are  shown  in 

..  ®  v,1  ^ere  ds  a  marhed  lack  of  uniformity  in  these  results  in  that 
4yl  Charac^eristics  of  each  curve  differ.  Thus,  Tool  #1  exhibits  an 

^Vi°rTVi^inCreaSing  f°rCe’  ha™6  a  P~~d Z  Tt 

^  ShoWS  a  moderate  but  continually  increasing 
Strutf-?ePt?  ^ u  increasinS  f°rce,  while  Tool  #3  shows  marked  irc^fases 

ahout^lUO^pounds^161,  f°rCeS ’  ^  ^  eXhlbits  *  Wea*  P°^  ab 

T,5eDnain  c°ncdusion  ^awn  from  these  results  was  that  static  force 
aPPear  to  be  a  highly  sensitive  controlling  parameter  of  the 
drilling  process.  Consequently,  in  most  later  drilling  tests  static 

^tT^SlSlfr  200  ?OUndS5  Was  the  largest  valua  obtainable 

the  available  air  supply,  pneumatic  cylinder  system. 


*Ibid. ,  p.  23,  29,  33,  34. 


67 


*  jj  u  ui^i  u;  h  i  MU41U 


(c)  Tool  #7 


(d)  Tool  #8 


(e)  Tool  #9 
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on  drilling  for  various  tools 


A  specific  set  of  tests  was  run  to  assess  the  effectiveness  of  the 
air -vacuum  system  on  improving  drilling.  A  specific  tool,  static  force 
value  and  drive  voltage  was  selected,  and  depth  of  penetration  versus 
time  measured  with  and  without  the  air-vacuum  system  operating.  The 
results  are  shown  in  Fig.  3-7*  Without  the  air-vacuum  system  to  remove 
drilling  debris,  it  is  seen  that  a  point  is  reached  at  which  penetration 
practically  ceases  (lower  curve  of  Fig.  3-7  at  about  90  seconds).  With 
debris  removal,  it  is  seen  that  a  steady  penetration  rate  is  achieved. 

In  both  cases,  initial  penetration  rates  are  quite  rapid  and  about 
comparable.  This  is  to  be  expected  in  the  early  stages  wlen  both  the 
lack  of  debris  and  the  aspect  of  cutting  into  the  initially  flat  sur¬ 
face  combine  to  speed  cutting. 

The  influence  of  variable  frequency  versus  fixed  frequency  drive 
to  the  transducer  was  also  assessed,  with  the  results  shown  in  Fig. 

3-8.  The  basis  for  conducting  this  test  was  the  known  difficulty  in 
achieving  a  transducer  resonance  frequency  at  the  fixed  power  supply 
frequency  at  high  voltage  operating  conditions.  Secondly,  there  was 
also  the  known  problem  of  shifting  transducer  frequency  during  operation 
due  to  temperature  effects.  A  drilling  test  0.0  fixed  frequency  was 
conducted.  Then,  a  drilling  test  was  conducted  during  which  the  drive 
frequency  was  varied  during  the  test  in  a  manner  to  maintain  peak  power 
into  the  transducer. 

As  is  evident  from  Fig.  3-8,  the  difference  in  penetration  becomes 
fairly  significant  for  increasing  time  (say,  beyond  30  seconds).  The 
probable  cause  of  fall-off  of  the  fixed  frequency  data  is  transducer  de¬ 
tuning  resulting  from  temperature  increase  within  the  transducer. 

Figures  3-9  and  3-10  show  penetration  tests  for  several  of  the 
cutting  tools  (Tools  #3 5  7,  8  and  6).  The  data  shown  was  obtained  with 
the  automatic  recording  equipment  previously  described  and  indicated 
in  Fig.  3-2.  All  of  the  results  show  the  similar  characteristic  of 
initial  rapid  penetration  (up  to  120  in. /min  initial,  rate)  and  then  a 
reduction  to  a  steady  cutting  speed.  The  best  results  obtained  were 
those  for  Tool  #6  (Fig.  3-10,  lower  trace),  which  had  the  following 
performance  characteristics: 


Material:  Indiana  Limestone 
Hole  diameter:  1  1/4  in. 

Drilling  rate  (steadv):  11  1/2  mi/min 
Volume  removal  rate:  14  in.3/min 
Specific  energy:  10,700  W-sec/in.3 

3-4  DRILLING  TESTS  IN  GRANITE 

As  in  layer  cutting,  difficulties  were  encountered  in  cutting 
granite.  Unlike  the  layer  cutting  work  discussed  previously,  where 
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Fig.  3-8  Depth  of  penetration  versus  time  for  fixed  gnri 
variable  frequency  operation 
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s-sae  success  was  finely  achieved,  the  efforts  at  drilling  in  granite 
were  essentially  negative  throughout.  The  various  cutting  tools  developed 
or  limestone  drilling  had  little  effect  on  granite  specimens.  Hardened 
steel  tools  dulled  in  matters  of  seconds  or  fractured  at  stress  concen¬ 
tration  points.  Drilling  tests,  while  frequently  ineffective,  we>*e 
nearly  always  spectacular.  The  region  of  the  tool-rock  interface  would 

r-^acn  a  red  heat  and  showers  of  sparks  would  be  thrown  from  the  impact 
area. 


Several  te.Jts  were  done  on  granite  drilling  using  slightly  unusual 
types  of  tools.  The  main  tools  were  carbide-tipped  steel  and  pieces  of 
ceramic  The  ceramic  tools  v/ere  made  of  zirconia.  It  was  possible  to' 
drill  about  1/2"  into  granite  boring  a  3/4"-diameter  hole.  The  ceramic 
would  chip  into  smaller  pieces.  It  appeared  the  drilling  was  achieved 
by  eating  and  fracturing  of  the  rock.  As  the  ceramic  was  impacted  by 
the  P-11,  the  interface  between  the  ceramic  and  granite  actually  became 
red  hot.  Next,  an  imoact  tool  with  two  carbide  inserts  was  tried,  with 
the  inserts  mounted  in  a  sloped  manner.  During  drilling  the  carbide 
tips  penetrated  the  granite  to  a  depth  of  3/16".  The  flat  surfaces  of 
the  tool  then  impacted  the  rock  ana  no  further  drilling  was  achieved. 

Drilling  tests  were  also  done  using  a  sharp  wedge  tool  made  of 
00  steel.  Before  using  the  transducer,  tests  were  done  on  a  static 
testing  machine.  With  a  few  thousand  pounds  of  force  it  was  possible 
to  force  the  wedge  tool  into  granite  about  3/16"  with  no  signs  of  tool 
deformation .  The  same  tool  was  then  used  with  a  P-11  transducer  drill¬ 
ing  into  granite.  The  tool  penetrated  about  3/l6"  with  considerable 
wear.  A  possible  method  of  reducing  this  wear  would  be  to  use  carbide 
tips  and  rotate  the  tool. 

The  end  of  the  project  work  period  arrived  before  tool  rotation 
schemes  or  additional  drill  configurations  could  be  tested.  For  example, 
*to°l  configuration  of  Fig.  3-5  was  not  completed  in  time  for  testing. 
Although  it  is  likely  that  these  modifications  would  have  somevrhat 
improved  granite  drilling,  it  is  questionable  whether  drilling  rates 
of  practical  significance  would  have  resulted. 
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4.  ANALYSIS  OF  ROCK  DRILLING 


The  feasibility  of  utilizing  sonic  energy  for  drilling  brittle 
materials,  such  as  concrete  or  rock,  has  been  under  study  for  some  time. 
The  basic  process  employs  a  sonic  transducer  mounted  on  a  support 
structure  and  a  slug  of  metal,  called  the  "tool",  as  shown  in  Fig.  4-1 
and  as  previously  described  in  this  report.  The  tools  are  of  a  general 
cylindrical  shape  with  various  tip  and  end  geometries.  Sonic  energy  is 
transmitted  by  means  of  impact  coupling,  or  in  other  words,  by  repeated 
impact  of  the  tool  between  the  transducer  and  load  at  rates  of  hundreds 
to  thousands  of  cycles  per  second.  The  tool  impacts  against  the  rock 
after  impacting  with  and  acquiring  energy  from  the  vibrating  transmission 
line .  Seme  energy  of  the  tool  is  transmitted  into  the  rock  and  results 
in  some  fracture  of  the  rock  surface.  The  tool  then  rebounds  from  the 
rock  surface  due  to  elastic  recovery  and  again  impacts  against  the 
transmission  line  tip.  Because  of  the  energy  withdrawal  from  the  trans¬ 
ducer  by  the  tool  during  tool-line  impact,  the  internal  energy  of  the 
transducer  is  at  a  lower  level  right  after  impact.  It  recovers  before 
the  next  impact  occurs . 

Although  the  impact  coupling  process  is  rapid  and  somewhat  random, 
it  is  a  continual  sequence  of  two  clearly  defined  events:  (l)  impact 
and  rebound  of  the  tool  against  the  vibrating  line,  and  (2)  impact  and 
rebound  of  the  tool  from  the  load.  Therefore,  a  rational  analysis  of 
impact  coupling  becomes  possible  and  can  be  divided  into  three  areas  of 
study  —  tool -line  impact,  tool-load  impact  and  transducer  analysis. 

Each  area  of  these  analyses  involves  several  controlling  parameters. 

The  mass  and  tip  geometry  of  the  tool,  the  rock  properties,  and  the 
energy  level  of  each  impact  affect  the  incident-rebound  velocity  rela¬ 
tionship  as  well  as  the  characteristics  of  the  impact  zone.  The  end 
geometry  and  mass  of  the  tool  and  the  geometry  and  vibration  amplitude 
of  transmission  line  tip  contribute  to  the  characteristics  of  tool¬ 
line  impact.  The  transducer  is  a  rather  complicated  energy  conversion 
system.  Many  parameters,  aside  from  the  supplied  terminal  voltage, 
affect  its  vibration  characteristics.  All  of  the  previously  mentioned 
parameters  interact  with  each  other;  any  one  parameter  can  easily  affect 
several  others  in  the  overall,  performance  of  the  sonic  system.  r.he 
static  force  applied  by  the  support  structure  is  another  factor  in 
determining  the  energy  transfer  of  the  impact  coupling  process. 

A  considerable  amount  of  work  has  been  done  previously  on  vaiv.ous 
aspects  of  the  sonic  impact  coupling  process,  with  a  fairly  complete 
summary  being  given  in  [l],*  pp.  184-186.  Similarly,  study  in  the 
area  of  transducer  characteristics  has  been  extensive ,  although  much 
is  yet  unknown  on  the  performance  of  these  devices.  McMaster,  Dettloff 
and  Minchenko  [2]  reviewed  the  basic  principles  of  piezoelectric 
materials,  while  Dettloff  and  Minchenko  [3]  studied  the  resonant  horn. 


■^Numbers  in  brackets  refer  to  references  at  the  end  of  this  chapter. 


The  overall  transducer  cnaracteristics  was  first  touched  on  from  the 
viewpoint  of  classical  circuit  and  transmission  line  theory  by  Hoffman 
and  Swartz  [4].  Following  this  line  of  approach,  Mahtan  and  Graff  [5] 
developed  the  transfer  function  of  a  sonic  power  system.  Graff  [6],  [7l> 
in  1969,  reviewed  the  problems  of  vibrations  of  crystals  and  stepped 
horn  resonators.  The  energy  dissipation  in  a  vibrating  transmission 
line  was  analyzed  by  Fretwell  [8]  and  Graff  [93*  Ma  [10],  as  part  of 
an  extensive  study  of  the  transducer,  has  measured  the  transient  response 
of  the  transducer  to  impact.  In  Section  5  t>f  this  report,  results 
obtained  by  Feng  on  transducer  analysis  are  given.  A  more  complete 
description  of  transducer  studies  appears  in  [11]. 

In  the  analysis  of  the  rock  drilling  process,  work  on  the  .tool- 
rock  and  the  tool -transmission  lines  has  been  previously  reported.* 
Additionally,  a  qualitative  explanation  of  the  general  interaction 
process  has  been  put  forth.**  In  the  first  section  of  this  chapter, 
certain  additional  information  on  tool-rock  impact  is  presented.  In 
the  second  section,  additional  results  are  given  in  the  tool-transmission 
line  aspects  of  the  process.  The  third  section  of  the  chapter  will  bear 
additionally  on  the  transducer  aspects  of  the  process. 

In  section  four  there  is  presented  an  analysis  of  the  effects  of 
various  parameters,  such  as  drive  voltage,  transducer  energy  storage, 
static  force,  tool-load  and  tool-line  parameters  on  drilling.  Finally, 

1  the  analyses  are  united  in  the  prediction  of  sonic  drilling  rate 
of  rocks.  A  numerical  example  is  given  at  the  end  of  the  chapter.  A 
comprehensive  treatment  of  all  aspects  of  the  problem  may  be  found  in 
[12]. 


4-1  TOOL  IMPACT  ON  ROCK- -SPECIFIC  ENERGY 

Previous  results  on  sonic  tool  impact  on  rock  have  appeared  in 
[l],  Chapters  1  and  4.  During  the  current  phase  of  work,  additional, 
tool-rock  impact  studies  were  conducted  using  apparatus  and  techniques 
described  in  the  cited  reference.  The  objective  of  the  work  was  to 
establish  specific  energy  data  for  rock. 

To  study  the  specific  energy  of  rock  removal  for  granite,  a  tool 
design  having  a  sharp  tip  made  of  carbide  inserts  was  used.  The  tool 
was  propelled  into  a  series  of  impacts  with  the  rock  surface  at  a 
specific  location.  The  total  energy  used  was  the  sum  of  all  the  indi¬ 
vidual  impact  energies,  and  the  volume  of  rock  removal  was  obtained  by 
dividing  the  total  weight  of  the  chips  and  dust  collected  by  the  average 
density  of  the  rock. 


*  Reference  [l].  Part  II,  Chapters  1,  4  and  5- 

**Ibid,  Chapter  6. 
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Figure  4-2  shows  the  relationship  between  incident  velocity  V0, 
specific  energy  eg,  and  volume  v  for  the  tool  impacting  on  granite. 

The  tool  design  is  also  shown.  A  linear  relationship  between  v  and  V0 
is  observed.  This  result  strengthens  the  belief  that  the  higher  the 
energy  level  per  impact  the  more  fracture  occurs  for  an  equal  amount  of 
work.  Rock  as  well  as  concrete  is  a  brittle  material .  It  is  weak  in 
both  tension  and  shear.  Fracture  occurs  when  the  tensile  and  shear 
stress  reach  a  certain  level.  High  energy  level  impact  would  cause  not 
only  more  fractures  but  also  cause  them  to  propagate  a  greater  distance. 
Therefore,  it  would  cause  larger  fractures  than  would  a  low  energy 
impact  and  save  energy  by  not  breaking  rock  into  smaller  fragments  than 
necessary. 


4-2  TOOL  IMPACT  ON  A  VIBRATING  TRANSMISSION  LINE 

Many  aspects  of  tool  impact  on  a  static  transmission  line  were 
described  in  [1],  Chapter  5,  including  the  effects  of  a  sligltly  curved 
tool  surface  on  the  rebound.  The  analysis  of  the  drilling  process  re¬ 
quires  consideration  of  impact  on  a  vibrating  line.  Feng  [13]  has 
solved  the  impact  problem  of  a  spherical  ball  against  a  vibrating 
transmission  line,  with  some  of  these  results  also  reported  in  [l4]. 

The  approach  used  by  Feng  will  be  applied  to  the  present  case  of  a 
cylindrical  tool  having  a  slightly  curved  contact  surface  impacting 
a  vibrating  line. 

The  tip  displacement  of  a  resonant  transducer  (or  transmission 
line)  can  be  expressed  as: 

Uo(t)  =  c^  sin(cot  +  op)  (4-1) 

where 

Op  =  phase  angle 

ci  =  amplitude  of  tip  displacement 
(o  =  angular  frequency  of  resonance 
t  =  time  measured  from  the  stai t  of  tool-line  impact 

The  displacement  3i  of  the  center  of  the  transmission  line  tip 
during  impact  is  a  combination  of:  (i)  the  displacement  due  to  trans¬ 
ducer  tip  vibration  [uo(t)  -  1^(0)],  (ii)  the  displacement  due  to  the 
stress  wave  traveling  down  the  transmission  line,  u1}  and  (iii)  the 
local  deformation,  wL .  The  positive  directions  of  pl5  ux  and  are 
shown  at  the  upper  left  hand  comer  of  Fig.  4-3.  Therefore,  we  have 

=  Ui  +  wx  +  Uq (t )  -  Uo(0)  (4-2) 

Similarly,  the  displacement  at  the  tip  of  the  tool,  02,  is 

02  =  u2  +  w2  (4-3) 
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where  u2  is  the  rigid  body  displacement  of  the  tool,  and  w2  is  the  local 
deformation  of  the  tool  tip.  The  positive  direction  is  assumed  to  be 
the  same  as  px . 


During  the  course  of  impact. 


Pi  =  P2 


v2  -  wx  =  ux  -  u2  +  Uo(t)  -  Uo(0) 


(4-4) 


The  impact  period  is  short  and  the  transmission  line  is  usually 
long  enough  that  the  wave  front  of  the  reflected  stress  wave  will  not 
reach  the  tip  during  this  period.  We  also  assume  that  the  stress  wave, 
due  to  previous  impact,  has  been  damped  out.  Then  the  displacement  u 
can  be  expressed  as 


u(x,t)  =  f(x  -  cDt) 


(4-5) 


Prom  eq.  (4-5),  we  can  obtain 


u  =  -  Cof'(x  -  Cot) 


(4-6) 


F  =  EA^=EAf'(x-  c0t) 


(4-7) 


Eliminating  f' (x  -  cQt)  in  eq.  (4-6)  and  eq.  (4-7),  we  obtain 


u  =  - 


(4-8) 


%  =uo(0,t)  =  -^|  x=0 


(4-9) 


By  considering  the  conservation  of  linear  momentum  of  the  tool,  one 
gets  the  following  equation 


1  rT 
=  VG  +  ±  /  F  dt 
0  m  Jn 


(4-10) 


Substituting  eqs.  (4-9)  and  (4-10)  into  eq.  (4-4),  and  defining  the 
"approach",  a,  as  (wx  -  w2),  we  have 


«  -  vo  +  |  J0  F  «  *eTK^  x= o  -  Ci“  cos(ut  + 


i 
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Herts,  contact  theory  states  that 

F  =  -Ka3/2 


(4-12) 


where  K  is  a  constant  depending  on  elastic  propei’ties  and  geometry  of 
the  contact  surface,  and  is  given  as 


K  =  2ERj2/3(l  -  v2) 


(4-13) 


for  a  flat  transmission  line  tip  where  it  is  assumed  that  tool  and  line 
are  made  of  the  same  material.  Also,  Rp  is  the  radium  of  curvature  of 
the  tool  end,  and  v  is  Poisson's  ratio/  With  eqs.  (4-11),  (4-12)  and 
(4-13),  a  differential  equation  for  a  is  given  by 


_  F  do3 A 

■tp  T  p  Ac0  dt 


d2a 


+  -  a3/2  = 
in 


sin(<jjt  +  q>) 


(4-14) 


A  computer  program  similar  to  Feng's  was  written  to  solve  the 
incident-rebound  relationship  for  tools  of  different  mass  and  various 
curvature  at  the  impact  end.  Figure  4-3  is  a  typical  example  of  the 
impact-rebound  relationship  for  various  tip  vibration  amplitudes.  ALL 
the  tools  are  one  inch  in  diameter  but  with  different  tip  and  end 
geometries.  The  mass  of  the  tools  is  designated  as  me.  In  the  figure 
a  tool  with  curved  surface  of  ten  inches  radius  is  used.  The  tool  has5 
a  one-half  inch  equivalent  length  (me  =  1");  i.e.,  the  tool  has  a  mass 
which  is  equivalent  to  a  flat  ended  rod  of  one  inch  in  diameter  and 
one-half  inch  in  length.  One  observation  that  should  be  noted  is  that 
the  rebound  velocity  is  calculated  in  the  sense  of  average  as  stated 
by  Feng  [13],  because  the  actual  rebound  velocity  varies  as  to  where 
the  tool  comes  into  contact  with  the  transmission  line  tip. 

Figure  4-4  shows  the  results  for  four  different  values  of  me.  The 
main  purpose  of  such  a  plot  is  to  examine  the  effect  of  tool  mass  on 
the  rebound  velocity.  It  is  easy  to  see  that  the  rebound  velocity  drops 
as  the  mass  increases.  The  value  of  rebound  velocity  not  only  depends 

upon  the  mass  of  the  tool  but  also  depends  unon  the  curvature  of  tool 
tip. 


4-3  TRANSDUCER  VIBRATION 

As  pointed  out  previously,  study  in  the  area  of  transducer  charac¬ 
teristics  has  been  extensive,  yet,  much  is  still  unknown  of  transducer 
performance.  This  holds  parti cularly  for  the  response  of  the  transduce: 
during  and  immediately  after  impact  from  the  sonic  tool.  Thorough 
theoretical  as  well  as  experimental  study  is  needed  to  get  a  complete 
understanding  on  this  matter.  However,  it  is  not  the  purpose  at  this 
stage  to  go  into  detailed  research  in  that  area.  Instead,  it  is  simply 
to  study  the  internal  energy  and  tip  vibration  amplitude  of  a  sonic 
transducer  for  the  purpose  of  analyzing  the  sonic  interaction  process 
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AVERAGE  REBOUND  VELOCITY  On./sec) 


INCIDENT  VELOCITY  (in./sec) 


Fig.  4-4  Variation  of  incident-rebound  velocity  relationship 
with  respect  to  equivalent  length  nig  (R  =  5") 


i  i.i.«  PJ^WBppil.-jPf. I'. ,M.,-U'Vr 


and  finding  the  power  transmission  of  the  transducer  system.  The  trans¬ 
ducer  response  to  impact  will  be  modeled  based  on  observed  transducer 
behavior  and  basic  vibrational  principles  governing  solid  materials. 

Ma  [10]  conducted  an  experiment  in  the  study  of  energy  removal  and 
transducer  recovery  during  tool  impact.  From  this  experiment,  an  inter¬ 
esting  phenomenon  can  be  observed.  The  recovery  of  the  tip  vibration 
amplitude  (velocity  as  well  as  displacement)  is  along  an  exponential 
curve  as  shown  within  the  dotted  region  of  Fig.  4-5.  Two  experimental 
results  of  such  removal  and  recovery  curves,  after  data  reduction,  are 
shown  In  Fig.  4-6.  The  ordinate  is  nondimensionalized  with  respect  to 
the  total  amount  of  amplitude  drop  immediately  after  impact  -  Some 
theoretical  background,  in  addition  to  the  previous  experimental  obser¬ 
vations  will  be  given  to  justify  the  resulting  modeling  of  the  recovery 
of  a  sonic  transducer  after  impact. 

For  convenience  of  analysis,  the  discussion  will  start  with  the 
mos'c  simple  and  basic  piezoelectric  resonator,  a  cylindrical  sandwich 
type,  as  shown  in  Fig.  4-7.  Then,  the  results  of  this  analysis  will  be 
expanded  to  the  study  of  the  vibration  and  damping  characteristics  of 
a  more  complicated  sonic  resonator. 

(A)  Forced  Vibration  of  a  Cylindrical 

Sandwich  Type  Resonator 

The  most  simple  and  basic  piezoelectric  resonator  is  the 
cylindrical  sandwich  type  as  shown  in  Fig.  4-7.  Because  of  its  simple 
geometry,  the  analysis  is  straightforward  compared  with  resonators  of 
more  complicated  shape. 

Usually  two  kinds  of  damping  exist  in  a  vibration  system, 
viscous  damping  and  Coulomb  damping.  In  a  sonic  transducer,  viscous 
damping  is  present  within  the  metal  and  ceramic  materials  while  the 
Coulomb  damping  exists  as  interfacial  friction  between  the  metal  and 
ceramic  parts.  The  viscous  damping  force  is  along  the  direction  of 
particle  velocity.  In  the  sonic  transducer,  radial  vibration  is  always 
coupled  with  longitudinal  vibration  even  though  the  first  radial 
resonant  mode  is  usually  a  higher  order  mode  than  the  first  longitudinal 
mode.  In  the  discussion  of  longitudinal  resonance  of  a  transducer,  the 
effect  of  radial  motion  is  usually  neglected  and  the  problem  treated  as 
one -dimensional.  The  direction  of  the  Coulomb  friction  force  in  the 
sonic  resonator  is  radial.  The  effect  of  interfacial  friction  on  longi¬ 
tudinal  vibration  is,  therefore,  indirect.  With  the  above  arguments, 
the  Coulomb  friction  and  viscous  dampirg  force  in  the  radial  direction 
can  be  neglected  if  the  approximation  is  made  to  discuss  longitudinal 
vibration  only. 

In  a  steady  state  resonance  condition,  the  amplitude  of  longi¬ 
tudinal  vibration  is  a  constant.  The  force  that  the  ceramic  applies  to 
the  metal  part  also  reaches  a  steady  state  constant  amplitude.  The 
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Fig.  4-5  Transducer  tip  vibration  under  impact 
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Fig.  4-6  Recovery  of  tip  vibration  amplitude 


Fig.  4-7  Cylindrical  sandwich  type  resonator 


following  discussion  analyzes  the  vibrational  characteristics  of  a 
simplified  model  of  a  cylindrical  sandwich  type  resonator  as  shown  in 
Fig.  4-8.  Two  equal  and  opposite  forces  P0  sin  tot  are  applied  on  a 
cylinder  of  length  Z  and  area  A  at  each  side  of  the  midpoint . 


Viscous  damping  in  metal  is  usually  of  the  Voigt  type  [15], 
[16  J.  The  stress-strain  relation  can  be  expressed  as 

ax  =  Y  ex  +  y’  (4-15) 

where  x  is  the  direction  of  stress  and  strain,  Y  is  the  Young's  modulus 
and  Y'  is  a  damping  factor  [173,  [18],  The  one-dimensional  equation 
governing  the  cylinder  is 


A(x)  =  cross-sectional  area  at  point  x 
u(x,t )  =  longitudinal  displacement 

For  the  simplified  model,  we  can  express  u(x,t)  as  an  infinite 
aeries  of  orthogonal  vibrational  modes  Un(x)  as  follows 

u(x,t)  =£  <pn(t)  un(x)  (U-i7) 

with 

(Jn(x,t)  =  cos(pnx) 

Pn  =  T  (11=1,2,3,....) 

Using  the  principle  of  virtual  work,  a  set  of  differential 
equations  for  (^(t)  can  be  obtained  as 


Mo<Pm  +  Wta  +  Em<Pm  =  -2Po  sin(^)  sin(™)  sin  wt 

where 

Mq  =  Ap.g/2,  rjn  =  AY'm?jt2/2^, 

Em  =  AYmV/2  i  Pm  =  -2P0  sin(~)  sin(— ) 

The  transient  solution  of  eq.  (4-l8)  is 

=  2xq)[-(Tjmt/2Mc)  ]  {Am  sin  qmt  +  Bm  cos  o^t  ) 

where 

jEm  nm  I'1 
"  jMo  " 


(4-18) 


(4-19) 
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The  steady  state  solution  is 


where 


^Pm  =  ^  sin(ut  “  *m) 


=*  (tt  -  «<>“) 


Mou )  +  1m 


i  — 

=  Em  -  Mot^ 


and,  the  steady  state  displacement  us(x,t)  is 


P„ 


^s(x,t)  =  r^~  sin(o)t  -  \| rm)  cos(i^pc) 

n  =  lUZm 


(4-21) 


When  the  exciting  frequency  u  reaches  the  value  f  (Ex /Mq  )  -  (rji/2M§)'j'/i, 
the  coefficient  in  eq.  (4-20)  beccmes  a  maximum,  and 


Pi  _  Pi  _  -2Pn  sin(na/ 

Tjiqx  niqi 


(4-22) 


In  this  condition,  the  system  resonates  at  its  lowest  mode  because  the 
terms  higher  than  one  in  eq.  (4-21)  can  be  neglected  in  comparison  with 
the  first  term.  Therefore, 


us(x,t)  =  cos(^)  sin(wt  -  sin(x) 

Since  Y'  (or  )  is  usually  very  small  for  metals,  we  have 


/IT  «  IT 
“  ~  '  VSb  -  IV  o 


Aj^CoY' 


It  is  known  that  [15 ] 


Q  (factor  of  merit)  = 

'  A F  AI  coY' 


where 


(4-23) 


(4-24) 


(4-25) 


F  =  resonant  frequency, 

AF  =  bandwidth  (between  half-power  points), 
I  =  internal  energy,  and 
AI  =  energy  dissipated  in  a  cycle. 


Equation  (4-24)  can  be  rewritten  as 

Ajt3cnY ' 


Ajt3cnY '  .  L  /  1  V 
^  2  f?  V1  "  V2Q^ 


(4-26) 
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Q  is  usually  very  large  (on  the  order  of  50  to  104),  so  we  can  neglect 
the  second  term  in  the  radical  sign  and  get 


Si 


li 


! 

«  * 


if* 

* 

I 

I 

I 

I 

I 

I 


111!  -  (**-27) 

With  eqs.  (4-19),  (*+-23),  and  (4-27),  we  have  for  the  general  solution 
for  the  vibration  of  the  resonator 


u(x,t)  -  ^2  exPt —  (■+! ) 3  (An  sin(ant)  +  Bn  cos(qnt))  cos(^j^) 

n  =  1  \  £  / 


sin(rta/l) 

AjtJCoY' 


cos^^j  sin(cot  -  ti) 


(4-28) 


In  actual  operation,  the  value  of  P0  depends  on  the  terminal 
voltage  applied  to  the  ceramics  and  prestress  in  the  transducer.  These 
two  variables  are  assumed  to  be  fixed  during  vibration  in  this  study. 
Another  variable  that  affects  the  value  of  Pq  is  the  previously  mentioned 
interfacial  friction.  Since  the  interfacial  friction  force  depends  upon 
the  dynamic  behavior  of  the  transducer,  the  value  of  PD  is  a  function 
of  time.  In  the  previous  analysis,  the  effect  of  interfacial  friction 
is  neglected  and  the  value  of  P0  is  assumed  to  be  constant  to  simplify 
the  problem. 


The  time  factor  Sn,. which  is  equal  to  2Mo/nn  in  t-q.  (4-28), 
is  a  function  of  Y * ;  the  exponential  coefficient  exp(-t/Sn)  represents 
the  effect  of  damping  in  the  system.  We  know  that 


Since  w2  «  Ex/Mq,  we  have 


S 


n 


§1 

n3 


(4-25) 


(4-30) 


where 


(B)  Transducer  Vibration  and  Recovery  of 
Transducer  Tip  Amplitude 


(4-31) 


In  the  previous  section,  the  forced  vibration  of  a  simplified 
cylindrical  sandwich  type  resonator  was  discussed.  It  was  found  that 
the  resonator  will  vibrate  in  a  natural  mode  if  the  exciting  frequency 
is  close  to  the  natural,  mode.  On  the  transient  part  of  the  solution, 
the  time  factor  Sn  depends  on  the  damping  characteristics  of  the 
materials  used.  The  amplitude  of  steady  state  vibration  is  controlled 
by  damping  as  well  as  the  voltage  appl  ied  to  the  ceramics . 

93 


rp.n„  .  The  results  obtained  for  a  simple  cylindrical  sandwich  type 
resonator  wall  now  be  expanded  to  a  more  complicated  resonator  suchas 
_7  4P"i1  transducer  wlth  body.  resonant  horn  and  transmission  line. 

differ entXmodP  the,  type  of  solution  will  exist,  except  with 

Afferent  mode  shapes,  tune  factor  and  steady  state  amplitude :  i . e . 

frpn6  Wl11  p6  a  !tea^y  state  solution  with  a  lowest  natural  mode  if  "the 
quency  o.  excitation  is  at  that  mode,  and  an  amplitude  which  increases 

solSioneS1^'in?ltafe  inCTeaSe?i  and  decreases.  The  transient 

?•  ?  f?  infinite  series  combination  of  all  the  natural  modes  with 
exponential  tune  functions.  The  time  factors  depends  on  all  the  damping 
characteristics  ot  transducer.  In  written  form,  we  have  the  following?6 


u(x,t)  =  Afi(x)  sin (wt  -  ti ) 


00 

n!C1  exP(~Vsn)  (An  sinCq^t)  +  Bn  cos(qnt))  fn(x) 


(4-32) 


where  fn(x)  is  the  nth  mode  shape.  The  lowest  mode  of  vibration  is 
*1 W « 


In  general,  since  the  mass  of  the  tool  is  rather  small  .  the 
nergy  loss  of  the  transducer  for  each  impact  is  small  compared  to  the 
total  stored  energy.  If,  for  example,  a  tool  of  mass  0.003439  slug 

5,;.V  a  transducer,  which  has  a  tip  vibration 

amplitude  of  1.15  mils,  and  the  incident  and  rebound  velocity  of  the 

"e  equai  &-6  to-/sec  and  111  in. /sec,  respectively^  the 
nergy  loss  of  the  transducer  is  equal  to  the  energy  gain  of  the  tool 

^31^17^5 'i^  lb’’1ThWhlle  thf  ±nternal  energy  is  approximately 
eqpal  to  17.29  in. -lb.  The  percentage  of  energy  loss  is  about  7.54$. 

n~  ^Hhth+-enerSr  l0SS  °f  thG  transducer  is  small  per  impact,  the  mode 

tl0n  ^er  XmpaCt  is  >  ssentdaAly  fi(ac);  i.e.,  the  value  of  A_ 
and  Bn  are  small  compared  to  A.  “ 


According  to  previous  analysis, 


Sn  =  Si/n* 


(4-33) 


the  hdShe>  -order  terns  in  the  transient  solution 
damp  out  much  faster  thin  the  lowest  mode.  Therefore,  another  assump- 

°a/  S,!Jf  aeTe  t0  n-6lect  the  effect  of  higher  order  terms  in 
eq.  ^4-32;,  giving 


u(x,t)  «  f  1  (x)  (A  sin  (u't  -  \|fx)  -  exp(-t/S1 )  [Ax  sin^t)  +  Bx  cos(q1t)]) 

(4-34) 


eqUation  rePresents  approximately  a  harmonic  vibration  with 
p-adually  varying  amplitude.  In  another  words,  the  recovery  of  the 
transducer  from  energy  loss  due  to  tool-line  impact  is  exponential. 

The  above  analysis  has  been  confirmed  by  the  experimental  results 

aarlJr,in  ^Jis  section.  Therefore,  we  may  model  the  recovery  of 
transducer  tip  vibration  amplitude  c(t)  as  follows, 
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(4-34) 


c(t)  =  C  -  A  exp(-t/S) 

where  C  represents  the  maximum  amplitude  at  the  no-load  condition,  A 
is  the  amplitude  drop  of  the  transducer  after  impact  and  is  dependent 
on  the  energy  drain  by  the  tool.  Immediately  after  impact,  t  =  0,  the 
tip  amplitude  is 


c(0)  =  C  -  A 

As  time  goes  on,  c(t)  increases  in  an  exponential  order  with  a  time 
factor  S  which  depends  on  the  total  damping  characteristics  of  a  sonic 
transducer,  i.e.,  the  damping  in  ceramic  as  well  as  in  metal  is  also 
considered. 

(C)  Internal  Energy  and  Tip  Vibration  Amplitude 

The  ability  of  a  sonic  transducer  to  do  work  is  dependent  on 
the  amount  of  energy  in  the  transducer.  It  is  known  that  internal 
stored  energy  in  a  longitudinally  vibrating  resonator  is  in  the  form  of 
strain  and  kinetic  energy.  The  amount  of  energy  is  proportional  to  the 
square  of  maximum  displacement  (or  amplitude)  at  any  point.  In  Section 
4-2,  the  tip  vibration  amplitude  was  used  in  discussing  the  tool  impact- 
rebound  relationship  on  a  vibrating  transmission  line.  The  amount  of 
tip  vibration  amplitude  will  again  be  used  as  the  measurement  of  tie 
internal  energy  of  a  transducer. 

It  is  not  difficult  to  find  the  stored  energy  in  terms  of  tip 
amplitude  in  the  resonant  horn  and  transmission  line  once  the  mode  of 
vibration  is  known.  It  is  more  difficult  to  find  the  energy  which  is 
stored  in  the  body  of  the  transducer.  This  difficulty  arises  from  the 
complicated  structure  of  this  part.  An  approximate  value  is  obtained 
by  assuming  the  approximate  mode  shape  in  the  body  and  neglecting  the 
energy  stored  in  the  ceramics.  The  proportional  constant,  K,  for  the 
P-11  transducer  with  a  10-inch  transmission  line  in  the  equation 

I  =  K  c2  (4-36) 

is  found  approximately  to  be  equal  to  11.2  in.-lb/(mil)2.  In  this 
equation  I  is  the  internal  energy,  and  c  is  the  tip  vibration  amplitude. 
The  energy  stored  in  the  body  of  a  stepped  transducer  has  been  calcu¬ 
lated  and  is  much  smaller  than  in  the  resonant  horn  and  transmission 
line.  Hence  the  approximate  value  of  K  thus  obtained  is  expected  to  be 
close  to  the  true  value.  The  reason  is  that  the  stored  energy  in  the 
resonant  horn  and  transmission  line  can  be  calculated  exactly  while 
approximation  is  made  only  in  calculating  the  stored  energy  in  the  body. 
If  the  energy  in  the  body  is  small,  the  approximation  just  made  would 
not  affect  the  total  stored  energy  drastically.  The  value  of  K  can  be 
considered  as  a  measurement  of  the  energy  storage  capacity  of  a  sonic 
transducer.  The  larger  the  value  of  K  the  more  the  energy  can  be 
stored. 
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The  value  of  the  tip  displacement  is  important  in  the  impact 
coupling  process  and  difficult  to  obtain  from  purely  theoretical  analysis. 
The  value  is  also  dependent  on  the  terminal  voltage  of  the  electric 
supply  and  temperature  of  the  transducer  system.  An  experiment  was 
conducted  to  find  the  effect  of  temperature  and  voltage  on  the  tip 
amplitude  of  an  unloaded  P-11  transducer  with  a  ten-inch  transmission 
line.  It  was  found  that  the  amplitude  of  vibration  increases  as  the 
applied  voltage  increases  as  shown  in  Fig.  4-9.  Figure  4-10  shows  the 
peak-to-peak  amplitude  and  power  supply  as  a  function  of  frequency  of 
vibration.  Tbs  voltage  is  apprOxiUHttely  1,200  V. 

To  find  the  effect  of  temperature  change  on  the  vibration 
characteristics  of  the  P-1 1  transducer  with  a  10-inch  transmission  line, 
a  set  of  curves  similar  to  Fig.  4-10  but  at  different  temperatures  is 
shown  in  Fig.  4-11.  The  resonant  frequency,  amplitude  and  input  power 
varies  as  the  temperature  changes.  At  room  temperature  the  amplitude 
is  largest  while  the  power  required  to  maintain  steady  state  vibration 
is  least.  As  time  goes  on,  the  heat  accumulates  in  the  transducer  and 
causes  the  temperature  to  rise.  The  amplitude  drops  drastically  while 
more  energy  is  required  to  maintain  resonance,  which  means  that  more 
energy  is  consumed  as  heat  and  acoustic  noise  at  high  temperature. 

Unde mbtedly  the  temperature  has  a  large  effect  on  the  properties  of  the 
piezoelectric  material  since  the  properties  of  metal  are  fairly  stable 
for  such  smaJI  ranges  of  temperature  change.  Another  important  fact  is 
the  drop  in  resonant  frequency  of  the  system  as  the  temperature  rises . 

This  is  not  particularly  desirable  if  a  single  frequency  power  supply 
is  to  be  used. 


4-4  DRILLING  ANALYSIS 

As  stated  previously,  the  sonic  rock  cutting  process  is  achieved 
y  using  an  intermittent  tool  between  the  transmission  line  tip  and  the 
work  surface.  The  tool  acquires  energy  from  impacting  with  the  vibrating 
transmission  line  tip;  it  then  impacts  against  the  rock  and  rebounds 
fran  the  surface.  During  impact  between  tool  and  rock,  the  tool  loses 
part  of  its  energy  and  causes  some  fracture  of  the  rock. 

The  resonant  frequency  of  a  sonic  transducer  can  be  affected  by 
many  factors,  such  as  voltage,  temperature,  and  static  force.  The 
direct  result  of  these  detuning  effects  is  to  change  the  tip  vibration 
amplitude.  Since  the  sonic  transducer  is  a  high  Q  system,  any  small 
amount  of  frequency  change  can  affect  the  tip  vibration  amplitude 
drastically.  This,  in  turn,  reduces  the  capability  of  the  transducer 
to  do  work.  Future  efforts  should  be  given  to  making  the  transducer 
resonant  at  all  times  during  operation.  In  other  words,  the  systan 
interaction  only  on  the  assumed  conditions  that  the  transducer  is  always 
in  a  resonant  state,  and  that  the  tip  vibration  amplitude  is  a  constant 
during  steady  state,  no  load  operation  will  be  discussed. 
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Fig.  4-10  Variations  of  amplitude  and  input  power  as  a 
function  of  driving  frequency 
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(A)  Internal  Energy  During  Steady  State  Operation 

The  recovery  of  the  transducer  after  tool-line  impact  has 
boen  discussed.  It  was  found  that  the  tip  vibrai  Ion  amplitude  will 
recover  in  an  exponential  manner.  If  no  further  impacts  occur,  the  tip 
vibration  amplitude  will  eventually  reach  the  steady  state  no-load 
amplitude .  In  actual  operation,  the  tool  will  bounce  back  and  forth  in 
the  gap  between  the  transmission  tine  tip  and  load.  The  impacts 
tool  against  the  tip  can  be  at  hundreds  and  even  thousands  of  times  per 
second.  The  tip  amplitude  has  no  chance  of  fully  recovering  back  to 
its  steady  state  no-load  condition  ttoausi.  the  flesci'.l  ^ 

the  first  one  in  a  short  period.  Hence,  the  internal  energy  during  the 
steady  state,  no-load  vibration  represents  a  maximum  value. 

As  discussed  in  Section  4-3 ,  the  tip  vibration  amplitude  is 
assumed  to  be  c(t)  =  C  -  A  exp(-t/S),  The  corresponding  internal  energy, 
according  to  eq.  (4-36)  is 

l(t)  =  K[C  -  A  exp(-t/S)]2  (4-37) 

Figures  4-12  and  4-13  represent  the  variation  of  tip  vibration 
amplitude  and  internal  energy  during  operation.  The  tip  amplitude  and 
internal  energy  at  the  steady  state  unloaded  condition  are  C  and  Iq, 
respectively.  When  a  load  is  applied  to  the  system,  the  tool  acquires 
energy  from  impacting  with  the  transmission  line  tip.  As  time  goes  on, 
the  amplitude  and  internal  energy  recovers  due  to  the  power  supply 
input.  Before  the  transducer  recovers  to  its  full  energy,  the  tool 
bounces  back  from  the  rock  surface,  and  again  drains  energy  from  the 
transducer.  This  process  continues  until  a  steady  state  condition  is 
reached.  The  dashed  lines  represent  the  recovery  curves  of  the  trans¬ 
ducer  if  no  further  drain  of  energy  from  the  tool  occurs. 

If  eq.  (4-37)  is  expanded,  one  obtains 
I(t)  =  K{C  -  A  +  A[1  -  exp(  t/S)]}2 

=  K(C  -  A)2  +  2K(C  -  A)A[1  -  exp(-t/S ) ]  +  KA2[1  -  exp(-t/S)]2 

(4-38) 

The  first  term  represents  the  internal  energy  right  after  impact.  The 
second  and  third  tenns  represent  the  recovery.  It  is  easy  to  see  that, 
in  the  second  term,  the  recovery  is  also  dependent  on  C  and  K  as  well 
as  the  time  factor  S.  The  larger  the  value  of  C  and  K  the  faster  the 
energy  recovery. 

In  the  steady  state  operating  condition,  the  energy  recovered 
equals  the  energy  drained  from  the  transducer,  assuming  Ti  is  the 
impact  interval.  Therefore, 

li  =  X (Ti )  =  K[C  -  A  exp(-t/S)]2  =  K(C  -  A)2  +  A  (4-39) 
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where  A  is  the  energy  that  the  tool  acquired  from  the  transducer;  i.e., 


A  =  i  m (Vi  -  vf )  +  % 


(4-40) 


with  Ejj  being  the  energy  loss  due  to  impact.  Vi.  and  V2  are  the  velocities 
of  the  tool  before  and  after  impact,  respectively.  If  the  transmission 
line  tip  and  tool  are  made  of  hard  material,  not  much  plastic  deforma¬ 
tion  should  occur.  Therefore,  is  considered  to  be  small  compared  to 
m^V|  -  vf)/2  and  can  be  neglected.  Thus, 


A  =  m(vf  -  vf)/2 


(4-41) 


luring  this  steady  state  condition,  the  internal  energy  is  at 
a  specific  value  before  impact  which  implies  that  the  tip  vibration 
amplitude  is  also  at  a  specific  value.  Now,  there  aie  several  questions 
to  be  answered:  What  is  the  value  of  the  interna]  energy  and  tip  vibra¬ 
tion  amplitude?  What  are  the  values  of  impact  and  rebound  velocity  of 
the  tool  against  the  transmission  line  tip  and  rock?  These  problems 
will  be  answered  as  the  discussion  procedes. 

In  actual  operation,  the  impact  and  rebound  velocity  of  the 
tool  is  somewhat  random  and  so  is  the  internal  energy  and  tip  vibration 
amplitude;  but,  in  the  sense  of  averages,  they  are  assumed  to  be  deter¬ 
ministic  in  the  above  analysis. 

(B)  Impact  Coupling 

The  general  idea  of  internal  energy  during  the  steady  state 
operating  condition  was  discussed  in  the  preceding  section.  Now  the 
tool  behavior  during  this  condition  will  be  studied.  The  tool-rock  and 
tool -line  phases  of  the  impact  coupling  process  have  been  studied  exten¬ 
sively.  The  present  task  is  to  combine  them  for  the  study  of  this  coupling 
process.  Feng  [13]  studied  the  impact  coupling  of  a  spherical  ball 
between  a  static  and  a  vibrating  transmission  line.  The  present  analysis 
follows  his  method  of  approach. 

In  [1],  Part  II,  Chapter  4,  the  incident-rebound  velocity  rela¬ 
tionship  of  a  tool  impacting  against  rock  was  represented  by  a  best-fit 
straight  line  while  in  Section  4-2  of  this  chapter,  the  incident- 
rebound  velocity  relationship  of  a  tool  impacting  against  a  vibrating 
transmission  line  was  obtained.  If  these  two  relationships  are  plotted 
on  the  same  figure  with  the  notation  of  incident  and  rebound  velocity 
as  shown  in  the  lower-right-hand  comer  of  Fig.  4-l4,  we  obtain  an 
example  figure  as  shown.  In  this  case,  a  tool  with  l/2-inch  equivalent 
length  (me  =  -=r")  and  10-inch  radius  of  curvature  of  the  tool-tip  contact 
surface  is  used.  The  rebound-incident  velocity  ratio  and  tip  amplitude 
ci  are  assumed  to  be  equal  to  0.6  and  1.0,  respectively.  The  steady 
state  velocities  happen  to  be  at  the  point  of  intersection  of  these  two 
relationships,  i.e.,  the  point  W(Vi  =  57-5  in. /sec  and  V2  -  96.0  in. /sec). 
To  illustrate  this,  let  the  tool  start  at  a  velocity  VG1  and  impact 
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against  the  vjbrrting  line  (point  Mi).  The  rebound  velocity  V2  from 
the  line  becomes  the  incident  velocity  V2  for  tool-rock  impact  (point 
Mp).  With  this  rebound  velocity  VS  from  rock,  the  tool  impacts  again 
against  the  transmission  line  tip  (point  M3).  This  process  continues 
until  the  point  of  intersection  W  is  reached.  The  same  result  will  be 
obtained  if  we  start  at  a  lower  velocity  V31 .  The  above  is  true  as 
long  as  the  tip  violation  amplitude  remains  unchanged.  the  Amplitude 

c^  changes  so  does  the  point  W  representing  the  steady  state  velocities. 
Therefore  ,  we  can  say  that  steady  state  velocities  VL  and  V2  are  func¬ 
tions  of  since  the  tool-rock  impact  relation  and  the  geometry  of  the 
tool  are  assumed  to  remain  unchanged  during  operation. 

It  is  interesting  to  note  that  the  situation  may  arise  where  the 
slope  of  the  rebound-incident  velocity  relationship  for  tool-rock  impact 
(slope  =  Vj./V2)  is  so  small  that  no  point  of  intersection  W  can  be  ob¬ 
tained.  This  means  that  no  steady  state  velocities  can  be  reached  a.id  the 
problem  then  becomes  more  complicated.  It  is  believed  that,  in  some  in¬ 
stances,  the  tool  tip  is  in  close  contact  with  the  rock  while  the  trans- 
mission  line  hits  the  tool  at  the  end;  i.e.,  the  tool-rock  and  tool-line 
impacts  happen  at  more  or  less  the  same  time.  The  s'  ope  of  the  rebound- 
incident  velocity  relationship  for  tool-rock  impact  (Vi/V2)  depends  on  the 
tip  geometry  of  the  tool  and  the  properties  of  the  rock.  Usually  harder 
rock  yields  higher  values  of  Vi/Va.  in  the  remainder  of  this  analysis, 
only  the  case  where  steady  state  impact  coupling  occurred  was  considered. 
Therefore,  our  study  was  restricted  to  hard  rock  drilling  in  which  inter¬ 
section  point  W  exists  for  a  specific  tool  geometry. 

Another  interesting  observation  in  previous  analysis  of  impact 
coupling  is  that  the  steady  state  velocities  do  not  change  so  long  as 
the  tip  vibration  amplitude  remains  unchanged;  the  static  force  Fs 
applied  by  the  support  structure  of  the  transducer  plays  no  direct  role. 
It  seems  to  have  no  effect  on  the  impact  process.  Actually  the  static 
force  affects  the  steady  state  velocities  via  affecting  the  value  of 
steady  state  tip  vibration  amplitude  This  will  dear  as  the 
discussion  of  impact  coupling  contimies. 

The  discussion  of  static  force  is  started  by  looking  at  the 
dyxwjtifs  cf  a  *Cnir  truieducvr  eofiaiiterei  by  UralT  [ly]  (Fig*  4-15). 

The  rigid  body  motion  is  governed  by 

M  x  =  Fs  -  f(t)  (4-42) 

where  M  is  the  mass  of  the  transducer,  f(t)  is  the  time  varying  force 
at  the  tip  of  the  transducer  due  to  the  impact  of  a  sonic  tool,  and  Fg 
is  the  static  force  applied  by  the  support  structure.  Force  f(t)  is  a 
combination  of  a  sequence  of  impact  forces,  each  impact  force  is  high 
in  magnitude  but  short  in  duration.  It  usually  is  less  than  50  nsec. 
Actually,  the  impact  forces  would  be  of  various  amplitudes  and  intervals, 
but  thinking  in  terms  of  averages  it  is  justifiable  to  consider  equal 
amplitudes  end  spacings.  The  interval  of  impacts  is  quite  long  compared 
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Transducer 


Mass=M 


F-ig.  4-15  Trandsucer  acted  upon  by  static  force  F_  and 
impact  force  f(t) 


with  each  impact  duration  in  sonic  processes.  For  example,  if  the  fre¬ 
quency  of  impact  is  1000  impacts/sec,  the  time  interval  is  in  the  order 
of  1000  nsec . 


Integrating  eq.  (4-42)  gives 

M(x  -  Xq)  =  Fst  -  f (t)dt 


(4-43) 


Since  it  was  assumed  that  each  impact  force  is  the  same  and  the  duration 
is  short,  it  is  justifiable  to  represent  the  second  term  as  follows 


C  =  im  £  H(t  -  tn) 

u  n  =  1 


(4-44) 


where  N  is  the  number  of  impacts  between  the  time  interval  0  and  t  and 
H  is  a  Heaviside  function.  Now,  rewriting  eq.  (4-43) 


M(*  -  Xq)  =  Fst  -  Im  T  H(t  -  t,,) 


(4-45) 


s 


This  equation  shows  that  the  velocity  of  a  sonic  transducer  has  a  saw¬ 
tooth  form  as  shown  in  Fig.  4-l6.  Thus,  the  constant  force  Fs  causes 
a  linear  increase  in  velocity,  but  this  is  periodically  offset  by 
sudden  velocity  decreases  due  to  the  impact  pulses.  If  the  proper 
balance  holds  between  the  static  force  and  the  impulse  Im  and  inpulse 
interval  T^,  the  increase  in  velocity  is  completely  offset  by  the  drop 
caused  by  impact,  so  that 


FsTi  =  In 


(4-46) 


This  is  a  statement  of  impulse  balance,  so  that  the  average  net  change 
in  momentum  of  the  transducer  is  zero. 


According  to  the  definition  of  impulse, 

Ijn  =  m(V!  +  V2) 

Therefore,  fran  eq.  (4-46),  we  have 


T<  = 


m(Vi  +  Vs 


(4-47) 


(4-48) 


The  value  of  Ti  consists  of  impact  times  and  the  travelling  times  of 
the  tool  between  the  transmission  line  tip  and  the  rock  surface.  T  e 
frequency  of  impact  coupling  is  just  equal  to  the  reciprocal  of  T±;  or, 


f  =  l/Ti  = 


m(Vx  + 


(4-49) 


Velocity  increase  due  to  Velocity  decrease 
static  force  /  due  to  impact 


TIME 


Fig.  4-l6  Rigid  body  motion  of  a  transducer  acted  upon  by  static 
force  and  impact 
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As  stated  previously,  steady  state  velocities,  Vx  and  V2,  are  functions 
o  tip  vibration  amplitude  c^.  Thus,  the  relationship  between  impact 
frequency  and  the  tip  vibration  amplitude  are  related  to  each  other 
according  to  eq.  (4-49).  Generally,  the  increase  or  decrease  of  Vx  and 
V2  follows  the  trend  of  the  value  of  c^ .  When  we  increase  the  impact 
frequency,  the  amplitude  drops,  therefore  reducing  the  internal  stored 
energy  of  the  transducer.  If  the  static  force  is  kept  constant  during 
operation,  it  is  easy  to  see  from  eq.  (4-48)  that  T,-  is  a  function 
of  c*.  1 

(c)  Power  Transmission 

Now  again  consider  the  internal  energy  of  the  transducer,  and 
rewrite  eqs.  (4-35),  (4-39),  (4-4l)  and  (4-47)  as  follows: 

Internal  energy  just  before  impact  is 

Ii  =  K(C  -  A)2  +  A  =  K[C  -  A  exp(-Tj^/S)]2  (4-50) 

The  energy  drain  in  each  impact  is 

A  =  m(Vf  -  Vf)/2  (4-51) 

Tip  vibration  amplitude  was  given  previously  by  c(t)  =  C  -  A  exp(-t/S) 
while  the  interval  of  impulse  is  '  ’ 


TjL  =  1/f  =  5^V1  +  Vp) 

Fg 


(4-52) 


Eliminating  A  is  eqs.  (4-39)  and  (4-4l),  and  solving  for  A,  we  obtain 

C  emp(Tj/S)  t  ^[1  -  exp(-Tj/S ) j2)  -  (m/2K)(Va-V1 )  [1  -  exp(-Tj/S)] 


A  = 


1  -  exp(-Tj_/S) 


(4-53) 


Substituting  (4-53)  is  eq. 
c  (t  .  )  =  C  exP(Tj/S)  t  ^ 


(4-35)  and  letting  t  =  Tj_,  we  get 

-  (m/2K)  (Vg-Vx)  [1  +  exp(Tj/S  )]/[-!  +  exp(T-j/s)l 
1  +  exp  (T  j_/'S )  ~  ~ 


(4-54) 


To  decide  whether  the  positive  or  negative  sign  in  the  above  equation 
should  be  chosen,  consider  a  transducer  with  large  K,  i.e.,  a  transducer 
whose  energy  storage  capacity  is  very  large.  Then,  by  neglecting  the 
second  term  in  the  radical  sign,  we  have 


ci  =  c(Ti)  ~  C[exp(Ti/s)  ±  l]/[exp(Ti/S)  +  1]  (4-55) 

If  we  choose  the  positive  sign,  we  have 
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Ci  =  c(Ti)  =  C 

This  means  that,  for  a  transducer  with  large  eneigy  storage  capacity, 
the  steady  state  vibration  amplitude  before  impact  is  close  to  the 
unloaded  vibration  amplitude  C.  This  is  consistent  with  the  statement 
made  in  Section  (A)  of  this  chapter  that  the  larger  the  value  of  K  the 
faster  the  energy  recovery.  If  we  choose  the  negative  sign,  we  have 

cCl’i)  C[exp (T-l/S  -  l)]/[exp(Ti/s)  +  l] 

This  has  no  special  physical  meaning.  Therefore,  we  choose  the  positive 
sign  in  eq.  (4-53)  and  (4-54),  i.e., 

ci  =  c(Ti) 

C  exp(Ti/s)  +  v/c2  -  (m/2K ) (Vjf-Vf ) [1  +  exp^/S ) ]/[-!  +  expC^/S)] 

(4-56; 

As  discussed  at  the  end  of  Section  (B),  the  static  force  Fs  is  kept 
constant,  and  Vx ,  V2,  and  T^  are  functions  of  cp  [or  c(Ti)].  Then,  the 
right  hand  side  of  eq.  (4-5b)  is  a  function  of  cp  only.  Therefore,  the 
value  of  cp  can  be  solved  numerically  if  not  analytically,  since  no 
mathematical  relationship  between  Vx,  V2  and  c±  exists. 

The  questions  stated  at  the  end  of  Section  (A)  become  clear 
once  the  values  Ci,  or  c(Ti),  is  obtained.  The  values  of  Vx  and  V2  can 
be  found  according  to  Section  (B).  The  value  of  Tp  can  be  obtained  by 
using  eq.  (4-52)  with  the  values  of  Vx  and  V2  just  found.  The  internal 
energy  just  before  impact  can  also  be  evaluated  according  to  eq.  (4-36). 

With  all  the  information  obtained  previously,  the  power  trans¬ 
mission  into  rock  can  be  calculated  by  considering  the  energy  loss  of 
tool  to  rock  in  each  impact  and  the  frequency  of  the  sonic  process . 

The  energy  absorbed  per  impact,  A,  is 

A  =  m(vf  -  Vf )/2  (4-57) 

Therefore,  the  power  transmission  P  is  the  product  of  A  and  f,  i.e., 

P  =  Af  (4-58) 

Substituting  f  of  eq.  (4-53)  in  the  above  equation,  we  have 

P  =  F  (V2  -  Vl)/2  (4-59) 

With  the  value  of  ci  known,  the  value  of  V2  -  Vx  (or  P)  can  be  obtained. 
Figure  4-17  shows  the  values  of  V2  -  Vx  for  R  =  10"  and  various  values 
of  equivalent  length  of  the  tool,  me.  The  ratios  of  the  rebound-incident 
velocity  relationship  in  tool-rock  impact  used  here  are  equal  to  0.7  and 

0.6. 


1  +  exp(Tx/S] 
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(D)  Dr  tiling  Rate  of  a  Sonic  Transducer 

As  we  mentioned  previously  in  discussing  the  specific  energy 
(es)  of  rock  removal,  the  time  required  to  remove  a  certain  volume  of 
rock  is  important  in  the  discussion  of  economy  of  a  drilling  system. 
Therefore,  drilling  rate  is  the  focusing  point  for  any  kind  of  rock 
drilling  system.  With  the  power  transmission  (P)  and  specific  energy 
known  for  a  sonic  transducer  system  and  rock,  the  drilling  rate  (R) 
can  be  found  very  easily  as  follows 


R  =  P/e£ 


(4-60) 


With  all  the  studies  on  the  sonic  rock  drilling  being  completed,  an 
example  problem  using  the  information  presented  is  given.  This  will 
aid  in  understanding  the  previous  studies  and  also  give  a  general  idea 
of  a  practical  situation  in  a  sonic  rock  drilling  process. 

(E)  Numerical  Example 

The  information  with  regard  to  the  transducer,  static  force 
applied,  the  tool  and  the  rock  are  given  as  follows: 


TRANSDUCER: 


P-11  transducer  with  1-inch  diameter  and  10-inch 
long  transmission  line 
Resonant  frequency  =  9750  Hz 
Applied  terminal  voltage  =  1200  V 


STATIC  FORCE:  100  lb. 


TOOL: 


ROCK: 


Dimensions  -  Type  F  tool 
Radius  of  curvature  of  the  tool  end  =  10" 
Equivalent  length  mg  =  \ 


in 


Granite 

Vi/V2  =0.6 

(Laboratory  experiment  shows  that  the  value  of  rebound- 
incident  ratio  is  between  0.54  to  0.60  for  Type  F  tool 
impacting  against  granite.  0.6  is  used  here  for 
analysis . ) 

Specific  energy  es  -  Fig.  4-2 


The  internal  energy  I  according  to  eq.  (4-36)  is  I  =  Kc2 
where  K  is  equal  to  11.2  in.-lb/(mil)2  for  P-11  transducer  with  a  10  in. 
transmission  line  [following  eq.  (4-36)]. 

Angular  frequency  oj  =  2nf  =  6l,260  rad/sec. 

No-load  amplitude  C  =  3.4/2  =  1.7  mils  according  to  Fig.  4-9. 
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Time  factor  S  is  dependent  upon  the  total  damping  characteristics  of  the 
transducer.  Assuming,  that  eq.  (4-31)  still  holds  after  the  modeling  of 
transducer  tip  vibration  amplitude  recovery  is  made,  then  S  =  2Q/oj  where 
Q  can  be  obtained  by  using  eq.  (4-25);  i.e.,  Q  =  F/aF,  and  data  from 
^8*  4-10.  The  value  of  Q  is  found  to  be  equal  to  76.5  for  the  terminal 
voltage  of  1200  V  and  temperature  of  75° F.  Therefore, 

S  =  2  x  76.5/61,260  =  2.498  m  sec. 

Tf  =  m(Vi  +  V2)/Fs  =  0.0055027  V2 
Ti/S  =  0.001147536  V2 


(m/2K)(v|  -  Vf)[l  +  exp(Ti/S)] 


-1  +  exr,l 


1  +  exp(Ti/S) 


=  8.1885  x  10-6  v| 

-1  +  exp  (Ti/S ) 


(mil)' 


With  the  above  information,  the  curves  for  the  left-hand  side  and  right- 
hand  side  of  eq.  (4-60)  are  <p  =  c(Tj_), 


C  exp(Tj/s)  +  JC2  -  (m/2K)(vf  -  yf)[l  +  exp(Tj/s) ]/[-!  +  exp(Ti75)T 


exp (Ti/S)  +  1  ' 

and  are  plotted  in  Fig.  4-18.  The  solution  of  eq.  (56)  is  the  inter¬ 
section  of  the  above  two  curves,  i.e.. 


Ci  =  c(Ti)  =  1.15  mils 


From  Fig.  4-17, 


Vx  -  V2  =  44.4  in. /sec 
Vi  =  66.6  in. /sec 
V2  =  111.0  in. /sec 
Ti  =  0.504  m  sec 

Ii  =  11.2  x  (1.15)2  =  17-25  in. -lb 
A  =  \  m(v|  -  Vf)  =  1.124  in. -lb 
A/Ii  =  1.124/17.25  =  7.54 % 


Power  transmission  P  is 


P  =  \  Fs(V2  -  V3.)  =  2,220  in. -lb/sec  =  11,100  ft-lb/min 

=  250.8 


|c-=  1.15  mils 


PLITUDE,  c *  (mils) 


(graphical  method) 


puggs1 


t. 


We  do  not  have  the  specific  energy  chart  for  me  =  j  in.  but 
do  have  it  for  mg  =  1.08  in.  (Fig.  4-2).  To  find  the  specific  energy 
for  me  -  0.5  in.  and  V2  =  111.0  in. /sec,  assume  that  equal  incident 
energy  of  the  tool  would  cause  equal  amount  of  damage  to  the  rock  if 
t.'p  geometries  are  the  same  for  tools  of  different  mass  (or  me);  then 
the  corresponding  incident  velocity  V2  of  the  tool  with  can  be  found 
an  follows 


2  me  v|  =  -| 


1  ™  1  V  ’  < 

2  me  *2 


me  =  \  in.  and  me'  =  1.08  in. 

V2’  =  75-5  in. /sec 

The  value  of  es  corresponding  to  V2  =  75*5  in. /sec  and  me  =  1.08  in. 
is  11,000  ft -lb/in. 3  according  to  Fig.  4-2.  Therefore,  the  drilling 
rate  in  volume  per  minute  is 

Ry  =  P/es  -=  11100 /llOOO  =  1.01  in.3/min 

and  the  drilling  rate  in  in. /min  is 

RL  =  Ry/A  =  1.01 /«(|)2  =  1.29  in. /min 

Drilling  rates  and  power  transfer  are  also  calculated  for  the 
case  when  the  terminal  voltage  and  static  force  applied  to  the  trans¬ 
ducer  are  equal  to  2000  V  and  200  lb,  respectively.  An  estimated  value 
of  Q,  is  used  (Q  =  60)  in  this  calculation  since  no  experimental  result 
is  available  for  a  transducer  driven  at  2000  V.  The  results  are  listed 
in  the  following  table  along  with  the  previous  example. 

Table  4-1  -  Sonic  Drilling  of  Granite 


1200  V 
100  lb. 

2000  V 
200  lb. 


Drilling 

Drilling 

Power 

transfer 

Ti> 

rate  Rl 

rate  Rv 

in. /min 

in. 3 /min 

W 

ft-lb/min 

m  sec . 

1.29 

1.01 

250.8 

11,100 

0.504 

8.58 

6.75 

819.0 

36,250 

0.416 
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The  drilling  rates  of  1 
conventional  percussive 
drilling  rate  with  twu 
apparatus  is  within  the 
piston  energy  consumed 
ft-lb/min.  Considering 
sonic  transducer  system 
low  range  of  percussive 


29  and  8.58  in. /min  obtained  are  lower  than 
drilling.  Bruce  [20],  in  a  study,  found  the 
special  bits  (four -wing)  and  percussive  drilling 
range  of  5  to  40  in. /min  while  the  maximum 
is  in  the  range  of  46,000  ft-lb/min  to  162,000 
the  power  consumed  and  the  drilling  rate  oi  a 
,  it  is  easy  to  see  that  the  system  is  in  the 
drilling  systems. 


2. 
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5.  TRANSDUCER  ANALYSIS 


Although  sonic  and  ultrasonic  power  transducers  axe  widely  used, 
the  design  of  transducers  is  still  more  an  art  than  a  science.  Con¬ 
cepts  such  as  the  half  wavelength  design  and  the  full  wavelength  design 
are  the  dominant  approaches  in  industry.  Theoretical  design  based  on 
more  exact  theories  is  lacking,  especially  for  the  high  power  applica¬ 
tions.  Possibly  this  may  be  because  of  the  complicated  geometry  of  an 
actual  transducer .  Here,  attempts  will  be  made  to  apply  the  more  "exact" 
theory  of  linear  piezoelectricity  to  predict  the  resonant  frequency  of 
a  transducer.  However,  the  lack  of  an  applicable  short  piezoelectric 
rod  theory  leaves  the  long  rod  and  large  plate  theories  as  the  only 
remaining  a3  +  ematives.  Thus,  these  two  models  will  be  followed 
throughout  the  report. 

Additionally  and,  possibly  more  importantly,  is  the  objective  of 
determining  the  causes  for  shifts  in  the  resonant  frequency  of  trans¬ 
ducers  operating  in  different  environments.  In  the  Sonic  Power  Labora¬ 
tory  of  The  Ohio  State  University,  the  resonant  frequency  of  a  transducer- 
in  operation  has  been  observed  to  shift  with  different  levels  of  pre¬ 
stress  on  the  ceramics.  A  change  of  the  transducer  temperature  or  the 
driving  voltage  will  also  produce  a  shift  of  the  resonant  frequency. 
Although  the  shift  is  not  necessarily  drastic,  it  is  sufficient  to 
throw  the  transducer  out  of  its  designed  operational  bandwidth.  This 
drastically  reduces  the  workipg  capability  of  the  transducer.  The  usual 
way  to  correct  this  problem  f  s  to  retune  the  transducer  to  its  original 
designed  frequency  by  changing  transducer  length,  or  to  employ  a  variable 
frequency  power  supply  to  change  the  operation  frequency.  Experience 
reveal o  that  thu  tv+-ui  tug  process  It  tfcliOMi  wd  time  ccnsyjir-g.  On 
the  other  hand,  variable  frequency  power  supplies  are  expensive.  This 
situation  presents  a  problem  in  attempting  to  use  the  less  expensive, 
fixed  frequency  power  supply  approach  to  sonic  processes.  In  order  to 
overcome  this  difficulty,  causes  of  the  shifting  must  be  understood; 
proper  measures  may  then  be  taken  to  reduce  the  shifting  or  to  make 
allowance  in  advance  for  the  shifting  and  to  offset  it.  In  this  chapter 
some  experimental  data  will  be  given  to  show  the  extent  of  the  shift 
with  the  prestress,  temperature  and  voltage.  Theoretical  analysis  of 
the  shifting  will  also  be  given. 

In  the  following  theoretical  formulations  for  various  resonators 
and  transducers  will  be  given  and  the  frequency  equations  will  be  solved 
numerically.  The  shifts  of  resonant  frequency  will  then  be  treated 
from  two  approaches:  (i)  from  the  effects  of  the  material  property 
Ganges,  and  (ii)  for  the  effect  of  prestress  from  determining  an  equiv¬ 
alent  spring  constant  for  the  ceramics  from  measuring  the  resonant 
frequency  of  the  simple  resonator  at  different  prestress  levels.  The 
theoretical  results  will  be  compared  with  the  experimental  values  ir 
the  final  section. 
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5-1  SONIC  RESONATORS 


consists^^several^iezoe^pof1^  the  f0rm  °f  a  sandwich  which 
cylinders  This  arrant  +•  i1"*0  dlscs  or  rin&s  clamped  between  metal 

resonato drlsUcSy °f  «* 
alone,  since  the  -ntlL  i  V  tlie  free  Piezoelectric  elements 

p  ezoelectrics  used  are  generally  thin. 


of  motio^Srea^h^^o^thf'a136115!?01'  °f  the  resonators  the  equations 
solutions  must  be  fmmrt  assembly  must  be  formulated.  The  general 

intofacL^t  VZLT  at  the  —  ®4  the 

CMrc  5  ssr 

“  included,  the  a-PUtudee“rv^bra?iona^n^foebeScSoSaet:d?erlalS 


thereto  ^  be  analyzed  asSumln8 

»«  (1)  the  sySetricS  1^irreao„SrrlSl-tb'he  reTators 

resonator,  and  fl)  the  p  -  +  so^°rJ  ( 2 )  th<_  asymmetrical  simple 

with  the  torn  rmived)  ThI (*,lc?  iS  “  0SU  ^  transducer 
models  will  be  used  -inth  thickness  plate  and  the  longitudinal  rod 

Will  be  compared  to  the  theorSS'onfs?'1115'’  eXperime'ltal  results 


(A)  The  Symmetrical  Simple  Resonator 


snnrhj-j  v.  A  symmetrical  simple  resonator  has  two  PZT  ceramic  rina<? 
sandwiched  between  two  steel  end  rods  connected  by  a  cento  boS  A 

■takTSK  voltoe^  *5  at  the  "“*«  °f  «»  resonator^to  serve 

this  resonator^"  FlSure  5'1  shows  the  schematic  diagram  of 


applied  to  the  SS  coders’  ™Ve  0,6017  ^  be 

SaiesStrlH«y  SPe8f ng>  ““  Mither  be  tiaasiSfafror^nmur 

than  the  rod  oJ’the'SSe  °tohthCOI,IPteftytnf  the  m°de:U  f°r  other 
one  or  the  other  r>n  the’  S  study»  they  be  approximated  as 

(less  tha^  'JV  Jer  hand’  the  copPer  strode  is  very  thin 

Viess  Phan  o.l  )  and  is  treated  as  rigid  and  massless  m  +h~  rJn  • 

the  resonator  will  be  analyzed  bv  model  W  th“  ?G  followine> 

rods.  y  Q  Dy  modeling  the  ceramics  as  plates  and 


the  eiS^hienfcircuiatofmthis:'Tting  adTtage  °f  the  s*™try  P^rty, 
techniques  [1].*  The  resultant ™  can  be  obtained  by  straightforward 

The  var\^s  sylholfhave  to^rSsf  **  Sh°™  “  FiS'  5"2' 


•Numbers  in  brackets  refer  to  references  at  the  end  of  this  chapter. 
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Z,10 

Z,g)  Z?  j 

9 - 1 

iBni 

vT& 

a  _  1 

Zi(i)  =  JaiZo^1)  tan  £  r^i  (5-1) 

Z2(i)  =  ~  jaizo^tysin  r^i  (5-2) 

>rtiere  i  =  1,  2,  3  for  elements  1,  2,  3>  respectively  (referring  to 
Fig.  5-1)  and 

Zq^1)  =  Zq(2)  =  (pgYg)7*  (metal  elements  1  and  2) 

Zq(3)  =  (pcc3§)'4  (ceramic  plates) 

ii  =  r2  =  w/cb,  Cb  =  (Ys/ps)'4 

r3  =  (n/(c^/pc)‘A  (5-3) 

The  other  quantities  are 

Z*g  =  &33&  3/2^3 

N  =  h33a3/p3|^3  =  h33/juZL|  (5-4) 

The  equivalent  circuit  can  be  simplified  to  that  of  Fig.  5-3  (a),  which 
in  turn,  is  equivalent  to  Fig.  5-3  (b). 


The  impedance  in  Fig.  5-3(a)  can  be  written  as 
Zt  =  Z^1)  +  Z^2)  +  Z2^2)  +  Za^3^  +  Z2(3) 

+  Zi^^Ze^VC2!^1'  +  z2^) 

=  daiZo^1)  tan  ri^l  "  ja2Zo^2Vtarl  r 2^2  -  ja3Z0^3Vtajl  r3^3 

Thus,  the  terminal  impedance  Z  is  as  follows: 

Z"1  =  i/V  =  (Z1^)"1  +  ("ZLe  +  ^t/N2)-1 
Z  =  ZL£(l  -  N2ZL.C/Zt) 

The  resonant  frequency  is  given  by  Z  =  0 ,  or , 


(5-5) 


(5-6) 


Q 

Zt  =  N2ZL£  =  -  jh332a3/u)J53p33 


(5-7) 


Substituting  the  definitions  of  Zt,  N  and  ZL|,  the  following  frequency 
equation  is  obtained; 

.  s 

a^A1)  tan  ri^i  “  a2Z0(2Vtan  ri^2  -  a3zQ(3)/tan  T 3^,3  =  -^33^3/^3^33 

(5-8) 

Since  z0(3)  and  z0(1)  are  comparable,  Ti^x  »  r2£2,  >>  r3^3>  ttms 

when  /cb  -»  n/2,  the  value  of  the  first  term  in  eq.  (5-8)  is  large 
and  increasing,  while  the  value  of  the  third  term  is  also  large  but 
decreasing.  This  implies  that  the  lowest  frequency  solution  to  eq. 
(5-8)  will  occur  near  o4x/cb  =  n/2.  Therefore,  for  the  lowest  resonant 
frequency,  the  approximation  tan  V2l2  -  r2£2,  tan  r3£3  =  r3f3,  can  be 
made  and  eq.  (5-8)  becomes 


tan  =  [  (a2  ^i/ai )  +  a3^i/ai  ^3)  (^33  -  e33h3l)/Ys]/riTi 


(5-9) 


One  can  show  that  c33  -  £3^133  =  c33>  Thus  eq.  (5-9)  i-s  simplified  to 

tan  T±l\  =  [(a2-2i/ai^2)  +  (a3-2i/ai^3)(c3l/Ys))/2i^i  (5-10) 

To  understand  the  physical  meaning  of  this  result,  eq.  (5-1°)  will  be 
compared  to  the  frequency  equation  obtained  by  replacing  the  ceramic 
rings  by  an  isotropic  elastic  rod  with  Young's  modulus  equal  to  Yc. 

It  is  easy  to  show  that  the  frequency  equation  for  this  system  is 

tan  Tx^x  =  [ (a2^i/ai^2)  +  (a3^i/ai^3)(YcAs)  V^i^i  (5-H) 
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Comparing  eq.  (5-10)  to  eq.  (5-11)  one  can  see  that  the  two  frequency 
equations  are  identical  if  one  substitutes  Yc  in  eq.  (5-10)  for  c3I. 

Thus  the  following  conclusion  is  made:  for  the  thickness  plate  model, 
the  resonant  frequency  of  a  simple  resonator  is  given  by  the  pure 
elastic  system  in  which  the  ceramic  has  an  equivalent  Young's  modulus 
equal  to  c33.  Stated  differently,  the  resonant  frequency  is  the  natural 
frequency  of  the  system  with  the  piezoelectric  rings  short  circuited. 


Equation  (5-10)  can  be  solved  numerically.  For  the  OSU  simple 
resonator,  the  following  numerical  data  were  used: 


Table  5  -1  -  Numerical  Data  for  the  OSU  Simple  Resonator 


For  steel  elements: 


al 

=  2.074  in.2 

a2  = 

0.1964  in.2 

£x 

=  5-0  in. 

£2  = 

0.265  in. 

Ys 

=  20.02  x  1010  N/m2 

Cb  = 

5050  m/sec 

=  29.0  x  10s  psi 

= 

19-88  x  104  in. /sec 

For  the  PZT-4  rings: 

a3 

=  1.325  in.2 

£3  = 

0.25  in. 

03I 

=  11.5  x  IO10  N/m2 

C33  = 

15.9  x  IO10  N/m2 

y! 

=  16.7  x  ICY3  psi 

= 

23.1  x  10s  psi 

=  6.45  x  IO10  N/m2 

Y?  = 

12.65  x  IO10  N/m2 

=  9*35  x  10s  psi 

= 

18.2  x  10s  psi 

After  substitution  of  the  above 
equation  becomes 

data  into  eq. 

■  (5-10),  the  frequency 

tan(w/x/cb) 

=  9.197/(1^ 

./c b)  (5-12) 

!!!iSOlUti°n  iS  ui::L//cb  =  thus  fPlate  =  0.226  CfcAx  =  8980  Hz. 

While  the  resonant  state  is  characterized  by  the  zeros  of  the  reactance, 
the  antiresonant  state  is  characterized  by  the  zeros  of  susceptance. 

For  this  problem,  the  antiresonant  frequency  is  given  by  Z 


-  '  - - 6xvui  uj  a  — >  00,  or 

ft  "  0*  ,f«ain  the  approximation  tan  r2i2  =  r2j>2,  tan  r3i3  =  r3£3, 

the  result  is  J 


tan  rxtx  -  [(a2^1/a1^2)  +  (a3^1/a1£3)(c3l/Ys)]/r1i’1  (5-13) 


120 


Physically,  from  this  result  it  is  equivalent  to  i: ay  that  the  ceramic 
has  an  equivalent  Young's  modulus  equal  to  c3l.  Hence,  the  antiresonant 
frequency  is  the  natural  frequency  of  the  system  with  the  ceramic  open 
circuited.  Substitution  of  the  numerical  data  yields 

tanCcoii/cb)  =  12.01/(0)^/%)  (5-14) 

Thus  fPlate  =  0.2309  %/h  =  9180  Hz 

The  longitudinal  rod  model — For  the  longitudinal  rod  model,  the 
only  modification  to  be  made  is  in  the  part  of  the  equivalent  circuit 
representing  the  ceramic.  The  equivalent  circuit  for  a  longitudinally 
vibrating  ceramic  for  the  rod  is  incorporated  and  the  new  equivalent 
circuit  of  this  model  is  shown  in  Fig.  5-4.  Some  of  the  symbols  appear¬ 
ing  in  the  figure  were  given  previously  in  eqs.  (5-l),  (5-2)  and  (5-3) • 
Additional  definitions  are 

z0(3)  =  (pcY3®),/»  for  the  ceramic  rod 
r3  =  «/<*<,  ,  cbc  =  (Y3D/pc)* 

ZLg  =  s33£3/jii)a3e33s33 

N  =  833/1^2  33Z^g  =  a3d33/i3s33  (5-15) 


Fig.  5-4  Equivalent  circuit  for  the  symmetrical  simple 
resonator-longitudinal  rod  model 


Again,  the  equivalent  circuit  is  simplified.  The  resulting  circuit  is 
exactly  the  same  as  that  of  Fig.  5-3  for  the  plate  model,  except  that 
the  symbols  have  different  definitions.  After  manipulation,  the 
impedance  is  found  to  be 


Z  =  ZL£  (l  .  N2ZLC/zt) 


(5-16) 


where 


Zt  =  ja1z0(1)  tan  -  ja2z0(2)/tan  rxZ2  -  ja3z0(3)/tan  r3f3 


(5-17) 


For  the  resonant  frequency,  Z  =  0  yields 


tan  =  a2z0(2)/a1z0(1)  tan  rxf2  +  a3zo(3)/a1z0(1)  tan  r3^3 

-  (aa/aiWf 3) (s3|  -  Ss^/ssIssaZoC1)  (5-18) 

With  the  same  argument  given  previously,  the  approximations, 

tan  Fxfiz  ±  rif2,  tan  r3i3  =  r3^3,  can  be  made.  Thus,  eq.  (c,-l8) 
becomes  ’  v'  7 

tan  w^/cb  =  [(a^/a^s)  +  (&3£i/ei1fi3)(Y3E/Yt,)]/{ul>1/ch)  (5-19) 

Equation  (5-I9)  can  be  regarded  as  representing  a  mechanical  system 
wi  h  a  short  circuit  across  the  faces  of  the  ceramic.  Substituting 
the  numerical  data  given  in  Table  5-1,  eq.  (5-19)  reduces  to 

tan(u)j>x/cb)  =  5-975/Mi/cb)  (5_20) 

This  has  a  solution,  (u^/cb)  =  1.3488,  or  f£od  =  8540  Hz.  The  anti¬ 
resonant  frequency  is  given  by  Zt  =  0,  which  yields 


tanoj^/cb  =  [(asix/a^)  +  (a^x/ax^aJCY^/YsJV^/cb) 


(5-21) 


Equation  (5-21)  can  be  regarded  as  representing  a  mechanical  system 
with  an  open  circuit  across  the  faces  of  the  ceramic.  Using  the  data 
of  Table  5-1,  eq.  (5-21)  becomes 

tan  uix/cb  =  9-937(ufa/cb)  (5-22) 

A  solution  to  this  equation  is  c^i/ck  =  1.4281,  or  frod  =  0  2273 
CbAi  =  9038  Hz.  a  •  (0 

Comparing  fPlate^  f£od,  fgl&te^  f^od  to  one  another,  one  can  see 

that  the  resonant  frequencies  of  the  above  two  models  differ  considerably, 
while  the  antiresonant  frequencies  do  not  differ  very  much. 
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The  spring  model- -The  third  model  being  formulated  is  the  "spring 
model".  For  this  model,  the  crystals  are  approximated  by  a  linear 
spring  with  its  spring  constant  remaining  to  be  determined.  A  simple 
resonance  test  is  required  to  assign  the  equivalent  spring  constant. 

It  is  hoped  that  this  constant  can  be  applied  to  other  resonators  of 
similar  configuration.  In  this  way,  the  complicating  factors  which 
result  in  discrepancies  between  experimental  resonant  frequencies  and 
the  one-dimensional  theoretical  results  can  be  disposed  of  simply. 
Furthermore,  the  complicated  piezoelectric  relations  can  be  eliminated. 
Figure  5-5 (a)  shows  this  model  schematically.  Taking  advantage  of  the 
symmetry  property  of  this  problem,  the  problem  is  simplified  as  shown 
in  Fig.  5-5(h).  The  equations  of  motion  are  similar  to  those  used 
earlier.  The  general  solutions  are 

ui(xi,t)  =  (Ai  cos  wxi/cb  +  Bj.  sin  cjXi/c^)  e*^ 

u2(x2,t)  =  (A2  cos  o^/c^  +  B2  sin  cox2/cb)  e*^  (5-23) 

The  boundary  conditions  are 

dui/dxi|  =  0 

lxx  =  0 
u2(0,t)  =  0 

=  -  u2(.02,t) 

aiYgdui/dxJ  =  a2Ysdu2/dx2  +  2K  u2(i2,t)  (5-24) 

I  xi  ~  x2  =  & 2 

Substituting  eq.  (5-23)  into  eq.  (5-24)  yields 

B]_  =  A2  =  0 

Ai  cos  w^x/cb  =  -  B2  sin  w^/ct,  -  A1a1Ys(u/cb)  sin  ^l/cb 
=  B2a2Y5(u)/ct,)  cos  “•fe/cb  +  2K  B2  sin  toi2/ Cb 

(5-25) 

The  frequency  equation  is  obtained  from  eq.  (5-25)  for  nontrivial 
solutions  of  kx  and  B2.  It  is 

tan  ti^/cb  =  (a2/ax  )/tan  u^2/cb  +  2K/(a1u>Ys/cb)  (5-26) 

Since  lx  »  iZi  tan  ui2/c^  ^  u£2/cb  so  that 

tan  to£i/cb  =  [a^/a^  +  2K/(Ysa1/jJ1 ) ]/(^x/cb)  (5-27) 

Substituting  the  numerical  data  of  Table  5-1  into  the  above  equation, 
the  solution  for  the  lowest  natural  frequency  can  be  plotted  against 


Ps 

x#- 

(b) 


Ps  »^S,a2*^S 


Fig.  5-5  The  spring  model  for  the  symmetrical,  simple 
resonator 


the  unknown  K  value  as  shown  in  Fig.  5-6.  Specifically,  if  one  uses 
the  strength  of  material  K  values,  i.e.,  if  2K  =  Y3Ea.,/X  =  5. 00  x  107 
lb/in.,  eq.  (5-27)  becomes 


tan  u}^/ Cb  =  g.Ol/Cw^i/cb) 


(5-28) 


The  solution  is  =  1.35*  Thus  the  natural  frequency  is 

fE  =  8545  Hz.  If  the  value  2K  =  Y3Dac/.0c  =  9-70  x  107  lb/in.  is  used, 
one  obtains 


tan  to.0! /cb  =  9.95/Mi/cb) 


(5-29) 


The  solution  is  w^/cb  -  1.4283,  thus  the  natural  frequency  is 
is  fD  =  9040  Hz.  The  above  two  natural  frequencies  are  indicated  in 
Fig.  5-6.  Theoretically,  comparing  eqs.  (5-19)  and  (5-21)  to  eq.  (5-27), 
the  frequencies  obtained  above  should  be  exactly  the  same  as  those  for 
the  longitudinal  rod  model.  However,  due  to  the  round-off  errors,  the 
values  are  slightly  different. 


(B)  The  Asymmetric  Simple  Resonator 


In  Fig.  5-7,  a  schematic  diagram  of  an  asymmetric  simple 
resonator  is  shown.  In  this  figure,  Part  1  is  considered  to  be  a  solid 
cylinder  which  consists  of  the  end  nut  and  a  portion  of  the  center  bolt 
screwed  together;  Part  2  is  a  hollow  cylindrical  steel  sleeve.  Part  3 
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Fig*  5-7  The  asymmetric  simple  resonator 


represents  the  remaining  center  bolt  while  Part  4  is  the  solid  steel 
cylinder  on  the  other  side  of  the  ceramic  rings,  and  Parts  5  and  6  are 
the  two  ceramic  rings . 


The  one-dimensional  wave  model--To  analyze  this  problem,  one¬ 
dimensional  wave  equations  will  again  be  employed  to  model  each  part 
of  the  resonator.  Boundary  conditions  at  the  ends  and  interfaces  must 
be  incorporated  to  solve  for  the  various  constants  as  well  as  the  fre¬ 
quency  equation.  The  equivalent  circuits  obtained  previously  can  be 
used  to  relate  the  forces  and  the  velocities  for  the  boundary  conditions , 
tor  this  problem,  this  will  yield  26  equations  for  26  unknown  terminal 
velocities,  forces  and  currents.  Numerical  calculation  of  the  resonant 
frequency  is  possible  by  looking  for  the  frequency  for  which  the  deter¬ 
minant  of  coefficients  of  the  above  unknowns  vanishes.  This  could  be 
determined  by  use  of  an  electronic  computer  but  would  consume  consid¬ 
erable  computational  time.  However,  by  using  the  following  scheme 
without  resorting  to  the  equivalent  circuit  method,  the  procedure  can 
be  simplified.  Thus,  from  this  point  on,  the  equivalent  circuit  will 
not  be  used  again.  The  general  solutions  presented  previously  will  be 
substituted  directly  into  the  boundary  conditions  and  the  frequency 
equation,  which  is  obtained  frcm  the  vanishing  determinant  of  the 
resulting  coefficient  matrix,  will  be  solved  numerically  on  the  computer. 


The  general  solutions  for  the  steel  members  1  through  4  are 
expressed  as  follows : 


ui(xi>t)  =  (%  cos  +  Bi  sin  fxj_)  e1- 


(5-30) 


where  T  -  w/c^,  c^  =  (Yg/ps)'4,  and  i  =  1,  2,  3>  4  for  members  1,  2,  3> 
4,  respectively.  The  velocities  and  forces  can  be  derived  as  follows: 
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Vi(xi,t)  =  jc^A*  cos  rx±  +  B*  sin  rx^  eJwt 

fi(xi,t)  =  Ysair(-Ai  sin  Txi  +  Bi  cos  Txi)  eJtl)t 
where  again,  i  =  1,  2,  3,  4.  For  the  ceramics 

Up(xp,t)  =  (Ap  cos  rcxp  +  Bp  sin  rcxp)  e-3wt 

VP^XP*^  =  J  (^p  cos  f  cxp  +  Ep  sin  TcXp)  e^'w^ 

^p(xp)>t)  =  ac[uzoc(-Ap  sin  TcXp  +  Bp  cos  TcXp)  -  £lp]  e*^ 
Va  =  -  h33[ (cos  rc£p  -  l)Ap  +  Bp  sin  rc£p]  +  ZL£  I® 
where  p  =  5  and  6  for  ceramic  rings  5  and  6,  respectively,  and 

rc  =  w/cbc 

UcaB/Pc)'4  for  thickness  plate  model 
Cfac  |  (Y3D/PC  )4  for  longitudinal  rod  model 

IQ 

0334c/jwac  for  thickness  plate  model 
Sal^c/j^ac^asSaB  for  longitudinal  rod  model 

^  _  |h33/juac  for  thickness  plate  model 

(g33/jwacs33  for  longitudinal  rod  model 


(5-31) 


(5-32) 


zoc  - 


I  (pcc33)y'  for  thickness  plate  model 


i (pcYaD  )v‘  for  longitudinal  rod  model 
The  boundary  conditions  are 


1) 

at  free  end  xx  =  0 

fi(0,t)  =  0 

(5-33) 

2) 

at  interface  xx  =  £lf 

x2  =  0  and  x3  =  0 

«i(ii»t)  =  u2(0,t) 

(5-3*+) 

%Ui,t)  =  u3(0,t) 

fi(fi,t)  =  f2(0,t)  + 

f3(0,t) 

(5-35) 

3) 

at  interface  x2  -  £2i 

o 

ii 

in 

X 

^2(^2,t)  =  u5(0,t) 

f2(i2,t)  =  f5(0,t) 

(5-36) 

*0 

at  interface  x5  =  £c, 

» 

II 

<** 

o 

u5Uc,t)  =  -u6Uc,t) 

f5Uc,t)  =  f6Uc>t) 

(5-37) 
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(5-38) 


5)  the  current  versus  voltage  relations 

V  =  /  Es  .  dx5 

V  =  /  E6  •  ox6 

6)  at  end  x4  =  0 

f4(0,t)  =  0  (5-39) 

7)  at  interface  x4  =  U»  x3  =  £3,  and  x6  =  0 

^3(i3,t)  =  -  U4(^4,t) 

u6(0,t)  =  -  u4U4,t) 

f3(^3,t)  +  f6(o,t)  =  f4(*4,t)  (5-40) 

In  susaary,  there  are  14  equations  for  the  14  unknown  A's,  B's,  and  I's 
resulting  from  substituting  the  general  solutions  into  the  above 
boundary  conditions .  This  is  quite  a  simplification  compared  to  the 
equivalent  circuit  approach  which  has  26  equations  for  26  unknowns. 

The  reason  for  this  simplification  can  easily  be  explained:  from 
minimum  energy  principle,  it  can  be  proved  that  one  cannot  specify 
force  and  velocity  at  the  same  end  simultaneously.  Thus,  al  1  the  end 
forces  and  velocities  appeared  in  the  equivalent  circuit  approach 
cannot  be  mutually  independent .  Therefore  there  must  be  one  equation 
relating  the  end  force  and  velocity  at  each  end.  For  our  problem, 
there  are  6  members  and  thus  12  ends .  This  is  why  there  are  12  ecr  Lations 
and  12  unknowns.  In  the  following,  the  general  solutions  will  be  . sub¬ 
stituted  into  the  boundary  conditions  and  then  the  relations  will  be 
simplified  to  give  a  frequency  equation  of  a  7  x  7  determinant. 

The  substitutions  yield, 

Bj.  =  0 

Ax  cos  r ix  =  a2 

Ai  cos  r&x  =  A 3 

-ajAx  sin  r £x  ~  s.2B2  t  a3B3 

A2  cos  F£2  +  B2  sin  T£2  =  A5 

a2uz0(-A2  sin  V t2  +  B2  cos  Y  Z3)  =  a5(uz0cB5  -  t;i§) 

Ascos  rc£c  +  B5  sin  rc£c  =  -(Ag  cos  Yc£c  +  Bs  sin  Tc£c) 

(“A5  sin  rc/,c  +  Bs  cos  T q£c  +  As  sin  Tc£c  -  Bg  cos  rc£c) 

-  (£/uzoc)(l§  -  Ig)  =  0 

V*  =  -  h33[(cos  rcec  -  1 ) Ag  +  Bg  sin  rc£c]  +  ZL{jj  1% 

V*  =  -  h33[(cos  rc£c  -  1)AS  +  Bs  sin  Tc£c]  +  ZL£  1%  (5-41) 
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A3  cos  YZ3  +  B3  sin  YZ3  =  -  A4  cos  r Za 
Ag  -  A4  cos  YZ* 

a3wz0(-A3  sin  r^3  +  B3  cos  r Z3)  +  ac(ujzocB6  -  (; Ig) 
=  a4u>z0(-A4  sin  r^4) 

The  simplified  relations  are 

Bi  =  0,  B4  =  0,  A2  =  A3  =  Ax  cos  rix 

A-*  =  -  (Ai  cos  cos  r^3  +  B3  sin  r^3)/cos  TZa 

As  =  Aj_  cos  r^i  cos  T Zz  +  B2  sin  TT2 

A6  =  -  (Ax  cos  TTx  cos  TT3  +  B3  sin  TT3) 


(5-43 


,  DiJ 

(a  7x7  matrix) 


(5-42 


where  =  0  and  =  0  except  for 

=  sin  VZxj  Bx2  =  a2/ax>  Dx3  =  a3/ax 

D21  =  (a2/ax )  cos  r Z\  <iin  YZ2,  D22  =  -  (a2/ax )  cos  VZ2 
B24  =  (ac/al  )  (zOc/zo)  >  D25  =  -  (ac/al)(C/u,zo) 

=  cos  r Z\  cos  Y  Z2  cos  rcZc  -  COS  rz±  cos  TT3  cos  rcTc 
Z>32  -  sin  rz2  COS  rcTc,  D33  =  -  sin  rz3  cos  TCTC 
B34  =  sin  rcTc,  D35  =  sin  r^Zc 

D4i  =  cos  r^L  cos  r^2  sin  FCTC  +  cos  TZx  cos  VZ3  sin  rcTc 
D42  =  sin  VZ2  sin  r cZCt  D43  =  sin  V Z3  sin  Y  CZC 
D44  =  -  cos  r c^C>  B4g  =  COS  r c^c»  D4g  =  -  £/u>zoc 
D47  =  ^/uJZqc 

D51  =  -  (a3/ax)  cos  r^x  sin  Y  Z3  -  (a4/ax)  tan  YZ4  cos  T Z\  cos  Y  Z 
D53  =  (a3/ai)  cos  T£3  -  (a4/ax)  tan  Y Z4  sin  YZ3 
D55  =  (ac/ax ) (z0c/zo)>  B57  =  (ac/ax ) (£/uzq) 


Dei  =  (cos  rc*c 
De2  =  (cos  Vcic 
D64  =  sin  Vcic, 


l)  COS  r ii  cos  Viz 
l)  sin  VI p 

D6s  =  "  Z^S/h33 


(l  -  cos  Tc£c)  cos  ^£1  cos  r £3 
(l  -  cos  rc£c)  sin  r^3 
sin  D77  =  ~  ZL£/h33 


=  -  Va/h: 


The  frequency  equation  is  given  by 


det[Pij]  =  0 


(5-44) 


The  spring  model--As  for  the  symmetric  simple  resonator,  the  spring 
model  for  the  asymmetric  simple  resonator  can  be  obtained  by  simply 
modeling  the  ceramics  as  an  elastic  spring  (Fig.  5-8).  The  equations 
of  motion  and  their  corresponding  general  solutions  are  the  same  as 
previously  given  in  (5-30)  to  (5-31)  for  the  steel  members.  Boundary 
"conditions  1,  2  and  6  are  the  same,  while  at  steel-ceramic  interfaces 
the  original  set  of  equations  are  replaced  by  the  new  set  of  equations 
involving  the  spring  constant.  These  new  equations  are: 


1)  at  x2  =  Izi  ,  %  ,  .  . 

f2U2,t)  +  K[u2Ui,t)  +  luU^t)]  =  0 

2)  at  x4  =  i 4,  and  x3  =  i3, 

U4(£4,t)  =  -  u3(i3,t)  .  _ 

f4(i4,t)  =  f3U3,t)  -  K[u2U2,t)  +  U4(i4,t)J 


(5-45) 


(5-46) 


Copying  down  the  rest  of  the  valid  boundary  equations  and  substituting 
the  general  solutions  into  the  above  three  boundary  conditions,  one 
obtains  the  following  relations  after  some  manipulation: 


Bi  =  0,  B4  =  0 
A2  =  A3  =  Ai  cos 

A4  =  (-Ax  cos  r^i  cos  r.03  -  B3  sin  Vi3)/cos  r,e4 


Dll  Di2  Di3 
D2i  D22  D23 
D3i  D32  D33 


(5-47) 


(5-48) 


where 


Du  =  sin  d12  =  a2/ai»  Di3  =  aa/ai 

D2i  =  (a2/ai)  cos  Vtx  sin  r.02  +  cos  f£i  (cos  V i3  -  cos  Viz) 
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Fig.  5-8  The  spring  model 

for 

the  asymmetric 

simple  resonator 

D22 

=  -(aa/aj.)  cos  r  I2  ■ 

■  Cl  sin 

r  is 

1>23 

=  £1  sin  T£3 

i>3l 

=  -  (a3/ aj. )  cos  Tii 

sin  Vl3 

-  (« 

i2/ai)  cos  rij.  sin  F£2 

-  (a4/a1)  cos  r^i 

cos  ri3 

sin 

r £4/cos  T£4 

i>32 

=  (a2/a1)  cos  ri2 

D33 

=  (a3/ai)  cos  r-03  - 

(a4/a1) 

sin 

n3  sin  r^4/cos  ri4 

In  the  above,  the  symbol  =  K  Cb/Ys&iW  has  been  used.  The  frequency 
equation  is  given  by 

det[Dij]  =  0  (5-49) 

Again,  this  can  be  solved  numerically.  However,  numerical  solutions 
were  not  actually  obtained  in  this  study. 

(C)  The  P-7  Resonator 

The  P-7  resonator  is  essentially  the  OSU  P-7  transducer  with¬ 
out  the  stepped  horn.  It  differs  from  the  asymmetric  simple  resonator 
only  in  its  dimensions  and  its  cross  sections.  It  does  not  have  a 
uniform  cross-sectional  area  along  its  length.  Also,  it  has  a  hexagonal 
cylindrical  nut  and  steel  sleeve.  Figure  5-9  shows  the  schematic 
diagram  of  the  P-7  resonator.  The  set  of  equations  are  identical  to 
those  of  Section  (b)  except  that  the  geometrical  dimensions  are  different. 
The  following  data  are  used  for  the  P-7  resonator. 


Fig.  5-9  The  P-7  resonator 


Table  5-2  -  Dimensions  for  the  P-7  Resonator 


-  1.25  in. 

^2  =  2.75  in. 

^3  =  3-35  in. 

**  =  O.97  in. 

*  £e  =  £c  =  0.25  in. 


ai  =  1.947  in.2 
as  =  1.505  in.2 
a3  =  0.442  in.2 
a4  -  2.405  in.2 
a5  =  a6  =  ac  =  1.325  in. s 


equation,  eq.  (5^49)hca^abersSvPd°Perty  COefficients  the  frequency 
Puter.  The  lowest  resonant  m  ^ranl  com- 

tudinal  rod  model  cf  tTcer^cUenCy  iS  fr  =  15’33  Wz  f°r  the  l°ngi. 

material  x'lo^lh/^’8  the  sJatl0  strength  of 

in  eq.  (5-49  ,  the’lLst  fr{™  “<1  *  =  '‘•«5  51  10'  lb/in. 

4  =  16.4  kHz,  fg  „  17  6  J.1  ™  we  found  to  be 

the  one  for  the  longitudinal  rod' ^dS”g1ff,ab°Ve  tU°  ^^^oies  wit 
model  gives  a  frequency  which  is  too  high.  aPl,aren+-  that  the  spring 


5-2  THE  P-7  TRANSDUCER 


stepped  horn.  FigSe^-l^isV^8^^  aS  E  P"7  resonator  plus  a 
Although,  as  a  fiSJ  appioxSation  tT  °f  the  P'7  transduce 

*ined  by  using  the  «**  ^  deter 

b  pproacn,  the  more  exact  solution 


Fig.  5-10  The  P-7  transducer 


should  be  calculated  from  the  more  exact  model.  Especially,  when  the 
transfer  function  between  the  electrical  terminals  and  the  tip  of  the 
transmission  line  is  to  be  decided,  the  piezoelectric  equations  must 
be  taken  into  consideration.  In  this  section,  the  method  used  in 
treating  the  P-7  resonator  (Section  5-1  (C))  will  be  extended  to 
include  the  stepped  horn. 

The  general  solutions  for  the  steel  members  1,  2,  3,  4,  and  7  are 
given  by  eqs.  (5-30)  and  (5-31),  while  for  the  ceramic  rings  5,  6, 
e<l*  (5-30)  applies,  in  addition  to  the  boundary  conditions  given  by 
eqs.  (5-33)  to  (5-38),  the  following  equations  arise  from  the  attach¬ 
ment  of  /the  stepped  horn  7i  They  are: 

1)  at  the  free  end  x7  =  0,  f7(0,t)  =  0,  so  that 

b7  =  0 

2)  at  interface  x7  =  ZT,  and  x*  =  0,  u7U7,t)  =  U4(0,t)  and 
f7(i!7,t)  =  f4(0,t);  thus 

A7  cos  m7  =  A4 

-a7A7  sin  V  Z7  =  a4B4 

3)  at  interface  x4  =  £4,  x3  =  Z3,  and  Xq  =  0,  114(^4, t)  =  Ue(0,t)  = 

-  u3U3,t),  and  f4(£4,t)  =  f6(0,t)  +  f3U3,t);  thus, 

A4  cos  r^4  +  B4  sin  T£4  =  A6 

a4wz0(-  A4  sin  T i4  +  B4  cos  T£4)  =  ac(u)zocB6  -  (;i§) 

+  a3'jz0(-  A3  sin  V i3  +  B3  cos  T i3) 

A4  cos  r^4  +  B4  sin  r £4  =  -  (A3  cos  V  Z3  +  B3  sin  TZ3) 

Thus,  there  are  16  equations  for  16  unknowns.  Again,  the  above  equations 
can  be  simplified  and  the  frequency  determinant  obtained  and  solved 
numerically  on  an  electronic  computer.  Substituting  the  numerical  data 
previously  given,  the  lowest  resonant  frequency  versus  a  5/8  in.  diameter 
horn  with  a  length  of  approximately  5  in.  is  obtained  and  shown  in 
Fig.  5-U  as  the  "1-D  Rod  Model." 
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resonant  frequency  versus  the  length  of  horn  of  a  P-7  transducer 


.the  spring  model --Figure  5-12  shows  the  spring  model  for  the  P-7 
transducer.  It  is  obtained  from  Fig.  5-10  by  modeling  the  ceramic 
rings  as  an  elastic  spring  with  spring  constant  K.  The  formulation  of 
the  equations  of  motion  are  similar  to  the  previous  cases  and  will  not 
be  repeated  here.  After  substituting  the  general  solutions  into  the 
boundary  conditions  the  frequency  equation  is  found  to  be 


det[Dij]  =  0 

where  Dj_j  is  a  3  x  3  matrix  with 

Du  =  sin  0^2  =  a2/aij  ^13  =  a3/a1 

D2i  =  (as/aj,)  cos  Y  Zx  sin  Y  Z2  +  cos  Y  Zx  cos  r_g3 

-  £i  cos  YZX  cos  YZ 2 
D22  =  -  {^2/^1)  cos  r^2  -  £1  sin  yz2 

°23  =  £1  Sin  r^3 

D31  =  (a4/a!)  [sin  ri4  +  (a7/a4)  tan  Y  l7  cos  Y  ZA]  cos  Y  Zx 
cos  YZ3/  [cos  YZ4  -  (a7/a4)  tan  YZ7  sin  r£4] 

+  (a3/ai)  cos  YZX  sin  YZ3  -  (;x  cos  Vlx  (cos  r.g3  -  cos  Y  l2) 

D32  ~  £1  sin  r  £2 

D33  =  (a4/a1)  [sin  Y  ZA  +  (a7/a4)  tan  Y  l7  cos  r^4]sin  Y l3/ 

[cos  r^4  -  (a7/a4)  tan  Y  Z7  sin  r^4]  -  (a.3/ix)  cos  Y Z3 
-  £].  sin  YZ3 

The  resonant  frequency  (lowest)  for  K  =  KE  =  2.5  x  107  lb/in.  was  found 
numerically  for  several  lengths  in  the  vicinity  cf  5  in.  as  shown  in 
Fig.  5-11  as  the  "Spring  Model."  For  a  more  direct  comparison,  the 
results  of  the  longitudinal  rod  model  and  the  fixed-free  horn  are  also 
shown  in  the  same  figure .  It  is  seen  that  the  slopes  of  the  one¬ 
dimensional  model  and  the  spring  model  are  the  same,  while  that  of  the 
fixed-fPee  horn  is  somewhat  steeper. 
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5-12  The  spring  model  for  the  P-7  transducer 


In  conclusion,  it  may  be  worthwhile  to  differentiate  the  physical 
meaning  <~f  the  resonant  frequency  from  that  cf  the  antlrescr.ani  frequency . 
At  the  resonant  frequency  the  terminal  current  is  infinite  for  a  finite 
driving  voltage,  since  the  terminal  impedance  is  zero.  Consequently, 

Hut  jvirl  inlf  vnljcH.y  is  i"  fil  ♦jJWrq'ifc  iit  tip  Wi'ljri'y  ,  ,\'nl  point** 
Contrarily,  the  particle  velocity  will  be  zero  at  the  antiresonant 
frequency,  since  the  terminal  impedance  is  infinite  and  the  terminal 
current  is  thus  zero  for  a  finite  driving  voltage.  However,  it  must  be 
noted  that  the  pure  mechanical  systems  (neglecting  the  piezoelectric 
effects  completely)  do  have  a  natural  frequency  at  the  antiresonant 
frequency.  This  is  evident  from  the  vanishings  of  the  pure  mechanical 
impedances  for  the  equivalent  circuit  systems.  Therefore,  in  order  to 
drive  the  system  into  its  resonant  state  at  the  antiresonant  frequency, 
the  mechanical  forcing  function  ( force  or  displacement  }  must  It-  pfwidfed 
at  one  of  the  mechanical  terminals.  An  electrical  driving  voltage  will 
not  be  able  to  drive  the  system  into  its  resonant  state  at  the  anti¬ 
resonant  frequency.  Therefore,  in  order  to  obtain  maximum  energy 
transauc' ion,  the  system  must  be  driven  at  its  resonant  frequency. 


5-3  STATIC  BIAS  AND  TEMPERATURE  EFFECTS  ON 

TRANSDUCER  RESONANT  FREQUENCIES 

As  mentioned  earlier,  ferroelectric s  are  highly  nonlinear  in 
material  properties.  Their  properties  are  also  subject  to  change  with 
passage  of  time,  a  phenomenon  called  the  aging  effect.  Experiments 
have  shown  that  a  major  disturbance  such  as  large  temperature  change, 
or  high  stress  level,  or  intense  electric  field,  will  not  only  change 
the  properties  themselves,  but  also  trigger  a  new  cycle  of  aging.  Thus 
a  complete  understanding  and  precise  prediction  of  the  behavior  of  a 
ferroelectric  material  is  quite  difficult.  However,  thanks  to  various 
researchers ,  sufficient  useful  information  has  been  accumulated  on 
PZT-4  to  make  possible  a  better  understanding  of  its  properties. 

In  the  OSU  Sonic  Power  Laboratory,  the  resonant  frequency  of  a 
transducer  has  been  observed  to  shift  with  different  prestress  bias 
levels.  It  also  shifts  with  the  change  of  temperature  of  the  trans¬ 
ducer.  Moreover,  a  higher  driving  voltage  will  always  decrease  the 
resonant  frequency.  These  effects  greatly  hamper  efforts  to  operate 
the  transducer  at  a  fixed  frequency.  Because  the  shift  of  resonant 
frequency  of  a  transducer  shifts  it  away  from  the  fixed  supply  fre¬ 
quency,  the  transducer  is  detuned;  i.e.,  it  has  an  amplitude  of  vibra¬ 
tion  much  less  than  the  designed  amplitude  at  the  power  supply  fre¬ 
quency.  Figure  5-13  illustrates  the  detuning  effect  with  the  shift  of 
resonant  frequency.  The  detuned  transducer  is  seen  to  have  a  very 
small  response  at  the  power  supply  frequency,  although  the  shift  of 
frequency  is  small  compared  to  the  original  resonant  frequency.  When 
a  transducer  is  detuned,  one  way  to  correct  this  problem  is  to  tune  it 
back  to  its  original  resonant  frequency  by  changing  the  length  of  the 
horn,  since  the  resonant  frequency  of  a  transducer  is  very  sensitive  to 


Fig.  5-13  The  detuning  effect  on  a  transducer 


the  horn  length  (see  Fig.  5-H).  However,  this  process  of  correction 
is  rather  tedious. 

The  above  problem  prompts  this  investigation.  In  the  following, 
the  shift  of  resonant  frequency  of  a  transducer  due  to  the  prestress 
bias  level,  the  temperature  change  and  the  voltage  supply  amplitude 
will  be  investigated.  The  aging  effect  will  be  compared  to  the  experi¬ 
mental  ones ,  where  they  aire  available . 

Another  aspect  of  the  work  will  be  an  experimental-theoretical 
approach  utilizing  the  spring  models  previously  proposed.  This  approach 
is  used  because  the  ceramic  rings  are  hard  to  model  more  precisely  since 
they  are  neither  a  plate  nor  a  rod.  Although  the  more  exact  two- 
dimensional  model  can  be  formulated,  the  solution  is  difficult  to  obtain 
and  will  not  be  used  in  this  study.  It  is  envisioned  that  the  spring 
model  which  has  a  spring  constant  K,  that  may  be  found  from  experiments 
with  the  simple  resonator,  will  enable  one  to  determine  the  shift  of 
resonant  frequencies  in  the  other  resonators  and  transducers  using  the 
same  ceramics. 

(A)  Summary  of  Prestress,  Temperature  and 
Voltage  on  PZT-4  Properties 

It  has  been  suggested  that  the  shifts  of  the  material  proper¬ 
ties  of  piezoelectrics  are  dictated  by  the  intrinsic  nonlinearities  of 
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PRESTRESS  (ksi) 


Fig.  5-16 


S33  and  S 


E 

33  versus  prestress  as  calculated  from 


Fig. 


5-14 


k33  =  ^33/ ( S33£33 ) 


the  elastic  compliance  s3§  or  the  Yorung's  modulus  Y3E  =  l/s3|  can  be 
deduced  from  Pig.  5-14.  The  result  is  shown  in  Fig?  5-15.  Elastic 


compliance  s33  can  be  obtained  by  the  relation 


S33  =  s3l 


(d33)?/033 


(5-50) 


Figure  5-16  shows  the  compliance  s3l  and  s3|  versus  the  parallel  com- 

annliertStre?S  ’  Later>  data  From  these  figures  will  be 

applied  to  calculate  the  effect  of  prestress  on  the  resonant  frequency. 


(B)  Effects  of  Pre stress,  Temperature  and 
Voltage  on  Resonant  Frequency 


aPPlyin§  the  changes  of  material  properties  given  in  the 
ixteratx.re,^  summarized  in  Section  (A)  ,  to  the  frequency  equations 
for  the  various  resonators  and  transducers,  the  new  resonant  frequencie 

tl  1fed'  In  th±S  Way  the  effects  of  the  Prestress,  tampera- 

ture,  and  voltage  may  be  evaluated.  In  the  following,  the  prestress 

e  lect  will  be  treated  and  then  the  temperature  and  voltage  effects 
will  be  discussed. 


,  lj?e  Prestress  effect  on  resonant  frequency— In  Section  (A),  Figs. 
7~XH’  •3~-L5,  and  5-I0  have  shown  the  material  coefficient  changes  at 


different  pre stress  levels  that  are  needed  for  calculating  the  resonant 
frequencies.  For  example,  from  Fig.  5-15  the  Young's  modulus,  Y3^,  at 
a  certain  prestress  level  can  be  obtained.  Substituting  this  value 
together  with  other  data  which  are  not  changed  by  the  prestress  into 
the  frequency  equation  (5-19),  one  can  calculate  the  resonant  frequency 
at  this  prestress  level  for  the  symmetric  simple  resonator.  The  results 
of  this  calculation  are  shown  in  Fig.  5-17-  Similarly,  for  the  P-7 
resonator  and  the  P-7  transducer,  one  can  obtain  the  necessary  material 
coefficient  variations  from  Figs.  5-14,  5-15  >  a*id  5-l6  at  various  pre¬ 
stress  levels  and  then  substitute  them  into  the  frequency  eq.  (5-43) 
for  the  P-7  transducer  to  get  Figs.  5-18  and  5-19-  From  the  above  three 
figures,  one  can  see  that  the  resonant  frequencies  increase  with  an 
increase  in  prestress  level. 

The  temperature  effect  on  resonant  frequency- -Not  enough  data  is 
available  in  the  literature  for  a  direct  evaluation  of  this  effect. 

The  only  material  coefficient  variation  given  [l4],  [15]  pertaining  to 
our  problem  is  €33  which  shows  little  increase  from  -60°F  to  200  F.  t 
Thus ,  more  experiments  must  be  conducted  on  the  variations  of  the  material 
properties  with  temperature  changes.  However,  it  is  possible  to  show  by 
a  fairly  simple  calculation  that  an  increase  in  temperature  will  tend  to 
release  a  portion  of  the  prestress  in  the  ceramics  and  to  decrease  the 
resonant  frequency. 

The  prestressed  ceramic  and  center  bolt  of  a  simple  resonator  is 
shown  in  Fig.  20.  Suppose  at  room  temperature  the  ceramic  has  a  total 
length  g, 2.  The  original  prestress  at  room  temperature  was  f0 .  Then, 
at  temperature  Ad  above  room  temperature,  the  decrease  in  the  prestress 
f 0  -  f  can  be  shown  to  be 

f o  -  f  =  (cM2  -  ac^3)A9/(l/Kc  +  l/K 2)  (5-51) 

where  a2  and  etc  are  the  linear  thermal  expansion  coefficient  for  steel 
and  the  PZT-4  ceramic,  respectively,  and  the  K’s  represent  the  static 
spring  constants  for  the  ceramic  bar  and  the  steel  center  bolt. 

For  the  simple  resonator  given  in  Section  5-l(A),  Kc  =  2.5  x  107 
lb/in.  and  K2  =  1.12  x  107  lb/in.  For  steel  a?  =  6.5  x  10-S/°F,  and 
for  the  ceramic  OLq  =  2.5  x  10-S/°C  =  1.4  x  10  D/°F.  Utilizing  the 
values  in  eq.  (5-51) ,  one  obtains 

f 0  -  f  =  20.3  A0  (5-52) 

This  relation  is  plotted  as  shown  in  Fig.  5-21.  Using  eq.  (5-52)  a 
theoretical- experimental  approach  will  be  introduced  to  evaluate  the 
resonant  frequency  shift  due  to  the  prestress  releasing  effect.  For 
a  P-7  transducer,  the  typical  prestress  on  the  ceramic  is  f0/ac  =5-4 
ksi  (or,  125  ft -lb  torque  on  the  3/4-in.  diameter  stress  bolt).  Thus, 
from  eq.  (5-52),  the  prestress  f/ac  at  any  temperature  can  be  calculated 
for  the  symmetric  simple  resonator.  After  resonant  frequency  versus 
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Fig.  5-17  The  prestress  effect  on  the  resonant  frequency  of 
a  simple  resonator 


Fig.  5-18  The  pres tress  effect  on  the  resonant  frequency  of 
a  P-7  resonator 
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Fig.  5-19  The  prestress  effect  on  the  resonant  frequency  of  a 
P-7  transducer 
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Fig.  5-20  The  prestressed  ceramic  and  center  bolt 
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Fig.  5-21  Prestress  release  versus  temperature  change  on  a  simple 
resonator 


prestress  relations  are  determined  experimentally,  the  resonant  fre¬ 
quency  versus  temperature  relations  can  be  deduced.  Particularly,  for 
the  symmetric  simple  resonator,  the  experimental  resonant  frequency 
versus  prestress  relation  is  given  later  in  Fig.  5-34  (Section  (d) ) ; 
the  result  of  this  approach  is  shown  in  Fig.  5-22. 

For  the  asymmetric  simple  resonator,  the  P-7  resonator  and  the 
P-7  transducer,  the  prestressed  portion  is  shown  in  Fig.  5-23.  This 
time  the  same  type  of  derivation  leads  to 

f o  -  f  =  &{a3£ 3  -  <fej2  -  aci5)/(l/K2  +  l/K3  +  l/Kc)  (5-53) 

Substituting  the  data  given  in  Table  5-2  and  the  material  properties, 
one  obtains  K2  =  I.587  x  107  lb/in.,  K3  =  O.3826  x  107  lb/in.,  and 
eq.  (5-53)  becomes 


f o  -  f  =  8.78  A0  (5-54) 

Again,  experimental  resonant  frequency  versus  prestress  relations  were 
obtained,  as  shown  later  in  Figs.  5-35  and.  5-36,  for  the  P-7  resonator 
and  the  P-7  transducer,  respectively.  The  stress  release  effects  are 
obtained  and  shown  in  Figs .  5-24  and  5-25  for  the  P-7  resonator  and 
the  P-7  transducer,  respectively,  by  employing  eq.  (5-54)  and  Figs.  5-35 
and  5-36. 
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Fig.  5-22  Theoretical -experimental  temperature  effects  on  the 
resonant  frequency  of  a  simple  resonator 
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Fig.  5-23  The  prestressed  ceramic,  center  bolt 
and  metal  sleeve 
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Fig.  24  Theoretical-experimental  temperature-prestress  release 
effects  on  a  P-7  resonator 


Frcsn  the  above  results,  one  can  see  that  the  prestress  release  is 
a  linear  function  of  the  temperature  change.  As  a  result,  the  tempera¬ 
ture  increase  always  reduces  the  resonant  frequency  of  a  resonator  or 
a  transducer.  In  a  later  section  one  will  see  that  this  is  in  accord 
with  the  experimental  observation. 

The  electric  field  effect  on  resonant  frequency--As  mentioned  in 
Section  (A),  there  is  not  sufficient  available  literature  to  evaluate 
this  effect.  The  only  plausible  reasoning  which  can  be  put  forth  here 
is  that  a  higher  voltage  tends  to  cause  a  higher  amplitude  of  vibration 
in  the  transducer.  A  bigger  dynamic  stress  will  be  realized  in  the 
ceramics.  This  according  to  the  literature  [15],  tends  to  reduce  the 
Young's  modulus  which  in  turn  causes  a  resonant  frequency  drop. 
Berlincourt  [15  ]  gave  data  which  showed  that  a  peak  dynamic  stress  of 
4  ksi  corresponded  to  a  5 %  decrease  in  Young's  modulus  Yi^.  However, 
data  on  the  more  useful  Y3^  is  not  available. 

(C)  Effect  of  Prestress  on  Resonant  Frequency, 

Based  on  the  Spring  Model 

The  ceramic  rings  used  in  our  resonators  and  transducers  can 
be  modeled  strictly  neither  as  a  plate  nor  a  rod.  It  was  suggested  in 
Section  5-1  to  model  the  ceramic  rings  as  a  linear  spring  leaving  the 
spring  constant  K  to  be  determined.  In  so  doing,  it  is  envisioned  that 
an  experimentally  determined  K  value  for  the  symmetric  simple  resonator 
at  a  certain  prestress  level  can  also  be  applied  to  the  other  resonators 
and  transducers.  It  is  envisioned  that  the  spring  models  so  applied 
will  account  for  the  resonant  frequency  shifts  due  to  the  prestress 
effect.  To  find  the  spring  constant,  K,  at  a  given  prestress  level, 
the  resonant  frequency  is  obtained  from  the  experimental  resonant  fre¬ 
quency  versus  prestress  relation  for  a  symmetric  simple  resonator  given 
by  the  later  Fig.  5-34.  Then,  the  spring  constant  is  determined  from 
Fig.  5-6  at  this  frequency  corresponding  to  this  prestress  level. 

Figure  5-26  gives  the  K  versus  the  prestress  relations  so  obtained. 

The  data  points  can  be  approximated  by  a  straight  line  which  does  not 
pass  through  the  coordinate  origin. 

Having  acquired  the  K  values  at  different  prestress  levels, 
those  values  may  be  applied  to  calculate  the  resonant  frequencies  of  a 
P-7  resonator  and  transducer  at  different  prestress  levels  utilizing  the 
spring  model  formulations  presented  in  Section  5-1.  For  convenience,  the 
numerically  obtained  resonant  frequency  of  the  P-7  resonator  and  the  P-7 
transducer  versus  the  spring  constant  K  is  plotted  in  Figs.  5-27,  5-28 
and  5-29,  respectively.  Now  at  any  prestress  level,  the  spring  constant 
K  can  be  obtained  from  Fig.  5-26.  The  resonant  frequencies  can  then  be 
obtained  from  Figs.  5-27,  5-28  and  5-29  for  the  P-7  resonator  and  the 
P-7  transducer,  respectively.  The  resulting  resonant  frequency  versus 
prestress  relations  are  given  in  Figs.  5-30  and  5-31.  It  is  seen  that 
the  resonant  frequency  increases  with  an  increase  in  prestress,  also 
that  the  rate  of  increase  of  the  resonant  frenuency  with  the  prestress 
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RESONANT  FREQUENCY  (kHz) 


Resonant  frequencies  of  P-7  transducers 
of  different  lengths  versus  K 


ig.  5-30  The  theoretical  prestress  effect  on  resonant  frequency  of  the  P-7  resonator 
the  spring  model 


is  quite  large  at  the  beginning  of  the  prestress  and  then  appears  to 
level  off  asymptotically  at  higher  prestress  to  a  certain  value  for 
each  resonator  or  transducer. 

(D)  Resonant  Frequency  Experiments 

Seme  simple  experiments  have  been  conducted  in  the  OSU  Sonic 
Power  Laboratory  to  measure  the  resonant  frequencies  of  various  resona¬ 
tors  and  transducers  for  different  envii  ^nments .  The  experimental  set¬ 
up  is  shown  schematically  in  Fig.  5-32. 

The  prestress  effect  on  resonant  frequency- -This  experiment  has 
been  conducted  on  the  symmetric  simple  resonator,  the  P-7  resonator  and 
the  P-7  transducer.  The  prestress  was  applied  by  a  torque  machine  on 
which  accurate  torques  on  the  end  nut  were  measured.  The  frequency  at 
every  torque  reading  was  then  measured  on  a  frequency  counter .  Although 
an  empirical  torque  -  prestress  relation  is  available,  which  reads  as  [16] 

Torque  =  0.2  x  (Prestress  on  ceramics)  x  (A r„  >.  of 
ceramics)  x  (Outer  diameter  of  threaded 
rod) 

an  actual  calibration  curve  was  obtained  by  placing  a  strain  gage  on 
the  steel  sleeve  of  the  P-7  transducer.  This  curve  is  shown  in  Fig. 

5-33.  It  is  seen  that  the  empirical  relationship,  which  is  a  straight 
line,  gives  a  prestress  which  is  slightly  too  high  at  large  torque 
values.  With  a  torque  at  125  ft -lb  (typical  for  a  P-7  transducer),  the 
experimental  prestress  given  by  Fig.  5-33  is  5.4  ksi,  while  the  empirical 
one  predicts  7-5  ksi. 

The  experimental  results  are  shown  in  Figs.  5-34 >  5-35  and.  5-36. 

For  comparison,  the  theoretical  results  are  also  included  in  the  above 
figures,  being  given  as  the  dashed  lines.  In  the  following  figures, 
"theory  // 1"  refers  to  one -dimensional  rod  theory;  while  "theory  #2" 
refers  to  the  spring  model.  It  is  seen  that  the  one -dimensional  rod 
theory  generally  gives  a  resonant  frequency  which  is  too  high  compared 
to  that  given  by  the  experiment.  Furthermore,  the  variation  of  the 
resonant  frequency  with  the  prestress  is  much  less  than  the  experimental 
results.  On  the  other  hand,  the  spring  model  gives  very  good  results 
for  the  P-7  transducer,  which  can  be  seen  from  Fig.  5-36.  Further,  at 
high  prestress  in  the  vicinity  of  the  operation  point  (5.4  ksi  for  P-7)» 
the  theoretical  results  and  the  experimental  results  are  quite  close. 

From  these  results,  one  may  conclude  that  the  change  of  the  material 
properties  due  to  prestress  effects  could  not  account  for  all  the  fre¬ 
quency  shift  observed  experimentally.  However,  at  high  prestress  level, 
material  properties  appear  to  account  for  most  of  the  change.  Since  a 
high  prestress  level  tends  to  increase  the  interface  friction  which 
may,  in  turn,  suppress  the  vibrational  modes  other  than  the  longitudinal 
ones  of  the  ceramic  rings,  the  ceramic  may  behave  more  like  a  longitudinal 
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torque  for  P-7 


Fig.  5-33  Calibration  and  empirical  curves  for  the  torque  versus  prestress 
relation  of  a  P-7  transducer 


Experimental 
Theory  number  I 
Theory  number  2 


rod.  At  low  prestress  level,  there  must  be  other  effects.  It  could  be 
the  friction  effect  at  the  interfaces  or  the  two-dimensional  effect  of 
the  ceramics. 


of  course,  both  the  theoretical  approaches  au  e  far  from  perfect. 
This  is  evident  frcm  Figs.  5-34,  5-35  and  5-36.  In  Fig.  5-3*+}  theory 
#1  is  off  in  both  magnitude  and  shape,  though  the  variation  of  the 
resonant  frequency  with  the  prestress  at  high  prestresu  level  is  com¬ 
parable  to  that  of  the  experimental  ones.  In  Fig.  5-35}  it  is  seen 
that  above  moderate  prestress  level  theory  # 1  agrees,  better  than 
theory  #2,  with  the  experiment  ad.  values.  In  Fig.  5-3t,  theory  #2  is 
in  accord  with  the  experimental  results,  while  theory  #1  is  again  off 
in  both  magnitude  and  shape.  Unfortunately,  there  seems  to  be  no  con¬ 
sistency  in  both  approaches.  However,  at  ceramic  prestress  levels  of 
approximately  5  ksi  or  greater,  the  discrepancies  seem  to  be  less 
severe  for  both  the  theories . 

The  temperature  effect  on  resonant  frequency- -Some  data  were 
obtained  previously  in  the  OSU  Sonic  Power  Laboratory  on  this  effect 
[17].  Unfortunate!.'/,  the  data  are  somewhat  rough  since  the  temperature 
was  not  even  along  the  transducer.  Therefore  it  is  difficult  to  make 
any  meaningful  analysis  on  these  data.  To  overcame  this  difficulty, 
the  transducer  (or  resonator)  was  placed  in  an  oven  and  slowly  hee.ted 
to  a  given  temperature  (130°F),  at  which  time  the  oven  was  turned  off. 
While  the  transducer  (or  resonator)  was  gradually  cooling  down,  resonant 
frequencies  were  measured  at  various  temperature  readings.  The  results 
are  given  in  Figs.  5-37,  5-38  and  5-39  as  solid  lines.  On  these 
figures,  the  theoretical  results,  as  given  previously  in  Figs.  5-22, 

5-24  and  5-25}  are  also  plotted  as  dashed  lines  for  more  direct  compar¬ 
isons.  Note  that  the  slopes  of  the  theoretical  curves  are  similar  to 
those  of  the  experimental  ones  which  would  imply  that  the  decrease  in 
the  theoretical  resonant  frequency  due  to  th  thermal  prestress  release 
effect  is  comparable  to  that  of  the  experimental  observation.  This 
may,  in  turn,  imply  that  the  main  cause  of  the  decrease  in  resonant 
frequency  due  to  temperature  increase  is  the  thermal  prestress  release 
in  the  ceramics.  However,  it  was  observed  that  when  the  temperature 
of  the  transducer  returned  to  room  temperature  after  the  test,  the 
resonant  frequency  was  higher  than  the  original  one  (before  the  test) 
by  a  finite  amount.  Evidently  the  heat  treatment  causes  some  permanent 
changes  in  the  PZT-4  properties.  Krueger  [12]  noted  that  this  effect 
can  be  used  to  stabilize  the  ceramic  properties.  The  heat-stabilizing 
treatment  consists  simply  of  immersing  the  specimens  in  200°C  oil  for 
one  hour  according  to  Krueger. 

The  voltage  effect  on  the  resonant  frequency--It  has  been  observed 
in  the  OSU  Sonic  Power  Laboratory  that  a  voltage  increase  will  cause  a 
resonant  frequency  drop  for  a  transducer  [18],  The  experimental  results 
in  this  investigation  reaffirm  the  above  observation.  Figures  5-40, 

5-41  and  5-42  give  data  for  the  symmetric  simple  resonator,  the  P-7 
resonator,  and  the  P-7  transducer.  They  also  reveal  that  the  paths  of 
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Fig.  5-37  Resonant  frequency  versus  temperature  for  the  symmetric 
simple  resonator 
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Fig.  5-40  Experimental  resonant  frequency  versus  voltage  for  the 
symmetric  simple  resonator 
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Fig.  5-4l  Experimental  resonant  frequency  versus  voltage  for  the 
P-7  resonator 
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Fig.  5-42  Experimental  resonant  frequency  versus  voltage 
for  the  P-7  transducer 


increasing  and  decreasing  voltages  do  not  coincide.  In  general,  the 
path  for  a  decreasing  voltage  is  lower  than  that  for  an  increasing 
voltage.  This  is  evidence  that  a  high  voltage  on  the  ceramic  causes 
some  irreversible  changes  of  the  ceramic  properties  as  well  as  certain 
reversible  effects. 


5-4  SUMMARY 

The  ultimate  objective  of  this  work  has  been  to  model  the  transducer, 
predict  resonant  frequency  and  to  account  for  frequency  shifts  due  to 
prestress,  temperature  and  voltage.  To  accomplish  the  goal,  the 
symmetric  simple  resonator,  the  P-7  resonator  and  the  P-7  transducer 
have  been  formulated  by  one -dimensional  wave  theory,  and  the  resulting 
frequency  equations  were  solved  numerically  on  a  computer.  Prediction 
of  the  resonant  frequency  shifts  has  been  attempted  from  two  stand¬ 
points:  l)  by  utilizing  the  one-dimensional  rod  model  to  consider  the 
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effects  of  piezoelectric  material  property  changes  due  to  environmental 
changes  and  2)  by  modeling  the  ceramic  as  a  massless  spring  and  leaving 
the  spring  constant  to  be  determined  by  resonance  experiments  on  the 
symmetric  simple  resonator  to  account  for  the  prestress  effects  quasi- 
theoretically .  Finally,  resonance  experiments  which  included  the 
prestress,  the  temperature  and  the  voltage  effects  were  conducted.  The 
experimental  results  were  compared  to  the  theoretical  ones  where  they 
were  available. 

Comparison  between  the  experimental  and  the  theoretical  results 
shows  that  at  low  prestress  level  on  the  ceramics,  the  one-dimensional 
longitudinal  rod  model  gives  resonant  frequencies  much  too  high;  while 
the  spring  model  has  a  better  agreement.  At  high  prestress  level  near 
the  operation  point  (5.4  ksi  on  ceramics),  both  the  one -dimensional  rod 
and  the  spring  model  give  results  close  to  the  experimental  ones. 
Although  both  theoretical  models  predict  an  increase  of  resonant  fre¬ 
quency  for  an  increase  in  prestress  level  in  accord  with  the  experimental 
observations,  neither  of  the  theoretical  models  accurately  predict  the 
rate  of  increase  in  resonant  frequencies .  The  temperature  and  voltage 
effects  on  the  transducer  (resonator)  are  not  thoroughly  treated  because 
of  the  lack  of  material  property  change  data  in  the  literature.  For  a 
more  complete  assessment  of  those  effects,  more  experiments  must  be 
performed.  Nevertheless,  the  prestress  release  effect  due  to  the 
thermal  expansion  of  the  transducer  elements  was  found  to  be  in  accord 
with  the  experimental  observations.  However,  the  finite  jump  of  the 
resonant  frequency  after  heat  treatment  can  only  be  accounted  for  after 
data  of  the  experimental  material  property  changes  are  available. 
Recycling  information  on  heat  treatment  and  voltage  exposure  also  needs 
to  be  investigated  in  the  future . 

Discrepancies  between  the  rates  of  increase  of  the  resonant 
frequencies  as  a  result  of  the  prestress  may  suggest  that  the  one¬ 
dimensional  rod  modjl  is  not  adequate.  A  more  exact  mathematical  model 
may  be  needed.  The  numerical  result  of  the  one -dimensional  rod  model 
reveals  that  a  30$  change  of  the  Young's  modulus  can  only  shift  the 
resonant  frequency  of  a  simple  resonator  by  about  200-300  Hz  (comparing 
to  about  1400  Hz  obtained  experimentally).  Possibly,  modeling  two- 
dimensional  effects  and  including  the  interface  friction  in  the  formu¬ 
lation  will  improve  the  theoretical  results. 

Nevertheless,  in  conclusion,  some  interesting  results  should  be 
mentioned.  First,  it  was  noted  from  the  symmetric  simple  resonator 
analysis  that  the  elastic  constants  of  a  piezoelectric  material  are  the 
dominating  factors  in  resonant  frequency  calculation.  Secondly,  the 
stress  release  due  to  thermal  expansions  does  contribute  to  frequency 
shift  for  a  temperature  rise.  This  effect  can  be  reduced  by  using  a 
longer  center  bolt  or  soft  (with  smaller  Young's  modulus)  metals  to 
prestress  the  ceramics  and  thus  to  reduce  the  spring  constant  or,  this 
can  also  be  achieved  by  employing  metals  with  thermal  expansion  coeffi¬ 
cients  near  that  of  the  ceramic. 
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6.  DEVELOPMENT  OF  A  RJLSE  REFLECTOR  FOR 
SailC  TRANSDUCERS 


The  ultrasonic  transducer  is  a  device  which  converts  electric  input 
energy  into  vibratory  mechanical  energy.  It  consists  of  piezoelectric 
ceramic  rings  sandwiched  in  a  resonant  horn  assembly.  The  horn,  which 
is  an  energy  concentrator,  is  in  turn  connected  to  a  transmission  line 
that  transmits  mechanical  energy  to  the  working  surface.  The  piezo¬ 
electric  ceramic  rings,  which  are  the  most  important  elements  in  the 
transducer,  have  two  interesting  properties:  (i)  For  a  given  elec¬ 
trical  disturbance  to  the  piezoelectric  ceramic,  mechanical  pulses  are 
produced,  and  (ii)  given  a  mechanical  disturbance  to  the  ceramics, 
an  electrical  pulse  is  produced.  Due  to  the  first  property  of  the 
piezoelectric  ceramic,  the  sonic  transducer  has  been  widely  used  for 
converting  electric  energy  into  mechanical  energy  and  thus  has  found 

application  as  a  mechanical  tool  for  drilling,  cutting,  and  welding.  A 
technique  used  for  the  transmission  of  energy  from  the  sonic  transmission 
line  to  the  work,  surface  is  impact  coupling,  where  a  small,  bouncing 
slug  of  metal  is  placed  between  the  end  of  the  vibrating  sonic  trans¬ 
mission  line  and  the  work  surface.  This  slug  of  metal,  usually  called 
the  "tool"  and  typically  a  cylinder  one  inch  in  diameter  and  about  an 
inch  long,  bounces  back  and  forth  between  line  tip  and  work  piece 
through  repeated  impacts  w?th  and  rebounds  frem  the  respective  surfaces. 

As  the  tool  rebounds  from  the  working  surface  and  consequently 
impacts  on  the  tip  of  the  transmission  line,  a  stress  pulse  is  generated 
and  propagated  along  the  transmission  line  into  the  transducer.  As  a 
result  of  the  second  property  of  the  piezoelectric  ceramic,  a  high 
voltage  electric  pulse  is  then  produced  across  the  face  of  the  ceramic 
ring.  The  high  voltage  electric  pulses,  especially  whenever  coupled 
with  the  input  electric  voltage  from  the  electric  power  supply,  may 
cause  an  arc-over  phenomena  across  the  face  of  the  ceramic.  This  elec¬ 
tric  arc  may  crack  the  piezoelectric  ceramic  or  eventually  cause 
dielectric  breakdown.  Therefore,  the  development  of  a  reflector  device 
to  prevent  stress  pulses  frem  entering  the  transducer  becomes  desirable. 

The  basic  concept  of  designing  a  pulse  reflector  originates  from 
understanding  that  a  stress  pulse  propagating  along  a  transmission  line 
will  be  completely  reflected  from  a  fixed  boundary  condition  and  will  be 
largely  reflected  by  a  sudden,  large  increase  in  cross-section.  There¬ 
fore,  a  desirable  configuration  for  a  pulse  reflector  would  be  a  stepped 
cylinder  which  closely  simulates  the  fixed  edge  boundary  condition. 

Early  studies  in  the  P-7  transducer  by  Feng  [l]*  and  Libby  [2]  observed 
that  P-7  transducers  did  not  fail  during  the  impact -coupling  applica¬ 
tion,  but  it  was  noted  that  this  was  not  the  case  for  the  P-11  trans¬ 
ducer.  This  observation  suggested  that  the  step  at  the  P-7  transducer 
may  reflect  a  large  amount  of  the  stress  pulse  and  thus  prevent  the 
failure  of  the  piezoelectric  ceramic. 


♦Numbers  in  brackets  refer  to  references  at  the  end  of  this  chapter. 
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,  P1  or^eJ  to  design  a  stress  pulse  reflector  without  affecting  the 

JSe  transducer -transmission  line  system,  the  pulse  reflector 
should  have  the  same  resonant  frequency  as  the  transducer  so  that  after 
the  reflector  is  connected  to  the  transducer,  the  whole  system  will 
vibrate  at  the  same  resonant  frequency. 

Based  on  previous  preliminary  studies,  two  basic  types  of  pulse 
reflectors  were  developed.  Theoretical  analyses  to  investigate  the 
resonance  and  reflection  characteristics  of  the  pulse  reflectors  were 
conducted.  The  theoretical  results  were  later  verified  by  experimental 
tests.  It  was  found  that  with  the  pulse  reflector  connected  to  the 

transducer,  80$  of  the  incident  pulse  was  reflected.  The  voltage 
generated  across  the  face  of  the  piezoelectric  ceramic  rings  by  the 
pulse  was  also  reduced  80$. 

6-1  REFLECTOR  DESIGN 

(A)  Theoretical  Analysis  of  Reflection 
Characteristics 

As  the  first  step  in  the  design  of  a  stress  pjlse  reflector, 

1  is  necessary  to  understand  the  reflection-transmission  characteristic 
of  a  stress  pulse  at  a  discontinuity.  A  one -dimensional  wave  theory  is 

used  for  this  purpose.  The  configuration  of  the  discontinuity  is  shown 
in  Fig.  6-1. 


Reflected  wove 


Tronsmitted  wave 


Incident  wave 


— ►; 


Fig.  6-1  Waves  across  a  discontinuity 
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Cross-sections  and  material  densities  for  section  I  and  II 
are  Ax,  A2,  and  pj_ ,  p2,  respectively.  An  arbitrary  incident  pulse 
fi(x  -  Clt),  may  be  considered  to  be  made  of  many  frequency  components 
in  a  Fourier  sense.  Thus,  a  superposition  of  many  waves  of  the  type 
uj  =  A  exp[i7i(x  -  c1t)]  may  be  used  to  obtain  a  pulse  shape  f(x  -  cit), 
where  uj  and  cj.  are  the  displacement  and  wave  velocity  of  the  incident 
wave.  tte  r*  flee  fed  wuve  5n  section  T  and  the  transmitted 

wave  in  section  II  can  be  expressed  as  up  =  gs(x  +  Cit),  and 
up  =  f2(x  -  c2t),  respectively. 

Boundary  conditions  at  the  discontinuity  are  based  on  conti¬ 
nuity  of  force  and  of  particle  velocity  and  are 


Aj,  Cj  +  ctr)  —  A20j  ,  x  —  0 


VI  +  VR  =  VI 


x  =  0 


(6-1) 


where 


Vj  =  -  Cj.fi1,  VR  =  Cl  g2  ’ ,  vT  =  -  c2f2 1 

^Ut  ^Ur  ^Urp 

ai  =  Sl  IbT  =  Eifl  ’»  °r  =  Ei  “ST  =  Elg2' »  ai  =  E2  "5T  =  E2f2' 

(6-2) 

and  where  Ex  and  E2  are  the  Young's  modulus  of  section  I  and  section  II, 
respectively.  The  prime  indicates  differentiation  with  respect  to  the 
argument . 

Substituting  eq.  (6-2)  into  eq.  (6-1)  and  solving  for  f2'  and 
g 2*,  yields 

«  »  _ SCjAjEj _  p  i 

fs  =  AiEiC2  +  A2E2Ci  fl 


e  »  =  1A2E2  ~  c3AiE1  f  , 

AiEic2  +  A2E2ci 


(6-3) 


Since  ctt  =  E2f 2 1 ,  aR  =  E1g2’  and  Ex  =  cf^,  E2  =  cft;2,  the  results  in 
eq.  (6-3)  may  be  put  in  the  form 


Ai£ici  +  A2£2c2  I 


tR  Ai^Ci  +  A2£2c2  I 


(6-4) 


Equation  (6-4)  represents  the  relationships  between  transmitted  stress 
Op,  reflected  stress  ctr  and  incident  stress  Oj. 


If  the  materials  for  both  sections  are  the  same,  i.e., 

El  =  e2  =  E,  U  =  £2  =  5,  and  cL  =  c2  =  c,  eq.  (6-4)  can  be  expressed  as 


2Ai 

91*  =  A,  +  A^  CTI  = 


x 


^2.-1 

Ap  -  A-l  At  _ 

=  AfTAt  ai =  frs  1 

Ai 


(6-5) 


Plots  of  the  relations  between  2/(1  +  A2/Al ) ,  (A2/Al  -  l)/(l  +  A2/A1)  8113 
A2/Ax  axe  presented  in  Fig.  6-2.  For  A2/Ax  =  9  it  is  seen  that  80f0  of 
the  incident  stress  should  be  reflected  at  the  discontinuity  and  onl^ 

2C$  of  the  incident  stress  should  be  transmitted  to  section  I. 

(B)  Analysis  of  Resonance  Characteristics 
of  Reflectors 

Two  types  of  pulse  reflectors  axe  considered  in  the  present 
investigation.  The  "Type  A"  reflector  is  shown  in  Fig.  6-3- 
one-dimensional  analysis  for  the  stepped  cylinder,  simple  strength  o 
materials  theory  and  Hooke's  law  are  assumed  to  apply.  The  governing 
differential  equations  of  motion  are 


(6-6) 


where  uL  and  u2  axe  displacements  in  section  I  and  II,  respectively  and 
c2  =  E/^.  Solutions  for  eq.  (6-6)  are: 


Uj.  =  (Bj.  sin  kxj.  +  B2  cos  kxx)  cos  ut 
u2  =  (B3  sin  kx2  +  B4  cos  kxx)  cos  ut 
where  k2  =  u2/c2  and  u  is  the  radial  frequency. 

Boundary  conditions  for  the  stepped  cylinder  are. 


(6-7) 


Ux  =  -u2  , 


at  Xx  =  0,  x2  =  0 
at  Xx  =  0,  x2  =  0 


at  r.x  =  Li  =  5  in. 


at  x2  =  L2  =  5  in. 


(6-8) 


172 


Fig.  6-3  Dimensions  and  coordinates  of  Type  A 
pulse  reflector 


Substituting  solutions  of  eq.  (6-7)  into  the  boundary  conditions  of 
eq.  (6-0),  the  following  set  of  equations  is  obtained; 

B4  =  0,  -  A2B3  =  0 

Bxk  cos  kLx  -  B2k  sin  kLx  =  0 

B3k  cos  kL2  -  B*k  sin  kL2  =  0  (< 

These  can  be  further  reduced  to 


B3  cos  kLi  -  B4  sin  kLx  =  0 


B3  cos  kLx  -  B4  sin  kLi  =  0  (6. 

The  only  nontrivial  solutions  to  eq.  (6-10)  are  those  for 
which  the  determinant  of  coefficients  is  equal  to  zero.  This  deter- 
min ant  is  given  by 


!  xip 

cos  kLi  sin  kLx 


cos  kL2 


-sin  kL; 


r 


which  expands  to  the  frequency  equation 


Ap 

5 *  cos  kLx  sin  kL2  +  sin  kLx  cos  kL2  =  0  (6-12) 

For  the  present  case  A2/Ax  =  9,  and  I*  =  L2  =  5  in. ,  the 
frequency  equation  becomes 

10  cos  kLi  sin  kLi  =  0  (6-13) 

The  solutions  for  eq.  (6-13)  are 

kLi  =  n(rt/2),  (n  =  1,2,3,...) 

Let  n  =  1  for  the  first  mode  of  vibration,  to  give 


«  _1 _ it 

“  2  ’  Li  "  10 


so  that  the  first  mode  resonance  is 


kc  _  10 
=  2k  =  2rt 


“  =  9989-4  Hz 


A  "Type  B"  reflector  configuration,  shown  in  Fig.  6-4,  was 
also  considered.  The  basic  idea  was  to  have  a  half-wavelength  device, 
with  the  discontinuity  in  the  center.  The  method  of  analysis  closely 
followed  that  used  for  the  Type  A  reflector.  Thus,  the  governing 
equations  of  motion  are 


yuT  _  _i_  yuj 

cixi2  c2  dt2  ’ 

dfua  -  1  d^Up 
6x3^  c1-  ~5t 


^^u  _  1  yu2 

dxs2  = 


(6-14) 


where  ux ,  u2,  and  u3  are  displacements  in  sections  I,  II,  and  III, 
respectively  and  c2  =  e/£  . 


Solutions  for  eq.  (6-l4)  are 


ux  =  (A  sin  kxx  -u  B  cos  kxx  )  cos  ut 

u2  -  (C  sin  kx2  +  D  cos  kx2)  cos  ut 

u3  =  (F  sin  kx3  +  G  cos  kx3)  cos  ut  (6-15) 

where  k  =  co  /c  ,  and  u  is  the  radial  frequency.  The  boundary  conditions 


Fig.  6-4  Dimensions  and  coordinates  of  the  Type  B 
pulse  reflector 
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Substituting  the  solutions  of  eq.  (6-15)  into  the  boundary 
conditions  of  eq.  (6-l6),  yields  the  following  set  of  equations; 

A  cos  kL-.  -  B  sin  kLx  =  0 


I 

I 


\ 

mm 


*» 


r 


i 

i 

i 

i 

i 

i 

i 


C  sin  kL2  +  D  cos  kL2  -  G  =  0 
A2(C  cos  kL2  -  D  sin  kL2)  -  A3F  =  0 
F  cos  kL3  -  G  sin  kL3  =  0 


(6-17) 


This  reduces  to 


Ai 


C  cos  kLj,  +  D  sin  kLi  -  0 

C  sin  kL2  +  D  cos  kL2  -  G  =  0 

C  cos  kL2  -  D  sin  kL2  -  ^  F  =  0 


F  cos  kL3  -  G  sin  kL3  =  0 
This  determinant  of  coefficients  is  given  by 


(6-18 


£2  cos  kl,  sin  kLx  0 

Ai 

sin  kL2  cos  kL2  0 

cos  kl2  -sin  kL2  -  A3/A2 

0  0  cos  kL3  -sin  kL3 

which  gives  the  frequency  equation 


0 

-1 

0 


=  0  (6-19) 


^2.  cos  k In/ —  cos  kL2  sin  kL3  +  sin  kLx  cos  kL3^ 

Ax  ^■\A2 

-  sin  kLx(^  sin  kL2  sin  kL3  -  cos  kL2  cos  kL3)  =  0  (6-20) 

For  the  present  case  A2/Ai  =  9,  A3/A2  =  l/9»  ~  ^3  ’ 

and  L2  =  3  in.,  so  that  eq.  (6-20)  becomes 

cos  3k  cos  3-5^  cos  3-5k  +  9  cos  3-5k  sin  3k  cos  3-5k 

-  *  sin  3.5k  sin  3k  sin  3-5k  +  sin  3-5*  cos  3k  cos  3-5k  =  0 

9  (6-21) 

The  solution  for  the  first  mode  is  k  =  0.48433,  giving  a 
resonance  frequency  of 

kc  0.48433  c  _  =  15,323  Hz 

2n  2n 
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6-2  MEASUREMENT  OF  RESONANCE  AND  REFLECTION 

CHARACTERISTICS 

Experimental  measurements  of  the  resonance  and  reflection  charac¬ 
teristics  is  desirable,  to  verify  some  of  the  analysis,  particularly 
since  there  are  sane  deviations  of  the  configuration  used  in  analysis 
from  the  actual  pulse  reflector.  Three  types  of  tests  were  conducted. 
The  first  test  was  the  measurement  of  the  resonance  characteristics  of 
the  Type  A  and  B  pulse  reflectors.  Also,  resonance  tests  were  made  on 
a  P-11  as  a  check  of  that  design.  The  second  test  was  measurement  of 
the  stress  reflection  from  the  Type  A  and  B  reflectors.  Additionally, 
pulse  reflection  tests  were  also  done  on  the  P-11  horn.  The  third 
test  was  a  comparison  of  the  reflection  characteristics  of  a  P-11 
transducer  with  and  without  the  Type  A  pulse  reflector. 

(A)  Measurement  of  Resonance  Characteristics 

The  basic  set-up  for  this  test  consisted  of  two  small  magnetic 
transducers,  a  Bruel  and  Kjear  1013  oscillator,  an  electronic  counter 
and  an  oscilloscope  as  shown  schematically  in  Fig.  6-5* 


B( BrUel  and  Kjear 
MM  -  000?.) 


Fig.  6-5  Schematic  diagram  of  apparatus 
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The  B  and  K  oscillator  generated  a  harmonic  signal  to  the  mag¬ 
netic  transducer  A  and  it,  in  turn,  vibrated  the  pulse  reflector  through 
magnetic  coupling.  The  magnetic  transducer  A  was  placed  approximately 
one-sixteenth  of  an  inch  from  the  end  face  of  the  pulse  reflector.  An 
identical  magnetic  transducer  was  used  as  a  pickup  at  B,  with  its  output 
signal  displayed  on  the  oscilloscope.  By  varying  the  driver  frequency 
frctn  the  oscillator,  the  resonance  frequency  of  the  pulse  reflector  was 
easily  determined  from  the  maximum  amplitude  of  the  vibration  shown  on 
the  oscilloscope.  The  results  of  the  tests  are  shown  in  Table  6-1. 


Table  6-1  -  Resonance  Frequencies  of  Reflectors  and  P-11  Horn 


Specimen 


Type  A  reflector 
Type  B  reflector 


P-11  horn 


Resonance  Frequency  Observed 


10,779  kHz 
14,309  kHz 
10,308  kHz 


(B)  Measurement  of  Reflection  Characteristics 

Incident  and  reflected  stress  pulses  along  a  transmission  line 
can  be  monitored  by  means  of  strain  gages  attached  on  the  transmission 
line.  Jhe  incident  stress  wave  is  produced  by  firing  a  one  inch  in 
diameter,  one-inch-long  steel  tool  from  a  spring  gun  as  shown  in  Fig. 
6-6.  By  comparing  the  amplitude  of  the  incident  stress  pulse  and  the 
reflected  stress  pulse,  one  can  determine  the  percentage  of  the  incident 
wave  being  reflected  from  the  pulse  reflector.  The  results  can  be  used 
directly  for  comparison  with  the  previous  theoretical  results. 


Oscilloscope 


Signal 

processor 


Strain  gages 


[tn*\ 


if] 


Test  of  Type  A  pulse  reflector — The  results  of  the  Type  A  pulse 
reflection  are  shewn  in  Fig.  o-7*  It  is  seen  that  the  joint  at  section 
where  the  reflector  attaches  to  the  transmissi ;n  line,  and  ex , 
causes  some  slight  degree  of  influence  in  reflection  of  the  stress 
pulse.  It  is  also  observed  that  the  amplitude  of  the  reflected  wave 
is  approximately  80 %  of  the  incident  wave;  these  results  are  quite 
close  to  theoretical  prediction. 

Test  of  Type  B  pulse  reflector — The  results  for  the  Type  B  reflec¬ 
tor,  shown  in  Fig.  b-o,  are  similar  to  the  case  of  the  Type  A  reflector; 
an  effect  due  to  the  connection  between  the  transmission  line  and  pulse 
reflector  is  rioted.  Again,  the  amplitude  of  the  reflected  wave  equals 
to  80%  of  incident  wave,  and  is  in  good  agreement  with  the  theoretical 
calculations . 

Test  of  P-11  horn--It  is  interesting  to  observe  the  reflection 
characteristics  of  a  P-11  horn.  As  shown  in  Fig.  6-9>  the  amplitude 
of  the  reflected  wave  equals  approximately  50%  of  the  incident  wave 
amplitude .  This  means  approximately  50%  of  the  energy  from  the  incident 
pulse  will  transmit  into  the  transducer. 

(C)  Measurement  of  Pulse  Reflector  Efficiency 

Tests  were  made  to  determine  the  difference  in  electrical 
pulse  amplitude  from  the  ceramic  rings  of  a  P-11  transducer  with  and 
without  the  Type  A  reflectors  in  place.  Figure  6-10  shows  the  setup 
and  results  fer  the  P-11  transducer  without  the  reflector.  The  top 
trace  of  the  photograph  is  the  familiar  stress  wave,  while  the  lower 
trace  is  the  voltage  generated  by  the  stress  wave  transmitted  into  the 
transducer. 

Figure  6-11  shows  the  results  with  the  Type  A  reflector  in 
place.  It  is  seen  that  the  induced  voltage  pulse  is  reduced  from  $.b 
kV  without  the  reflector  to  a  level  of  1.0  kV  'with  the  reflector. 


6-3  SUMMARY 

The  theoretical  analysis  of  the  reilection  characteristics  indicated 
that  a  pulse  reflector  design  with  a  stepped  cylinder  configuration  and 
an  kp/h-i  ratio  equal  to  9  would  provide  80%  reflection.  Both  theoretical 
calculations  and  experimental  tests  verified  these  findings. 

Of  the  two  types  of  pulse  reflectors  considered,  it  was  found  that 
the  resonance  frequency  for  the  Type  A  pulse  reflector  was  much  closer 
to  the  natural  frequency  of  the  P-11  transducer  than  the  Type  B  reflec¬ 
tor,  so  that  only  the  Type  A  reflector  was  considered  in  final  tests 
with  the  P-11.  Final  comparison  tests  indicated  that  the  voltage 
generated  at  the  ceiamic  rings  by  the  tool  impact  was  reduced  more 
than  80%  compared  with  the  results  without  the  pulse  reflector. 
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Vert  div.  =  0.5  V 
Horiz.  div.  =  50  nsec 


*►  Reflected  wave  from  section  ex 
-Reflected  wave  from  section  dx 
— Reflected  wave  from  section  cx 

•Reflected  wave  from  section  bx 
Incident  pulse  at  section  ax 
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Fig.  6-7  Typical  stress  pulse  output  observed  for  Type  A 
reflector  with  interpretation  of  the  wave  form 
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Fig.  6-8  Typical  stress  pulse  output  for  Type  B  reflector 
with  interpretation  of  the  waveform 


Horiz.  div.  =  50  nsec 


L*— ►  Reflected  wave  from  f3 
Reflected  wave  from  e3 
‘—►Reflected  wave  fram  d3 
'—►-Reflected  wave  from  c3 
—►Reflected  wave  from  b3 

oecond  impact  wave  due  to  misalignment  between  tool  &  transmission  line 
—^..Incident  pulse  at  section  a3 


Fig.  6-9  Stress  pulse  output  for  P-11  horn  with  interpre¬ 
tation  of  the  waveform 
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Typical  stress  pulse  and  voltage  output  for  I -11  trans 
ducer  with  the  Type  A  p-  t:e  reflector  in  pl^  .e 
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